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PRÓLOGO 


La presente monografía titulada Dos 
Temas de Mecánica de la Construcción re- 
coge la traducción al español de dos mo- 
nografías publicadas en el Instituto de 


Ingeniería de la Construcción de Moscú 


titulada Cálculo de Cáscaras Prismáúticas, 


de V.N. Pastuchijin y el Método de los 
Elementos Finitos en las Tareas de la 
Teorta de las Estructuras, de N.N. Ieon- 
tiev e 1.1. Diomin en los años 1974 y 


1979 respectivamente. 


Estas dos publicaciones mantienen toda 


su vigencia y constituyen temas de actua- 


lidad para los ingenieros civiles cuba- 
nos dedicados al análisis y proyecto de 


las estructuras. 


En la obra de V.N.Pastuchijin se pre- 
sentan las kases del método de V.Z. Vla- 
sowv para el cálculo de estructuras ple- 
gadas de paredes delgadas como folded 


plate, bóvedas cilíndricas, y otras. 


El método de análisis utilizado es muy 


familiar a todo ingeniero civil, pues se 


parte de los esquemas de análisis planos 


tradicionales, y de la aplicación del mé- 


todo de los desplazamientos en su forma 


clásica. Sin embargo, la esencia del mé- 


todo consiste en que el análisis estruc- 
tural para obtener el gráfico de M, V, N 


de los desplazamientos se descompone en 
la superposición de dos problemas o eta- 


pas de trabajo: 


A En el primero de ellos (luego de estar 
establecido el sistema base del método 
de los desplazamientos) se resuelve el 
gráfico de momentos flectores de una 
estructura de nudos indesplazables 


bajo la acción de las cargas actuantes. 


H En la segunda etapa se obtienen todos 
los gráficos de momentos flectores que 
corresponden a los casos de un asenta- 
miento unitario de todas, y cada una 
de las ligaduras lineales introducidas 
para convertir la estructura en un caso 
de nudos indesplazables, y se multipli- 
can después estos gráficos por el des- 
plazamiento que el sistema de cargas 
exteriores produciría sobre la estruc- 
tura, en la dirección de cada una de 
estas ligaduras teniendo en cuenta la 
naturaleza de los vínculos en la di- 
rección perpendicular al plano de pór- 
tico, y las condiciones de borde de 
los elementos de paredes delgadas en 
dicha dirección. 


De esta forma se pasa a un análisis 


del comportamiento espacial de la estruc- 


tura, párticularizádo en este caso para 
estructuras de dimensiones en plantas 
con relaciones e/b > 5, pero cuyo proce= 
dimiento se puede fácilmente extender a 
OtrOs Cásos. 

En la obra de N.N. leontiev a J.l. 
Diomin se presentan los conceptos fiuda- 
mentales del método de cálculo más uni- 
versal y extendido hoy en día en el 
mundo. 

El método de los elementos finitos se 
aplica a algunas tareas Simples de tras= 
misión del calor, teoría de la elastici- 
dad y solución de sistemas de ecuaciones 
diferenciales ordinarias. 


Á Su vez, te presentan los principios 
fundamentales de la programación por el 
método de los elementos finitos, los cua- 
les han alcanzado un gran desarrollo 
desde los años anteriores a la publica- 
ción de N.N. Leontiev hasta nuestros 
días. El conocimiento de estos prinei- 
pios constituye una necesidad ineludible 
para todo ingeniero que pretenda mante- 


nerse medianamente informado de las posi- 
bilidades y los métodos de cálculo exis= 
tentes para el análisis de las estructu- 


ras en la actualidad. 


Ambas monografías fueron elaboradas 
como libros de textos para estudiantes de 
la Facultad de Ingeniería Civil que estu= 
diaban en aquella época la especialidad 


de Teoría de las Estructuras. 


Al mismo tiempo se consideraba la ne- 
cesidad de utilizarlas en estudiantes de 
las especializeciones de Estrueturas mu 
tálicas, ue Proyecto y Construcción de 


Obras Civiles e Industriales. 


Las tradujimos al español pensando que 
este material será de sumo interés para 
la mayoría de los ingenieros civiles de- 
dicvados al análisis y al proyecto de las 
estructuras, al ampliar sus posibilidades 
de trabajo con el conécimiento de métodos 
de actualidad que no formaron parte de 
nuestros planes de estudios en aquellos 


años. 


DrSc. Angel E, Castañeda llevia 


El MÉTODO DE LOS ELEMENTOS FINITOS 


TEMA 1 


EN LAS APLICACIONES A LA TEORÍA 


DE LAS ESTRUCTURAS 


INTRODUCCIÓN 


Para realizar la selección de la mejor 
solución constructiva de los elementos de 
las obras de ingeniería es necesario 
tener la pesibilidad de pronosticar el 
comportamiento de los campos físicos en 
dependencia del cambio en las condiciones 


de existencia de las construcciones, 


Para la determinación de un cuadro cua- 
litativo y cuantitativo de los campos £Í- 
sicos se utilizan métodos analíticos y mé- 


todos numéricos de cálculo. 


La ventaja fundamental de los métodos 
analíticos consiste en que los mismos per- 
miten obtener una forma general de la so- 
lución, la cual determina el comportamien- 
to de la característica buscada (por ejem- 
plo: temperatura, desplazamientos) en 
cualquier punto de la región y en cual- 


quier momento dado. 


Esa propiedad de continuidad de la so- 
lución en la región dada o en el interva- 
lo de tiempo analizado permite reflejar 
la influencia de todos los factores que 
actúan sobre la construcción, realizar un 
análisis comparativo de los mismos y se- 
leccionar los fundamentales entre todos 


ellos. 


La deficiencia fundamental de estos 


métodos consiste en que la forma general 


de la solución reviste, en muchos casos, 
una forma tan voluminosa que la misma no 
puede llevarse hasta un resultado numéri- 
co concreto, sin la ayuda de métodos nu- 


méricos especiales, 


Esa dificultad en la mayoría de los 
casos líquida las ventajas señaladas de 
los métodos analíticos, pues los programas 
de tabulación de las fórmulas generales de 
solución no son menos complejos que los 
programas realizados por métodos numéri- 
cos, pero su región de aplicación es, ade- 
más, mucho más estrecha que la de estos 


últimos. 


Además, los métodos analíticos son apli- 
cables, por regla general, a estructuras 


de forma constructiva simple, y son total- 
mente no aplicables para el cálculo de 


regiones no homogéneas de geometría arbi- 
traria, como es el caso de las construc- 


ciones hidrotécnicas. 


Los métodos numéricos permiten encon- 
trar el valor de la función buscada en un 
número finito de puntos de la superficie 
de cálculo; En estos casos, no se obtie- 
ne una forma general de la solución, lo 
cual constituye la desventaja fundamental 


de estos métodos. 


Sin embargo, brinda la, posibilidad de 
analizar regiones, de forma arbitraria, 


bajo una distribución arbitraria de las 


'cargas y de las propiedades físicas del 
material, lo cual compensa por completo 


la desventaja anteriormente señalada. 


Además, el desarrollo intensivo de las 
técnicas de computación de las últimas 
décadas, ha llevado a una considerable ex- 
tensión de los métodos numéricos en los 
cálculos de ingeniería y ha desplazado 


los métodos analíticos en muchos casos, 


En el presente, los métodos numéricos 
recibieron una amplia aplicación en la 
mayoría de las organizaciones de proyecto 


y de investigación científica del país. 


Esto ha planteado la necesidad de es- 
tudiar los métodos más difundidos en los 
cursos de preparación de los especialis- 
tas de la rama de la construcción pon los 


centros de educación superior. 


El presente libro de texto está dedi- 
cado a uno de los métodos numéricos más 
universales dentro de la teoría de las 
estructuras: el método de los elementos 
finitos (MEF) que constituye el logro más 
significativo de la mecánica de la cons- 
trucción en el período posterior a la 


Gran Guerra Patria. 


El método ha influido de forma signi- 
ficativa en el desarrollo posterior de 
esta rama de la ciencia, la cual, en nues- 
tra época, ha variado sustancialmente su 


contenido. 


Actualmente, no es solo la rama de la 
ciencia dedicada a la construcción e in- 
vestigación de los modelos matemáticos 
sobre el comportamiento de los cuerpos de- 
formables, sino también la rama de la 
ciencia dedicada a la elaboración de com- 
plejos computacionales de diversa funcio- 
nalidad para la solución de diferentes 
tareas de la teoría de las estructuras, 
la cual se construye sobre la base del mé- 
todo de los elementos finitos. 


Además, y conjuntamente con todo ello, 
se produce el surgimiento de la mecánica 
de la construcción en regiones del conoci- 
miento que hasta este momento le eran ajé- 
nos, como son la física de transmisión del 
calor, la teoría del flujo potencial, la 
aerohidrodinémica y otras. 


Los conceptos tradicionales de la mecá- 
nica de la construcción como esquema de 
cálculo (o esquema de análisis), matriz 
rigidez, sistema base y otros, se encuen- 
tran ahora en las más diversas publicacio- 
nes sokre la aplicación del método de 10% 
elementos finitos en otras ramas de la 


ciencia. 


La idea de discretizar el medio conti- 
nuo para la simplificación de los cálcu- 
los técnicos surgió ya en el pasado 


siglo XIX, 


Sin embargo, los primeros intentos 
para realizarlo en una forma sistemática 
ligados a la representación de cierta re- 
gión en forma de un conjunto de elementos 
discretos como se señala en el artículo 
de Zienkiewicz, o.c pertenecen a 
A. Jrenikov (1941), D. Max Genri (1943) y 
N. Niumark (1949)', los cuales aproxima- 
ron el comportamiento de un medio conti- 
nuo a través de un conjunto de elementos 


y barras. 


Posteriormente, este procedimiento re- 
cibió el nombre de método de las aproxi- 
macion=s con barras (método de los ele- 
mentos discretos) con ayuda del cual fue- 
ron Fesueltes un amplio grupo de tareas 


aplicadas. 


La descripción aproximada de un cierto 
medio elástico con ayuda de las uniones 
en los nudos donde se interceptan los ele- 
mentos discretos, recibió el nombre poste- 


riormente de método de los elementos fini- 


tos, lo cual fue propuesto por M, Turner” 
Ry Clark” y De Argiris'., 

En esos trabajos se formuló por prime- 
ra vez el método de los elementos finitos 
en una forma muy parecida a como nosotros 


lo conocemos actualmente y se determina- 


ron sus propiedades fundamentales. 


Después del año 1958, el método de los 
elementos finitos se desarrolló de forma 
intensiva y ya en el año 1964 determinó 
prácticamente su aparato matemático en 


Forma de un método variacional. 


la cantidad total de publicaciones 
sobre el método de los elementos finitos 
Y sus aplicaciones se puede calcular en 
unos cuantos miles, 


51657 Ya 
222% on diferentes 


En las monografías, 

. e R ; 
resúmenes y artículos se da una descrip- 
ción general de las bases del método de 
los elementos finitos y de los ejemplos 


de su utilización en diferentes ramas de 


la ciéncia y la técnica. 


En nuestro país el método de los ele- 
mentos finitos comenzó a extenderse desde 
el ano 1966. El motivo de la demora en 
la aplicación del mismo lo provocó la de- 
mora en el perfeccionamiento de las máqui- 
as computadoras de producción nacional, 
sin las cuales no era posible la realiza- 
ción de los algoritmos de numerosas eta- 
pas que requiere el método de los elemen- 


tos finitos. 


Desde el ano 1967, un trabajo activo 
para el desarrollo y la introducción del 
método de los elementos finitos se lleva 
a cabo por el VNIIG (Instituto de Inves- 
tigación Nacional de la URSS), B.F. Vene- 
diev. 

Desde el comienzo de la década del 70, 
el método de los elementos finitos recibe 
una amplia y multifacética extensión y 


aplicación, la cual permite, ya en nues- 


tros días, constatar la existencia de una 
escuela soviética del método de los ele- 
mentos finitos, la cual se distingue por 
una serie de rasgos característicos a la 
vora de tratar diversos aspectos del mé- 


todo, 


En los trabajos de los científicos so- 


Dl Es9310 E ON 
viéticos >? se áa úna descripción total 
de las bases rel rétodo de los elementos 


finitos para los proll=o:as de la teoría 


jo en el estalleciniento ¿e los planes de 
estudio en diversas disciplinas. 
1 EOxtengio GN E iz la presente 


publicación no permite arbúrcar todas las 
Yamas de la ciencia en las cuales se apli- 
ca el mátodo de los elenontos finitos, 
losotros expondreros las operaciones 
fundamentales del método en tareas linea- 
les (y en el plano) de transmisión del 
calor y de teoría de la elasticidad, ade- 
más, de detenernos en los problemas de la 


programación del método. 


Después del estudio del libro de texto, 
el estudiante se encuentra en posibilida- 
des de orientarse dentro de la literatura 
especializada con relación al método de 
los elementos finitos y resolver algunas 
tareas simples en las máquinas computa- 


doras'', 


1. CONCEPTOS BÁSICOS 


Todos los fenómenos físicos transcu- 
rren en el espacio y en el tiempo. La 
investigación analítica de los mismos nos 
conduce al estudio de los cambios que se 
producen en una cierta magnitud fundamen- 


tal ligada con dicho fenómeno. 


El conjunto de los valores de la mag= 
nitud fundamental en todos los puntos del 
espacio estudiado se conoce con el nombre 


de campo de dicha magnitud. 


Así por ejemplo, durante el análisis 
del proceso de la transmisión del calor 
analizamos un campo de temperaturas, mien- 
tras que en las tareas de teoría de la 
elasticidad analizaremos un campo de des- 


plazamientos, deformaciones y tensiones. 


Por Cuanto, en cada punto del campo de 
temperaturas se puede establecer un solo 
valor de la misma, el campo de temperatu- 
ras se llama escalar, por su parte, los 
campos de desplazamientos y de tensiones 
son campos vectoriales por cuanto a cada 
punto corresponde no sola el valor de 2 
característica correspondiente, sino, ade- 


más, una cierta dirección en el espacio. 


Se llaman campos no estacionarios uue- 
llos campos cuyas característica: no son 
solamente función de las coordenadas del 
punto en el campo, sino además, del tiem- 
po. Si por el contrario, las caracterís- 
ticas no dependen del tiempo, el campo se 


denomina como estacionario. 


Los puntos del campo que se caracteri- 
zan por tener un mismo valor de la carac- 
terística buscada y que forman líneas 
continuas se llaman isolíneas. En el 
caso de las líneas de igual temperatura, 
estas líneas reciben el nombre de iso- 


termas . 


En el caso general, cualquier fenómeno 
físico posee un infinito número de parti- 
cularidades de carácter geométrico, del 
material desde el punto de vista de las 
influencias externas y de otros factores, 
las cuales resulta imposible tener en 
cuenta en su totalidad durante el análi- 


sis del mismo y, además, no es necesario. 


Debido a ello, el fenómeno en estudio, 
deberá ser sustituido por un cierto mode- 
lo físico, el cual posee un número finito 
de particularidades, las cuales influyen 
de una forma decisiva y determinante 
sobre las características, cuyo estudio 


interesa al investigador. 


De esta forma, para la selección del 
modelo físico de una armadura real on la 
mecánica de la construcción se utiliza 
la hipótesis de la unión articulada de 
las Larras de la armidina en los puntos 

. > . - * 
de unión en lugar de Zas uniones rígidas 


de la armadura real. Así opbteneros un 


cómodo modelo físico lo lo que en este 
caso es le mismo, 41 esquema de cálculo 
ámodo) para la determinación de lor us- 
ruerzos interiores cu las arras. 

El cálculo de las armaduras teniendo 
en cuenta las uniones r1Íqidas en los 
nudos de la armadura complica de forma 
considerable el cálculo, tiene una peque- 
ña influencia sobre el valor de las fuer- 


zas normales en las harras de la armadura, 
pero constituye también un modelo físico 

* . e 
que representa una aproximación más pre- 


cisa al trabajo real de la armadura, 


Una mayor complicación del modelo 
puede lograrse teniendo en cuenta los es- 
fuerzos durante el montaje, el carácter 
no ideal de las ligaduras, las influen- 
cias por temperatura, por corrosión y 
otros aspectos. Todas estas magnitudes 
pueden tener una pequeña influencia sobre 
la magnitud de los esfuerzos, pero compli- 
can de forma significativa el modelo fí- 


sico. 


La selección del modelo físico depende 
considerablemente de la intuición del in- 
vestigador, la cual a su vez es resultado 
del conocimiento y la comprensión del ob- 


jeto de investigación y de sus caracte-. 


riclicar, depende conjuntamente con ello 
de su habilidad para seleccionar los fac- 
tores fundamentales, que actúan sobre el 
objeto estudiado y de la potencialidad de 


los métodos de cálculo disponibles. 


En calidad de modelo físico” de un 
medio real, el cual se adoptará en esta 
monografía, tomemos la hipótesis del 
medio continuo, es decir, ignoramos la 


estructura molecular de la sustancia. 


Supondremos, además, que todas las 
cargas e influencias exteriores son apli- 
cadas estáticamente, y qué componentes 
separados del medio, son en sí mismos ho- 

e . - . 
mogéneos e isótropos, y el comportamien- 
to físico de los materiales que los com- 


ponen está sometido a la lay de llooke. 


Por modelo matemático de la tarea de 
investigación entendemos el conjunto de 
correlaciones, que describen los procesos 
que ocurren en el modelo físico, y de las 
ecuaciones que describen sus condiciones 
de borde. Cada modelo físico puede co- 
rresponder á una cierta cantidad de mode- 


los matemáticos. 


Así por ejemplo, analizando el compor- 
tamiento de una viga sobre una base elás- 
tica, nosotros podemos obtener diferentes 
modelos matemáticos en dependencia del 
tipo de base adoptada: si está compuesta 
por un medio de Wincler, si es un semies- 
pacio elástico, si está formado por un 


medio de una sola capa y otros. 


En cada uno de esos casos, el sistema 
de ecuaciones que describe la tarea, cam- 
bia su forma, a pesar de que el modelo 


físico se conserva.» 


Analicemos dos formas fundamentales de 
obtención de modelos matemáticos. Plan- 
teémoros la tarea de describir la línea 


elástica de una viga de una sola luz y de 


sección transversal constante (fig. 1.1). 


—>X 
EI 


Figura 1,1 


Del curso de Resistencia de Materiales 


sabemos «(que : 


Ely" =M Gb1) 


O» 


teniendo en cuenta que: 


¡EY = UU 0 (120 
donde; 


b: tiódulo de Young. 
iz iomento de inercia de la sección 
transversal de la viga, 


Refiriéndonos al ceso de bordes liga- 
dos en la forma mostrada en la figura 1.1 


podemos escribir que: 


ti 
$“ 


y (0) =y* (0) = 0 


(143) 


II 
ú 


y (L) y" (2) 0 


eS 


De esta forma, la ecuación (1.2) y las 
condiciones de borde (1.3) determinan en 
forma única y por completo el modelo ma- 
temático (o dicho en otra forma, el plan- 
teamiento matemático de la tarea investi- 


gada) . 


La formulación rigurosa de esta tarea 


sería la siguiente: 


Encontrar la función de la línea elás- 
tica de la viga y(u)p) la cual debe ser 
determinada y continua en el intervalo 
cerrado 0 «zx <“l conjuntamente con 

todas sus derivadas (hasta la derivada 


cuarta incluida) la cual satisface den- 

tro del intervalo abierto 0<x<la 

la ecuación diferencial (1.2) y a las 

condiciones de borde (1.3) en los ex- 

tremos del intervalo.' 

La forma presentada de conformación del 
modelo matemático la llamaremos modo ope- 


racional... 


Plantéemonos esta misma tarea desde un 
punto de vista energético. El trabajo de 
los momentos flectores en la viga se de- 


termina por la expresión: 


L 


A 3 
. M du A » 1) e IS 
== 4 El" ae = 


de Mi 


a sra (1.4) 


ll tratajo de las fuerzas exteriores: 


a iS CA 
es = J qydx = - q J ydx (1.5) 


La energía potencial de deformación de 
toda la viga se determina por la expre- 
... 
sión: 


Ply(2)] =4y +4, =S 14" -qylds 
0 


(1.6) 


Por la conocida propiedad de la ener- 
yÍía potencial, la cual tiene en la lite- 
ratura el nombre de Principio de Lagran- 
ye, la verdadera ecuación de la línea 
elástica de la Larra será aquella, que le 
brinda un valor mínimo a la expresión de 
la energía potencial [en este caso a la 


ecuación (1.6)]. 


De esta forma la tarea planteada por 
este otro modo de resolución la podremos 


tormular en la siguiente forma: 


Encontrar la función de la línea elás- 
tica de la viga y (uc), la cual debe ser 
determinada y continua en el intervalo 
cerrado 0 <x <l conjuntamente con su 
primera y segunda derivada, aplicando 

en los extremos del intervalo los va- 

lores dados en (1,3) y la cual le pro- 


porciona a la energía potencial de de- 
formación (1.6) un valor mínimo. 
Esta forma de construcción del modelo 
matemático lo llamaremos modo variacio- 
nal. Destaquemos que, en este caso, se 
exige la satisfacción de las condiciones 
de continuidad hasta la segunda derivada 
de la función y no hasta la cuarta como 


en el modo operacional. 


Eso en muchos casos simplifica la 

4 .. - - , pa 
aplicación de metodos numéricos, Lasados 
en la minimización aproximada de expre- 


siones del tipo de la !órmula (1.6). 


En un planteamiento más riguroso de 
la tarea variacional es necesario tuner 
en cuenta su estreclo vínculo con el mé- 
todo operacional y exiyir la satisfacción 
de las condiciones de continuidad de de- 


rivadas de orden superior. 


El enlace de la forma variacional y 
la forma operacional del modelo matemáti- 
co se puede ¡¡Lantear en la siguiente 


forma: 


Se tiene en el caso general una expre- 


sión de la forma: 


Al da dz 
Vatesy)) =5 S FlosYs ds 
D y 
y? A 
I : 0" z DA oz dll 
Ou 0w0y ay 


(1.7) 


La función 2(%,4) proporcionándole un 
valor extremo a la expresión (1.7) deberá 


. 1 . 
satisfacer según'? la ecuación: 


aF 9 28 
TETAS A A 
9% (F.) a? (7,) 
E t 
Lo AA AKÁKÁ = 0 1.8 
0x0 y 5 Ena 


ao 


a? rs 
t=- Z 5 Po O (1, = P,qsPs8,€) 
dy? + 


La ecuación (1.8) no es otra cosa cue 
la formulación operacional de la tarea 
variacional. Comprobemos esta afirmación 


para el caso (1.6). 


simb5licamente la ecuación (1.6) se 
escribe en la forma: 
| Sel E 
Ply(m)1 = F(, Ly; 
0 da? 

2 ,. 
y Ly yy - ay 

da* 


AO e 
D(y”) 


Entonces en el lugar de la ecua- 


ción (1.8) tendremos: 


2 
oF A E, 
dy da? 


o teniendo en cuenta que: 


E 
dy q 
y: 
(7) 
2 AR > 
da* 
obtenemos : 
1V 


Ely” -q=0 


o lo que es lo mismo, se obtuvo la ecua- 
ción operacional (1.2) que era lo que nos 


proponíamos demostrar. 


De esta forma es evidente que las for- 
mas variacionales y Operacionales del mo- 
delo matemático están estrechamente liga- 
gas entre sí y hablando de una manera 
rigurosa, a la función y(w) en amos cax 
minos es necesario aplicar las mismas 


exigencias, 


Sin embargo, si se analizan ambos 


modos por separado uno del otro (como co- 
múnmente se hace), se puele ntilizar una 
forma menos exigente de la: condiciones 
de continuidad en el caso viriacional. 


La expresión concreta de la función 
v (1) puede obtenerse tanto «die (1.2) como 
de (1.6) por diferentes métodos. Por 1mé- 
todo de solución llamaremos al conjunto 
de procedimientos con ayuda de los cuales 
a partir del modelo matemático adoptado 
se pueden determinar las funciones de las 


características buscadas. 


Por solución de la tarea entenderemos 
el proceso de tránsito desde modelo mate- 
mático adoptado hasta los valores numéri- 
cos concretos de las características bkus- 


cadas. 


Es necesario señalar que a cada modelo 
matemático pueden corresponder varios mé- 


todos de solución. 


De esta forma, por ejeiplo, la ecua- 
ción (1.2) puede resolverse por integra- 
ción directa, por el método de las dife- 
rencias finitas, por el método de las 


colocaciones y OtroSa 


En el caso (1.6) se puede utilizar el 
método de Ritz, el método de los elemen- 
tos finitos, el método de las diferencias 


de variaciones y otros. 


En las aplicaciones de la teoría de 
las estructuras la forma más extendida 
corresponde al modo operacional de plan- 


teamiento de las tareas. 


Sin embargo, para la aplicación del 
método de los elementos finitos la forma 
variacional es las más cómoda y ventajosa, 
aunque el método puede ser aplicado a la 


forma operacional (ver el epígrafe 8). 


En lo que sigue utilizaremos el plan- 
teamiento variacional del modelo matemá- 


tiCÓs 


Le MODELO MATEMÁTICO DE LA TAREA 

ESTACIONARIA DE LA TRANSMISIÓN 

DE CALOR 

Con el objeto de hacer más cómoda y 
evidente la explicación de los postulados 
teóricos fundamentales del método de ele- 
mentos finitos analizaremos los mismos en 
el ejemplo de la tarea estacionaria de la 


transmisión del calor, 


En este caso, a Cada punto del medio 
se encuentra ligada una sola variable 
desconocida: la temperatura en dicho 
punto. Además, todos los postulados del 
método, en este caso, tienen un contenido 
físico preciso y una representación evi- 


dente y clara. 


Para medios isótropos compuestos por 
un solo material la tarea estacionaria de 
la transmisión del calor se formula de la 
siguiente forma; encontrar la función de y 
la temperatura 0 (x,y) la cual debe ser 
determinada y continua en la región ce- 
rrada D, conjuntamente con sus derivadas 
parciales de primer y segundo orden, y 
que satisfaga dentro de esta región a la 


ecuación diferencial: 


AÑO + Q=0 (2.1) 


y en la frontera S de la región D a la 


condición de primer género: 


0(x,yWl = p(,y) (2.2) 
Ss 


donde : 


A: Coeficiente de transmisión del 


calor del material, 


Y: Operador de Laplace que es 
y? y? 
Et_—, 

da? a 


igual a 


Q (1,4): Potencia de las fuentes distri- 


buidas de calor. 
Para simplificar la prosentac“ión del 
método no analizaremos las condiciones de 


segundo y tercer gyónerd. Nos limitaremos 


a encontrar la teuporatura en el interior 
de la región en el caso ue conocer La 
distribución de la tempezatura en la fron- 
tera. 

En el cazo en que la región está com- 


puesta de varios componentes, la ecua- 

ción (2.1) se camlia por un sistema de 

ecuaciones, que posee el mismo número de 
ecuaciones que de materiales componentes, 
y las condiciones de borde (2.2) se com- 
plementan con las condiciones de contac- 
to en la frontera MM que separa las capas 


adyacentes (1,/J): 


90 . (4) 00 . (41) 
9, (01) =0 .(M) 5 A —Í— 22) y 
¿ 97 yy on ,, 
(2.3) 


donde: 


Ty? Normal común a la frontera que se- 
para los componentes en el punto de la 


frontera M. 


El planteamiento variacional correspon- 


diente a esta tarea sería el siguiente! 


Encontrar la función de la temperatu- 
ra 0 (x,y) la cual debe ser determinada y 


continua en la región cerrada DP conjunta- 
mente con. sus derivadas parciales de pri- 
mer orden y que adopta en la frontera 5 
de la región los valores dados (2.2) y 

le proporciona un valor mínimo a la inte- 


gral: 


A 00 
IO 


D 


Pas rl )=40 Jdady 


(2,4) 


En el caso unidimensional del calenta- 
miento de una barra, la integral (2.4) 
adopta la forma: 


90 
J (0) ap? [£ > 
o L E 


e 0.0] dx (3.5) 
donde : 


Ll: Longitud de la barra, 


3. POSTULADOS FUNDAMENTALES DEL MÉTODO 
DE LOS ELEMENTOS FINITOS 


Analicemos las etapas fundamentales de 
solución de una tarea por el método de 


los elementos finitos. 


1. La región dada o la barra se divide 
por líneas imaginarias o puntos en una 
cierta cantidad de elementos fini- 
tos (Liga 3 1Me 


Figura 3.1 
a) Tarea plana; b) Tarea lineal. 


La división de la región dada en ele- 
mentos separados constituye una de las 
operaciones más delicadas y de mayor res- 
ponsabilidad del método, del desarrollo 
de la cual en mucho depende la precisión 
y la economía de la solución de toda la 


tarea. 


Como es evidente de la figura (3.1), 
la disposición de los elementos en la 
región, puede ser en el caso general arbi- 
traria, lo cual es muy ventajoso y rinda 
su rasgo distintivo del desarrollo de la 
solución por el método de los elementos 
finitos con relación a los esquemas rigu- 
rosamente regulares del método de las di- 
ferencias finitas. Además, en una misma 
región es posible la aplicación de diver- 


sos tipos de elenentos. 


La cantidad de elementos en una tarea 
concreta se determina, generalmente, por 
el grado de detalle deseado y por la poten- 


- 


cialidad del programa utilizado. 


Como regla general en las regiones de 
mayor interés para el investigador se 
utiliza una malla muy densa, y en otras 


regiones una malla menos densa. 


En el caso general es muy difícil dar 
cualquier tipo de recomendación yeneral 
sobre la designación de la malla debido 
a que la solución de este asunto depende 
en mucho de la forma de la región dada y 
de la intuición del investigador. 

2. Los elementos se consideran unidos 


entre sí en un número finito de puntos 
o nudos, 


En el método de los elementos finitos 
se inserta una hipótesis (figura 3.1) 
relacionda con el hecho de que en caso 
de coincidir las características busca- 
das (temperatura, desplazamientos) en 
los nudos de ynión de los distintos ele- 


mentos, coinciden a su vez las caracte- 


rísticas (temperatura, desplazamientos) 
de todos los puntos, dispuestos en las 
líneas que unen dichos nudos, (línea de 


contacto de los elementos). 


3. Para cada elemento finito se establece 
una. función, que de forma unívoca de- 
termina la característica buscada (tem- 
peratura; desplazamientos) en el inte- 
rior del elemento a través de los 
valores de las características en los 
nudos » 


Aclaremos este aspecto con un ejemplo, 
tomemos una karra recta (figura 3.2) de 


longitua (/) con dos nudos en los ex- 


tremosS . 
e- OS XXX XXX XX XX) Xx 
O 
| 2 
L 
HA A>—]—]á- 
Figura 3.2 


Se exige encontrar la función que per- 
mita determinar la temperatura de todos 
los puntos de la barra a través de la tem- 


peratura 0, y 0, en sus extrenos. 


Para poder determinar esta función es 
imprescindible, evidentemente, seleccio- 
nar una cierta ley de variación de la ten 
peratura a través de la longitud del ele- 
mento, Supongamos que la temperatura 


varía de acuerdo con una ley lineal. 


0 (1) = a, + 4,2 (3:01) 


donde: 
A, y 4, : Parámetros indeterminados. 


Para determinar la magnitud de los pa- 
rámetros debemos tener en cuenta las con- 
diciones de borde, esto es, que para 


<c=x,+0(x,) =0,, y parax=<,; 
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0 (x,) = 0. 


mos el siguiente sistema de ecuaciones: 


A partir de lo cual obtene- 


0, = q, +0,%, 


Resolviendo este sistema con relación 


a los parámetros 4, y 4, tendremos que: 


2 


1 i 


(13,3) 


y entonces la expresión (3.1) puede es- 
cribirse en la forma: 


0 (a) = + (9 0d] 


1 ; 
1014 = 0 
(3.4) 

la cual después de una serie de transfor- 


maciones elernentales se puede escribir 


CONO : 
0 (a) = —Á—[-0, (6209) +0, (0 -2,)] = 
Lalo ss L 1 Le da q A = 
= 8,0, + 4,0, =[M0 13.5) 
donde : 
> 
N, == AA 
A 
Nh = 5 ; 
1] = [0 , Nal; 
Y 
E 9, 
0 = 
0 


De esta forma, la temperatura de un 
punto cualquiera de la barra puede deter- 
minarse a partir de la temperatura en sus 
nudos dentro'de. los límites de la ley de 


distribución tomada como aproxima- 

ción (3.1) con ayuda de la fórmula (3.5). 
Las Cinciónes [M tienen el nombre espe- 
cial de funciones de la forma del ele- 


mento, 


El tipo de estas funciones se determi- 
na de antemano de acuerdo con la ley aúop- 
tada de variación de la característica 


buscada a través del elemento. 


Esas funciones son invariantes en re- 
lación con el sentido físico de la carac- 
terística buscada (de esta forma, la fun- 
ción Ú puede ser, por ejemplo, el vector 
de los desplazamientos en lugar de la 
temperatura) y dependen solamente, Como 
es evidente a partix de las expresiones 
para NC = 1,2),de los parámetros geomé- 


tricos del elemento finito seleccionado. 


El sentido físico de las funciones /'; 


es el siguiente (figura 3.3). 


a) 


Figura 3.3 


La función de la forma es la función 
de la distribución pesada de la caracte- 
rística buscada a través del elemento 
finito. Ella vale la unidad en el nudo, 
en el cual coincide el número de la fun- 
ción con el número del nudo y vale cero 


en los restantes nudos. 


La necesidad de establecer de antemano 
la función de distribución de la caracte- 
rística buscada a través del elemento nos 


lleva al hecho de establecer como solu- 


ción de la tarea a investigar el compor- 
tamiento de una determinada ley dentro de 
los límites del elemento finito (como, 


por ejemplo, la dada por la fórmula (3.1). 


En relación con esto dentro de los 
límites del elemento, la solución verdade- 
ra puede diferenciarse fuertemente de la 


solución numérica | fig. 3.4 aj]. 


En estos casos, es necesario tener en 
cuenta, que en los nudos del método de 
los elementos finitos, como regla general, 
se obtienen valores de las característi- 
cas que se diferencian muy poco de los 


valores exactos. 


Si este mismo intervalo se divide en 
una serie de invervalos más pequeños el 
error de aproximación, o lo que es lo mis- 
mo, la interpolación lineal entre los va- 
lores de la característica en los nudos 
se disminuye considerablemente, como se 
indica en la figura 3.4 b). A partir de 
esto, se comprende la necesidad de hacer 
más densa la red de elementos finitos en 
aquellos lugares en que se espera puedan 


producirse fuertes variaciones de los va- 


lores de la característica buscada. 


Solución del M.E. F 


a) 


La TE 
a De 
.olución 


verdadera 


Solución del M.E.F 


b) He 


92 


solución 
ver dadera 


El método de los elementos finitos 
tiene, además, otra posibilidad muy im- 
portante para elevar la precisión de la 
solución, la cual consiste en utilizar 
elementos que posean un grado mayor de 
aproximación, así para un elemento con 
tres nudos [figura 3.5 a)] se puede esta- 
blecer la siguiente ley de variación de 


la temperatura. 


0 (2) =a, +42 +0, (3.6) 
a) 
Ye 
pe a aca 
Solución del M.E.F 
b) Su Solución 
verdadera; 


Figura 3,5 


En este caso, la interpolación lineal 
a través de dos puntos del elemento 
[figura 3.4 a)] se cambia por una inter- 
polación cuadrática.a través de tres pun- 
tos del elemento, lo que permite para 
las mismas dimensiones obtener una preci- 
sión considerablemente superior en los 


valores de la característica buscada. 


Además de elevar la precisión, la uti- 
lización de elementos de un grado mayor 
(lo que quiere decir con un mayor número 
de nudos) permite utilizar en una misma 
tarea, elementos de mayores dimensiones 


[ compare las figuras 3.4 b) y 3.3 b)]. 


Ello permite disminuir el número total 


de los elementos utilizados y lo que es lo 
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mismo, disminuye el tiempo total que re- 
quiere la máquina computadora para obtener 


la solución de la tarea. 


Es evidente a partir de esto la posibi- 
lidad de utilizar elementos no solo de se- 
gundo grado, sino también de tercero, cuar- 


to y aún de grados superiores. 


Sin embargo, en estos casos, la comple- 
jidad del cálculo de las funciones de la 
forma aumenta, pues se hace necesario re- 
solver sistemas de ecuaciones de grados 


superiores. 


Al lector se le recomienda obtener de 
forma independiente la matriz indicada 
[MN =10,, MM, N,] para el elemento de 


segundo grado. 


4. Todos los factores de fuerzas interio- 
res y exteriores (potencialidad de las 
fuentes de calor distribuidas, carga 
exterior, tensiones, etc.) actuantes 
sobre el elemento se llevan a la forma 
de esfuerzos en los nudos. Todos los 
factores de fuerzas (tensiones, defor- 
maciones, cantidad de calor) del inte- 
rior del elemento se expresan a través 
de las características fundaméntales 
(desplazamientos, temperaturas, etc.) 
de sus nudos. 


Vayamos nuevamente al estudio de la 
tarea unidimensional de la transmisión 
del calor. Aceptando que la temperatura 
del interior del elemento dado se deter- 


mina por una función del tipo (3.5); 


0 (x) = [1]0, (3340 


'y poniendo la función en el integral (2.5) 


por transformaciones sucesivas obtenemos 


para el elemento que: 


y 2 3(M0,) d bi 
IA pe E 0 di= 
1,0) Tf AA da Y AMO de 
+y / y / e 
= e 0 ( 414) ue de a - 
2 E vu PU da € 
, l y YES 
- a? | . | 19 + 0 P 
E € € e “€ 


donde: 


YU . 
Lo, 98N]" 90 > 
¿A as 


Matriz de ri- 


gidez del ele- 
mento. 
= E ALO 
La =- J Q|M' du : Vector de carga o 
de los términos in- 


dependientes. 


De acuerdo con lo visto anteriormente 
el signo 7 indica transposición de la 


matriz. 


Durante la deducción se utilizaron las 
siguientes reglas del análisis matricial: 
Pedo area ia? =ta tel 
cuadrado de la matriz fila [c]. Además, 
por cuanto el vector de la temperatura en 
los nudos no depende de x, el mismo se 


puede sacar fuera del signo de integral. 


El vector real de los esfuerzos en los 
nudos del elemento debe satisfacer la con- 
dición de darle un valor de mínimo a la 
función (3.8), o lo que es lo mismo, la 


candición: 


ÓN (0) (3.9) 


3. 10) 


pa 
S 
= 
+ 
ra 
8 
o 


Aquí se utilizan las reglas 


e E 00 
de A 20 
as A) ad 


De esta forma, los esfuerzos nodales 
E, se expresan a través de las temperatu- 
ras nodales A con ayuda de la matriz ri- 


gidez del elemento [4]. 


Para el caso lineal (ver figura 3.2) 
la expresión de la matriz rigidez del ele- 
mento finito (MRE) se determina de la si- 


guiente forma: 


ESTAS 


e ENE ES (3.11) 
be qee 


Senalemos algunas propiedades de la 
matriz rigidez del elemento (MRE). El 
orden de la matriz se determina por el nú- 
mero de grados de likertad del elemento, 
esto quiere decir, por la dimensión del 
vector de las características nodales 
(características en los nudos). A pesar 
de que en las tareas de transmisión del 
calor cada nudo se encuentra ligado solo 
a un valor de la temperatura, el elemen- 
to presentado en la figura 3.2 tiene dos 


grados de libertad, el de la figura 


3.5 a) y b) tiene tres y otro elemento 


podría tener aún mayor número de grados 
de libertad. 


La matriz rigidez (3.11) es simetrica 
con respecto a la diagonal principal, lo 
que es una consecuencia de la ley de re- 
ciprocidad de la naturaleza, de la cual 
el teorema de Betti (en la mecánica de la 
construcción) constituye un caso parti- 


cular. 


La MRE caracteriza las propiedades fÍ- 
sicas del medio dentro de los límites del 


elemento finito. 


Es necesario señalar el carácter singu- 
lar de la MRE, pues el determinante dde la 
mi'sma es igual a cero. Esto es una conse- 
cuencia del hecho de que la exigencia de 
un valor mínimo de la integral de la ener- 
gía es la condición de equilibrio del 


sistema físico. 


En la mecánica de la construcción este 
hecho recibe el nombre de principio de La- 
grange para la energía potencial de De- 
formación. 

5. Se determina la matriz suma de la ri- 
gidez del sistema (MRS) de elementos, 
que expresa las fuerzas nodales á tra- 
vés de las características nodales. 
Después de la total determinación de 

las propiedades de un elemento aislado 

pasaremos al estudio de todo el sistema. 

Analizaremos sobre un ejemplo los princi- 


pios fundamentales de esta traslación. 


Se tiene una barra con una longitud de 
7 m y el extremo izquierdo de la misma se 
calienta a una temperatura de 10 E y el 
extremo derecho se enfría a una tempera- 
tura de (-4C). El coeficiente de trans- 


misión de calor del material es: 


BT 
Me 


Dividamos la barra en siete elementos 
finitos de primer orden (figura 3.6) y nu- 
meremos los nudos y los elementos. Es 
necesario señalar que se podría haber usa- 
do una malla no regular, la cual no hubie- 
ra complicado en nada el análisis que 


vamos a realizar, 


USE END E z 
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La matriz rigidez de todos los elemen- 


tos en las condiciones A/l = 1, Tendría 
la forma: 

1 -1 0, 
Ll, = , 

E EN 

1 1 al 
(A, = Po 

=y 9 0 

1 -1 0, 
[A], = 


O escrito en forma simbólica: 


A 


d 


lod 1 
4% 1 h, 2| | 
| 


la), 


ñ 


LAl, 


1 
i 
y 
1 
2 4 a 
a? h y 
22 23 ¿ 
..o., 
0 


[A1, = 


Es necesario prestar una atención espe- 
cial sobre el hecho de que los Índices de 
los elementos de las matrices rigidez de 
los elementos (MRE) en las últimas fórma- 
las se numeran en rigurosa corresponden- 
cia con los números de los nudos. 
Analicemos ahora el proceso de tránsi- 
to de las MRE separadas, cada una de 13 
cuales está ligada solo con el vector d, 
hacia la matriz que la liga con el vector 


5 T 
resultante 0 = 0, 9,0, -0,) 


0 10) [0] 
11 12 e 
A 0 0 
21 2 212 33 
, e? +? yn 0 
32 33 187 SA 
e HR en yA 
0 (0) 43 44 314. 45 
0 0 0 al (ne + »/ 
54 ss 
o 0 0 h* 
65 
10) 0 0 0 0 
0 0 0 0 0 


Arialicemos por ejemplo el tercer nudo: 
El mismo es común para los elementos 2 y 
3. Entremos por las filas de las corres- 
pondientes matrices de rigidez LAl, y LAál, 


en la matriz común (2.12), 


Por cuanto el elemento ft, ¿ Pertenece 
al mismo tiempo a ambas matrices de rigi- 
dez de los elementos, en la célula corres- 


pondiente de la matriz se coloca la suma 


De una forma similar se colocan todos 
los elementos de las MRE en la matriz ge- 
neral (3.12). En las células no comple- 
tadas se colocan ceros debido a que las 
.Correspondientes temperaturas nodales no 
se incluyen en los vectores o. del ele- 


mento. 


Así por ejemplo, el punto (4) no se 
incluye en el vector 0, para el primero 
CA 
y el segundo elemento, consecuentemente, 
el correspondiente valor del elemento de 
la fila en la matriz rigidez del sistema 


debe ser igual a cero, 


De igual forma se determinan los com- 


ponentes del vector F. Por ejemplo, 


0 0 0, 

0 0 o, 

0 0 Y 

= 0 

0 0 o, 
+ %h / 0 0 

Pride ie 

6 7 7 

h +2 h 0 

( 77 3%) 18 + 
E Ar 7) 

87 al ' B 

3 (3,12) 


para el nudo 3, F, = 
1 


Py + Fyj para el 


nudo 8, £, =PF,. 
donde: 
A - 
Po Valor de los vectores Fo para cada 
u 


uno de los elementos. 


De forma definitiva, el vector '' para 
el caso de la ecuación (3.12) tiene la 


forma : 


pos 1 
E E 
1 1 de de , 3 4 4 5 


PS A 
5 6 7 7 3 


En forma simbólica la ecuación (3.12) 
tiene la forma: 
(3. 13) 


[2] A 


donde: 


M: Número de elementos. 


La ecuación (3.13) de forma semejante 
a la ecuación (3.10) establece la condi- 
ción de mínimo de la energía total del 
sistema, la cual en el caso general posee 


la forma: 


15 (3.14) 


J0) LF 19 d+ G 

Las propiedades de la matriz rigidez 
del sistema (MRS) son las mismas que las 
de las matrices rigidez de los elementos 


(MRE) . 


Además, prestemos atención al hecho de 
que el ancho de banda en la ecuación 
(3.12) es finito, en las condiciones de 
una numeración consecutiva de los nudos 
obteniéndose una matriz en banda, de 


ancho en banda igual a tres. 


La propiedad de la banda de la matriz 
rigidez del sistema (MRS) para cualquier 
tipo de elementos finitos, en condiciones 
de una correcta numeración de los nudos 
de la matriz se utiliza ampliaente para 
la confección de programas económicos de 
solución de sistemas de ecuaciones del 
tipo (1,139. . 

Regresemos a nuestro ejemplo (figura 
3.5) y supongamos que dentro de los lími- 
tes de la barra Q = 0. Introduciendo en 
la ecuación 3.12 los valores numéricos ; 
correspondientes, tendremos de forma de- 


finitiva que: 


e 0 o! o) 

1. 2 - S 5. ma, 
par Pp 20 o | LO, 

ED 0d ar e $ 0 odo 
Ba Qt 2 oo E a A 
E E E O lo 
0- DS Y A == 2 
bs Y a... 0 .t MES 

(3.15) 


6. Después de satisfacer las condiciones 
cinemáticas de borde se determina el 
vector resultante de las funciones mo- 
dales desconocidas (incógnitas en los 
nudos). 


Analicemos de una forma más detallada 
la ecuación (3.15). Por cuanto el deter- 
minante de la matriz (14) es igual a cero, 
el sistema no tiene solución, o lo que es 


lo mismo, la matriz no tiene inversa. 


Sin embargo, prestemos atención al 


hecho de que conocemos los valores de la 


Su S e 
temperatura 0, = 10 "Cy é, = -4 "CC, Por 
este motivo, en lugar de ocho incómnitas 
en el sistema de ecuaciones (1.15) exis- 
ten solamente seis incógnitas. 

En esas condiciones no es necesa Lo 


escribir la primera y la octava ecuación 
y por ello, de una manera formal pondre- 
mos en esas filas un valor unitario en el 
coeficiente de la diagonal principal y el 


resto los haremos igual a cero. 


Adewás, se hace necesario transformar 
la segunda y la séptima ecuación en la 
cual están incluidos los valores 0, y 0,. 
Para una mayor claridad, realizamos esto 


en forma algebraica y detallada. 


ps 2 ss. 
Segunda ecuación: 


Sustituyendo en estas ecuaciones los 
valores concocidos de la temperatura y 
después de pasar los valores de las cons- 
tantes obtenidas para la parte derecha de 


la ecuación, las mismas pueden ser escri- 


tas como: 


Segunda ecuación: 


0 0, +2 O, P 0, 4 s 6 


Aquí, de una manera formal fueron con- 
servados los miembros 0, y 9, con el ob- 
jetivo de no romper la estructura de la 


matriz del sistema de ecuaciones. 


El proceso descrito constituye el pro- 
ceso de satisfacción de las condiciones 
cinemáticas de borde (ei cual corresponde 
a las incógnitas fundamentales). En forma 
matricial obtenemos en lugar de (3.15) la 


ecuación matricial 


+10 4 EEE e o 0 
are 0 00 10 0 10 10 
0-1'2-1 0 0 O 0]| j0, 0 
0r 00 2he=t- 0.0 10 0, 0 
0 -0_0-+ 2-1 0: 0 e Po 
(0.010 Dier Dt. 01718, 0 
O 0, 0,8. 0:17 2 O 8 -4 
O 1014020: :0 010.3 4 0 


(3.16) 


Es fácil comprobar que la transforma- 
ción realizada de esta forma sobre la 
matriz rigidez del sistema elimina la sin- 
gularidad de la matriz, y luego de resol- 
ver el sistema obtenemos el vector solu- 
BTS US, Bf 6d, de O, +25 dr el 
cual coincide completamente con la solu- 


ción exacta de la tarea planteada. 


De esta forma, han sido descritas todas 
las etapas de solución de una tarea por 


el método de los elementos finitos en 


forma de método variacional. 


Todos los postulados del método son 
sumamente simples y pueden utilizarse sin 
mayor dificultad en la solución de cual- 
quier tarea que permita un planteamiento 


variacional. 


Es necesario señalar que, luego de la 
determinación del vector resultante de 
los valores de las incógnitas en los nudos 
en algunos casos, surge la necesidad de 
obtener los valores de algunas caracte- 
rísticas complementarias. La forma de 


hacer esto la veremos más adelante. 


de TAREA PLANA Y ESTACIONARIA 
DE LA TRANSMISION DEL CALOR 


Tenemos una región dada, dividida en ¿! 
elementos finitos, cada uno de los cuales 


tiene » nudos (ver figura 3,1). 


hamitamos que la temperatura en el in- 
terior de un elemento aislado se determi- 


na por una relación del tipo (3.7). 


Sustituyendo esta relación en la inte- 
gral (2.4) y después de realizar algunas 
transformaciones, por completo equivalen- 
tes a las realizadas en (3.8) obtenemos 
la siguiente expresión para la matriz ri- 
gidez del elemento y para el término in- 


dependiente. 
7 
CAE A 
D 
€ 
7 
am am l 
+ dad 4.1 


F =$ 5 Q[N] dxdy (4.2 
D 
€ 


donde: 


D ¿"Área del elemento. 


5, TIPOS FUNDAMENTALES DE ELEMENTOS 
FINITOS 


La forma de los elementos utilizados 
en las tareas de mecánica pueden ser las 
más diversas, además de que en cada tipo 
de elementos el número de nudos (grados 
de libertad) puede ser diferente. 


Analicemos y describamos las reglas 
para establecer las funciones de la forma 
de algunos tipos de elementos finitos de 


uso muy común. 


ELEMENTOS DE FORMA TRIANGULAR 


Constituyen la forma más simple y uni- 
versal de los diferentes tipos de elemen- 
tos finitos, ampliamente utilizada en los 
casos de aproximación de regiones no ho- 
moyéneas, de varias capas y de forma ar- 


bitraria. 


Para la descripción de estos elemen- 
tos utilizaremos un sistema normado /, 


>6 , el cual se encuentra' 


cono coordenada! 
relacionado con un sistema cartesiano 
por las siguientes correlaciones (figu- 


Ea. Sd 


2(0,1,0) 


Tb : + Rbo 
1 1 2 2 33 
=S Ly +Ly +1 .1 
y A 2% s%, ES ; 
l=al +L +1 
1 2 3 
donde: 
Tis Y¿? Coordenadas cartesianas del 


vértice 1 del triángulo: 


£ ¿"coordenadas triangulares 
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Cada conjunto de coordenadas e. Le 
E? corresponde a una pareja Única (x,y) ] 
de coordenadas cartesianas. El nudo 1 
tiene coordenadas Le =1, L, = A, =.0 y. 


de forma semejante los restantes nudos. 


El vínculo o relación lineal estable- 
cida entre los sistemas de coordenadas 
en las fórmulas (5.1) significa que las 
líneas ÓN = Constante representan rectas, 
paralelas al lado opuesto al ángulo 1, 
en la cual Li = 0 (ver las líneas L, 
sobre la figura 5,1 para los valores 1/3 


y 0,5). 


Desde un punto de vista geométrico bi 
representa la relación entre las Éáreas o 
: : 6 
las distancias hasta el lado opuesto”. 


Por ejemplo para el punto p (fig. 5.1). 


L =Jistancia Ap _ Área del triángulo 2p3 
1 distancia 41 Area del triángulo 123 
(532) 


Resolviendo el sistema (5.1) con rela- 


ción a bi obtenemos que: 


(5.3) 


donde: 


-” 
y L, ¿| ): Area del ele- 


Y, Yy Y mento. 


il valor abs indica que lo que intere- 
sa es el valor absoluto de la magnitud ob- 
nida: UU =x YU 4 1 b = - si 
. 1 5S 4, 4 y 1 Es q 


ce. =£ -= 2 , y los restantes coeficien- 
2 


1 3 
tes se okLtienen por un cambio cíclico de 


los índices 1,2,3. 


Para el cálculo de las integrales del 


producto, las coordenadas del tipo 


Sue A, Ley detarróllan sórmulas cónos 


á : 
das”, 


En el caso de tareas lineales y planas 
la fórmula generalizada de integración se 


puede presentar en la forma: 


A 48078: mn (o UL 0% < M0 CS y 
A a ÓN 


donde en dependencia de la dimensión de 
la tarea los parámetros du, F, y 5 tienen 


el siguiente significado: 
Para el caso unidimensional: 


p= d),= S 


tj due = de; 
Para el caso bidimensional: 


rs = 25 du) = uxcdly y lU= A 


vz Longitud del lado del elemento trian- 
gular con respecto a la cual se reali- 


za la integración. 


Para el elemento de primer orden con 
tres nudos en los vértices (fig. 5.1) las 
funciones de la forma coinciden con las 
coordenadas L: 


[WS EVA NV E AD a 
o a 


. 


La matriz [/] tiene tres miembros, así 
como cada elemento tiene tres grados de 
libertad. Es fácil comprobar que las 
funciones Li tienen las mismas propieda- 
des de las funciones de la forma, descri- 


tas anteriormente. - 


Para el elemento de segundo orden, el 
cual tiene seis nudos (fig. 5.2) las fun- 
ciones de la forma se expresan a través 
de las coordenadas L en la siguiente 


forma: 
MW = j m7 j 
[M A + Ye Ns. Yin Y. Ni 
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E 


[M =1 (2L 02, ALL ¿(21 -1)L, «ALL, 6 


(2L -1)L ¿41 L] (5.6) 
3 3 NS | 
3 
2 4 
| 5 
6 


«Figura 5.2 


De forma similar se obtienen las fun- 
ciones de la forma de elementos triangu- 
lares de orden superior. Una fórmula re- 
currente que permite obtener las funcio- 
nes de la forma en elementos triangulares 
de cualquier orden se presenta en el tra- 
bajo! . 

En la práctica, sin embargo, comúnmen- 
te se utilizan elementos de primer y se- 


cqundo orden. 


Es conveniente señalar que las funcio- 
nes de la forma para estos mismos elemen- 
tos se pueden obtener sin la utilización 


del sistema de coordenadas L. 


Para el elemento triangular de primer 
orden podría adoptarse “a ley lineal 
para el campo de las temperaturas en la 
forma: 

5.74) 


0(x,y) =4a +az+al 
(uy) : 5 pS 


En esta fórmula a, ( = 1,2 ..) son coe- 


ficientes indeterminados. Satisficiendo 
las tres condiciones de korde MER Y) = 
= o, (1 = 1,2,3) obtenemos un sistema de 
tres ecuaciones, el cual después de una 
serie de transformaciones, semejantes a 
las realizadas para el elemento lineal 
(3.1-3,5) se puede obtener las mismas 


¡funciones (5,5). 


Para el elemento de segundo orden, el 
campo de temperaturas se expresa a través 


de un polinomio de segundo grado; 
0 = HAD RAY AAA amy 
(1,4) e e eE pa E 


2 [el 
+ ay (5.8) 


a partir del cual, y con el objetivo 

de obtener las funciones (5.6) sería ne- 
cesario resolver un sistema de seis ecua- 
ciones, con seis incógnitas en el cual 
las funciones expresadas algebraicamente 
tendrían una forma complicada, como es 
fácil comprobar sustituyendo (5.3) en 
(5.6). 


De esta forma, la ventaja fundamental 
del sistema de coordenadas triangulares 
consiste en la posibilidad de brindar una 
expresión compacta y explícita de las fun- 
ciones de la forma para los elementos 
triangulares de cualquier orden en forma 
de una combinación de los productos de 


las coordenadas L;. 


Además, la utilización de la fórmula 
(5.4) simplifica de forma considerable el 
proceso de obtención de las expresiones 
finales para las matrices rigidez y los 
términos independientes de los elementos 
correspondientes, lo cual representa una 
simplificación en comparación con las re- 
presentaciones polinómicas (5.7),(5.8) 
que: aún se encuentran en algunas publica- 


, An 11 
ciones soviéticas? *'*'; 


ELEMENTOS RECTANGULARIS 


Para la solución de problemas en re- 
giones no homogéneas de forma rectangu- 
lar, resulta muy cómodo utilizar el ele- 
mento finito mostrado en la figura 5.2. 
Si se introduce un sistema normalizado de 


coordenadas tenemos que: 


-1<SE<1, -159<1 
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Figura 5.3 


De los anteriores valores ligados al 
centro de gravedad del rectángulo se 
puede obtener la siguiente fórmula de las 


funciones de la forma: 


1 ; . 
a O O 


4) 
E” 
- 
[3 
- 
uu 
- 
> 
— 


Sustituyendo aquí los correspondientes 
valores de MN; Y Es s fácil obtener las 
expresiones concretas de No En coordena- 
nas cartesianas, estas mismas funciones 


pueden obtenerse del polinomio: 


0(x,y) =a +az+a y +azy (5.10) 
1 2 3 4 


- 


Es fácil darse cuenta que la variación 
de la función 0 («,y) a través de los lados 
del rectángulo es lineal, y significa que 
el elemento mostrado (fig. 5.3) tiene por 
completo la misma significación que el 


elemento triangular de primer orden. 


La utilización simultánea de elementos 
triangulares y rectangulares en el cálcu- 
lo brinda cierta economía de tiempo de 


máquina en los algoritmos. 


ELEMENTO EN FORMA DE TETRÁGONO* 


Este elemento (fig. 5.4) se utiliza 
mucho menos que los elementos de formas 
triangulares. 

*Tetrágono: Elemento cerrado, de cuatro 


lados rectos, dispuestos arbitraria- 
mente. 


Figura 5.4 


Sin embargo, él posee la comodidad del 
triángulo y en muchos aspectos es similar 
al rectángulo durante la división de una 


región cualquiera. 


De forma similar al caso anterior, in- 
troducimos un sistema de coordenadas nor- 
malizado (f.,n) con ayuda del cual las fun- 
ciones de la forma pueden ser descritas 
por la propia fórmula (5.9), pero en este 
caso, y para poder determinarse por com- 
pleto este elemento, es necesario esta- 
blecer el vínculo entre los sistemas de 
coordenadas, el cual puede ser escrito en 


la forma: 


(5.11) 


La relación entre las derivadas con 
relación a las coordenadas (£ ,n) y (%,4) 
se determina por las siguientes expre- 


siones: 
d 00 d7 a 
Pr artatiy 


0 dr 3 dy 3 - 
An” on da * 3n ae) 
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o en forma matricial: 


2 el 
E A dE 
dy 7) 


donde: 


.. 
o 
o 


[7] 


ox 0 > 
+ 

-(1)  - (18) (19) (1-7) +A, 
a e-y 
= 8 Pr 
=(1-E) (1-£) (14€) - (1+É) ¿e 
w y 
4 4 
(5.14) 


El diferencial del área se determina 


por la expresión: 


dxdy = det (| 71) d ¿du (5.15) 


Elementos de otro tipo, como por ejem- 
plo, elementos curvilíneos, se describen 
en los trabajos *”. Sin embargo, en la 
mayoría de las tareas técnicas de análi- 
sis de campos en el plano son suficien- 


tes los tipos de elementos descritos. 


De forma similar se construyen las fun- 
ciones de la forma para elementos tridi- 
mensionales en forma de tetraedros y 
prismas, de elementos de placas y cásca- 


ras' s4 = 


En todos los casos es necesario conser- 
var el procedimiento indicado de obten- 


ción de las matrices [/]. 


UL, RELACIONES DEL MÉTODO DE LOS ELEMENTOS 
FINITOS PARA EL ELEMENTO TRIANGULAR 
EN LAS TAREAS DE TRANSMISIÓN DEL CALOR 


Analicemos los aspectos del cálculo 
práctico de las matrices rigidez y de los 


términos independientes de los elementos 
bidimensionales. Las derivadas de las 


coordenadas L se determinan por las rela- 


ciones: 
ob, b; aL, Co : 
Pa AA 


(6.1) 


entonces, calculando consecutivamente las 
expresiones «ue aparecen en la fórmula 


(4.1) obtenemos: 


b 
e : 
A 0". 1. ¿DA b VMbbb] = 
Jd z e 2 > 
b 
3 
E SIM 
A 1 
ES : b, d, b, 


DB: bibi Pe 
3 a; 3 3 


De forma similar se tiene: 


e? SIM 
a ; 
A 014] 22 De oÑ 5 da 
dy y ae 17 2 


donde : 


SIM: Simetría de la matriz, su valor 


és igqúal. da es E dos 
z 1 JT 


Admitamos que el coeficiente de trans- 
misión del calor Á es igual a una cons- 
tante en el interior del elemento finito, 
por lo que todos los coeficientes A, Bin 
e¿AÁ serán independientes de las coordena- 
das. Teniendo en cuenta que Sfdxdy = A 

D 


3 
de forma definitiva obtenemos que: 


CE, 
:b +0; SIM 
A NA 
= > " bi +e 
LA] 44 OS 2. 2 


bbrelbbrece be 
E O A E E e 


(6.2) 
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Pasemos ahora al cálculo del vector 


Fo (452) 


de los límites del elemento, tenemos que: 


Suponiendo (Q) constante dentro 


AE O 
le dxdy = y L; dudy = cd ARTO TAS 
y: 

1 
F. 0 (6.3) 


La cual es equivalente a las expresio- 
nes Obtenidas para las matrices [4] y A 


de elementos triangulares de mayor orden. 


/, TAREA NO ESTACIONARIA DE LA TRANSMISIÓN 
DEL CALOR 


Analicemos como varía el procedimiento 
por el método de los elementos finitos 


cuando la función investigada depende del 


tiempo. En este caso en el lugar de la 
ecuación (2.1) tenemos la ecuación: 

90 ; ; 
Co sp = 400 + 4 (2,Y,€) (7.1) 
donde: 

C + Calor específico volumétrico del 


material. 


La función de la temperatura U(EjUintde 
además de satisfacer las condiciones de 


borde: 


0(2,yt) l =p(0,y,t) (7.2) 


debe de forma complementaria satisfacer a 
la condición inicial 


0 (2,4 ,0) = Y (c,y) (13) 


Además a la integral (2.4) se le adi- 
ciona un miembro que tiene en cuenta el 


cambio de la temperatura en el tiempo. 


1 90 2 00 ; 
As dr le 


= (Q - Cp 00 Jatedy = da + 


+S50, $20 dedy (1.8) 
D 


donde: 
¿y : Integral (2.4). 


Tenemos así que dentro de los límites 
del elemento finito aislado D, la tempe- 


ratura se determina por la expresión: 
0 (eya t) =[N (2,y)) o. (£) (7.5) 


Y) 
transformemos el sumando con ER en (7.4) 


o (uo, ) 


110 qe 0 dudy OS —= A M0, dudy= 


gl SÍ ye E AE 0, ES : d. 
A 
e 


donde: 


1Cl = 550,1 Mé LW] dxay Matriz del 
€ 
calor específico del elemento (MCE). 
Entonces la forma discreta de la inte- 
gral (7.4) para toda la región de una 


forma definitiva toma la forma: 


q ¿de 0 


(7.6) 
Aquí [Cl] a [Cc] ¿ es la matriz del calor 
específico del sistema de elementos (MCE). 


La condición de mínimo de la integral 
(7.6) vendría dada por la ecuación di fe- 
rencial matricial de primer orden. 

7 
od 


= [€] 24 1495 +F=0 (7.7) 


A partir de la solución de la cual se 
obtiene el vector 0 en los nudos de la 
malla desarrollados a través de todo el 


intervalo de tiempo. 


Es necesario señalar que a diferencia 
de las tareas estacionarias el sistema 
resultante (7.7) constituye un sistema 
de ecuaciones diferenciales y no un sis- 


tema de ecuaciones algebraicas. 


Determinemos la matriz [Cl para el 
elemento triangular de primer orden en la 
suposición de que C,, = constante dentro 


del elemento. 


L0l, =SSC XL LL] dady = 
€ 2 1 2 3 
De 
7) 
5 l sim 
A | 
== CS VE dedy = 
D | 
diia 
DAA RS. 1 S 
' 1 
c,A 2 SIM 
a (7.8) 
A: 


> Diuxdr 2 1 2 Á MN A 
1 Q+2N 6” 


% Ae ra BA A 
4 - ua = zz —— 
JS BE cuy OFTEN 12 


de 


En lugar de la matriz completa [Cl 
calculada a través de la fórmula (7.8), 
en los cálculos se puede utilizar una 
matriz diagonal, formada a partir de la 
matriz (7.8) por la suma de todos los ele- 


mentos de cada fila. 


DAFO (7.9) 


Numerosos cálculos numéricos, realiza- 
dos por Diomin a MOStrao como la 
matriz (7.8), al igual que la matriz (7.9) 
dan resultados, prácticamente de igual 
grado de precisión. 


Sin embargo, la utilización de la ma- 
triz (7.9) permite fácilmente transformar 


el sistema 7.7 a la forma normal: 


00 PR 
E 1H 9 + E =0 (7,10) 
donde ; 
-1 iS se 
Lal, =[C1, 14M E [el 7 


y elaborar un programa de computación 
económico haciendo uso del método de in- 


tegración de Runge-Kutta. ' 


En el caso de la matriz (7.3) realizar 
la normalización del sistema no resulta 
económico pues se exige resolver el sis- 
tera (7,7) por el método de las diferen- 
cias finitas y el programa que se obtie- 
ne trabaja de una forma considerablemente 
más lenta que en el caso de la matriz dia- 
gonal. Por 6só; durante la resolución de 
tareas prácticas se recomienda la forma 
diagonal de formación de la matriz calor 
específico del elemento, la cual a pesar 
de su insuficiente fundamentación teóri- 
ca, permite, en general, obtener precisos 


resultados. 


8, MÉTODO DE LOS ELEMENTOS FINITOS 
APLICANDO EL MÉTODO DE GALERKIN 
(MEFG)L 5 


La forma variacional del método de los 
elementos finitos expuesta en los párra- 
fos anteriores constituye hoy en día la 
forma más extendida y cómoda de aplica- 
ción del método. 


Sin embargo, en algunos casos surge la 
necesidad de obtener las relaciones del 


"método de forma directa a partir de las 


conocidas ecuaciones diferenciales no re- 
curriendo al desarrollo del planteamiento 
variacional, el cual puede resultar en 


general desconocido para la tarea dada. 


Analicemos los postulados fundamenta- 
les del método de los elementos finitos 
utilizando el nétodo de Galerkin* en el 
ejemplo de la tarea bidimensional de la 
transmisión del calor (2.1), 2221. Es= 


Ccribamos la ecuación (2.1) en la forma: 
¿(0) =0 (8.1) 


donde: 


É : Forma simbólica de escribir el co- 
rrespondiente operador diferencial. 
Tenemos la región desarrollada en 
élenentos finitos y en el interior 
de cada elemento aislado la tempe- 
ratura se determina por la expre- 


sión (3.7). 


Por cuanto la ecuación (8.1) es apli- 
cable también para los elementos aislados, 
después de sustituir (3.7) en (7.9) te- 


nemos ; 
EM 0,)=RFO0 (8.2) 


Así como la sustitución de la solución 
aproximada en la ecuación exacta es nece- 


sario dar tolerancia R. 


El método de Galerkin demuestra cque la 
mejor solución será aquella que dé en 
todos los puntos de elementos un valor 
mínimo de la tolerantia R. SiR=0en 


toda la región, entonces: 


erddy = 0 (4.3) 
€ 


donde: 


W: Función arbitraria de las coordena- 


das (función pesada). 


POr cuanto el Hiánaxo de parámetros des- 
conocidos en el elemento es feñal «1 núme- 
ro de sus nudos, en calidad de función pe- 
sada es cómodo seleccionar las funciones 
de la forma [WN]. Entonces en correspon- 
dencia con (8.3) tenemos: 


LT años quirdy =00 (8.4) 


Disminuyamos el orden de las integra- 


les en la ecuación (8.4) integrando por 


partes. 
00 
gm! LL dea Y AMÓ gr y 
De E a 
yl 
SS do 7 dy e 
D, 
qm L de - 1 AU ass 


y de forma similar tenemos que: 


ír 90 
¿Mr 


7 : 
dy Lar] dedy= $ Las = 
[e , : 


Aa ¿nn >,. 
ia y dy duty 


Lo 
apnea e: 


L. y ly + ad directores de la nox- 


mal A a la frontera es 


Entonces de la ecuación (8.4) tene- 


ros que: 
15 e Ml o h, 


dE 1) 
- [MM "Qdwdy - $ ¿1 MAL, + 
A de “z 


00 


+ 


L,yd8 = 0 (8.5) 


Sustituyendo (3.7) y teniendo en cuen- 


ta que la expresión siguiente se repre- 
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senta así: 


00 00 00 
dx La, + L, == 
$ dn 
tenemos que: 
— = - P. OU 
[410 + ES y [NM pd as =0 (8.6) 
€ ón 


Después de la sumatoria a través de 
toda la región obtenemos un sistema de 
ecuaciones semejante a las ecuaciones 
(3.10), con la excepción de que le adi- 
ciona la integral a través de la frontera 
de la región. Por cuanto nuestro análi- 


sis se limita a la primera tarea de 


torde (2.2), entonces AA O y el 
ón 

sistema obtenido es por completo equiva- 

lente al sistema (3.10). 


Y, IPLICACIÓN DEL MÉTODO DE LOS ELEMENTOS 
FinrTOS (EFG) A LA SOLUCION DE 
ECUACIONES DIFERENCIALES ORDINARIAS 


Analicemos en calidad de ejemplo la 


ecuación diferencial ordinaria: 


ay? + by +. =:0 (9 1) 


donde: 


d, b y e: Constantes. 


Analicemos el intervalo 0<x<€ l, en 
los extremos del cual se conocen los va- 
lores de función (y), y desarrollemos 
este intervalo en una serio de elementos 
finitos lineales (fig. 9.1). 


en cuenta que (y) en el interior del ele- 


si tomamos 


mento varía linealmente, después de una 
serie de transformaciones por completo . 
semejantes a las realizadas en las fóxmu- 
las (3.1) y (3.5) tenemos: 


= |] y (9.2) 


donde: 
e z ÁAÁhx-ox Ez 
[Wi =[M,. Ni =[ Az , ET 
ón Yo 
Ax -2.z/,ys. 


SON 


La derivada de esta función se calcu- 


la en la siguiente forma: 


(9.3) 


El valor inicial de yy Para x =0 se 
supone conocido, de forma tal que en las 
relaciones (8.6) la única incógnita es y, . 
Por lo tanto, se puede utilizar una sola 


función pesada N, . 


dl 
taz ' 


a O 
A 


a lay” + by + cldx =p 


an, A e 
e del +b[N y +c) de =0 


Donde después de sustituir y realizar 


la integración se obtiene que: 
1 


CÁ z 


a DA x 
2 Ye 48d AH A) + Ela 
Si se resuelve esta expresión con re- 
lación a Y,» Obtenemos la fórmula recu- 
rrente: 


: l (a - 22 


)y,7 CA al 
(a + —ÑbA 2) 2 


(9.4) 
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Este resultado es semejante al que fue 
obtenido con ayuda del método de las di- 
ferencias finitas, usando el esquema de 
la derivada central, sin embargo, la de- 
ducción de la ecuación (9.4) es otra. 
Esto permite, por ejemplo utilizar otras 
ecuaciones de interpolación. En parte se 
puede analizar el elemento del tipo mos- 
trado en la figura 3.4 a) y obtener fór- 
mulas más precisas, lo cual no es posible 
a través del método de las diferencias 


finitas. 


De forma totalmente semejante el méto- 


do MEFG se aplica a la solución de siste- 
mas de ecuaciones diferenciales. De 
igual forma, para el sistema (7.7) des- 
pués de la realización de estas mismas 
transformaciones en el pequeño interva- 


lo Á € obtenemos consecuentemente que: 


- 28 , 00 1 1,1] : 
A O E ai e 
1 1 
y. E 
uN ZE yy > AL 
ON ON z 
At t ? 0 1 
[ ll 7 + 


+ [A4J1 + 4) + [Fhdt =0 


donde, después de las sustituciones y la 


integración se obtiene que: 
ie pa es +0. J+ 
EA AA St cd. 


At= 
a 
ap 0 
* 
o según la referencia | 5] : 


DA 2 hal. 4 
Y MAA SES LESA 
¿ LO 57 2_ ¡A 
O, + rr A 


Esta fórmula es la más cómoda en el 
caso de utilizar la matriz completa de 


calor específico del elemento. 


10, TAREA PLANA DE LA TEORÍA 
DE LA ELASTICIDAD EN FORMA MATRICIAL 


El campo de los desplazamientos 
u(x,u), de las tensiones 0 (u,4) y de las 
deformaciones € (1 ,/) en cualquier región 


plana se determina por los vectores si- 


guientes: 

mM u (x.y) a 9, (y) 

u(uy) = o (ey) y 
V (uy) "ey ey) 

q E Y) 

€ (UY) =9 €, (1,4) (10.1) 
Yazy rio) 


Aquí, u(c,y) y V(<,Y) representan res- 
pectivamente los desplazamientos horizon- 
tales y verticales de un punto arkitrario 
de la región. Todos los indicadores del 
estado tensional y deformacional (FTD) de 
la región se encuentran ligados por cono- 
cidas relaciones'? las cuales presen- 


taremos aquí en su forma matriciol. 


ECUACIONES Di LQUILILIO 


(a]* 0+p=0 (10.2) 
donde: 
d 0 
0x 
[41 =l 0 ST : Matriz de operadores 
A a diferenciales, 
O “O 
po.= : Vector de las fuerzas mási- 
y 
PY cas. 


CONDICIONES DE BORDE ESTÁTICAS 


[4,1% 5 =E (10.3) 


4 


donde: 
cos (ma) 0 

(4, =| 0 cos (ny) : Matriz de 
cos (my) cos (ree) | los cosenos 

directores 
R 
“y 2% 
R= : Tensiones en la frontera de la 


y 


KR ) región, 


'Az Vector normal a la frontera SÍ 


de la región /'. 


- nl . 
mu y ny: Ángulos que forman la nornal 
y Py g 4 


con los ejes coordenados . 


FÓPMULAS DE CAUCHY (ECUACIONES GEOMÉTRICAS) 


= [AJu (10.4) 


m1 


LEY DE HOOKE (WCUACIÓN FÍSICA 
O CONSTITUTIVA) 


o =[0]€ (10.5) 
donde: 


(A +26) A 0 
dz A (A + 2G) O 
G las relac. 


de la 


: Matriz de 


0 0 


elastici- 
dad. 


A : Coef. de 


e 


Lame . 


en el cual: 


Ey / (1-2 1) (1 + 17) 
Ne= :* Estado deforma- 
EXI - Y”) cional plano y 
estado tensional 
plano generali- 


zado. 


E: Módulo de Young. * Y: Coeficiente de 
Poisson. 
AN 
2(1 +7) 


[e 


CONDICIÓN DE SAINT-VENANT PARA EL CASO 
PLANO (ECUACIÓN DE COMPATIBILIDAD 
DE LAS DEFORMACIONES) 


ate. e Y. 
ME + A AA (10.6) 
ay da? du dy S 


En calidad de grupo de incógnitas fun- 
damentales, seleccionaremos el vector de 
los desplazamientos U(ey). 

Introduciendo en la ecuación (10.2) los 
componentes de la tensión (10.5) y cam- 
biando el vector de las deformaciones por 
la fórmula de Cauchy (10.4) obtenemos las 
ecuaciones de Lamé en forma matricial. 


A LLOL 141) +. =0 (10.7) 


Entonces, para medios isótropos planos, 
la tarea estacionaria de la teoría de la 
elasticidad se plantea en la siguiente 
forma: encontrar el vector de los despla- 
mientos U (uy) determinado y continuo en 
la región cerrada D junto con sus deriva= 
das parciales de primero y segundo orden, 
las cuales satisfacen en el interior de 
esta región a la ecuación diferencial 


(10.7) escrita en forma matricial, así 


como las condiciones de borde matriciales. 


[A,1 101 Lal Z) =2 (10.8) 


11, ENERGÍA POTENCIAL DE DEFORMACIÓN 


Analicemos una cierta región plana D 
(figura (11.1) soretidas a la acción de 
fuerzas distribuidas q = tq, 0d y fuer- 
zas concentradas P,. Las fuerzas másicas 
se determinan por el vector PP. Calcule- 


mos el trabajo de todas las fuerzas exte- 


=> 


riores e interiores aplicadas en la re- 
gión dada. El trabajo elemental de las 
fuerzas exteriores distribuidas. 


(0) 


Figura 11,1 


(11 1) 


/ =i n + ) SS = 
dA, (qu dy vja 


El trabajo elemental de las fuerzas 
másicas: 
Up daudy 


dA = 


% 017.2) 


o + a, Wdzdy = 


El trabajo de las fuerzas concentradas 


(exteriores) R 


(11.3) 


donde: 


5 .: Desplazamiento en el punto de apli- 
% 
cación de la fuerza P, en su propia 
dirección. 
J: Número de fuerzas concentradas apli- 


cadas en la región. 


El trabajo elemental de las tensiones 


interiores es: 
1 -T- 
dA, DE de o dxudy 


La influencia de la temperatura provo- 
ca en la región tensiones de temperatu- 
ra [13]0, =- [D]e,. 
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AA " 


donde : 


E Vector de las deformaciones provo- 
cadas por el cambio de tempera- 


turá. 


A partir de lo cual el trabajo total 
de las tensiones interiores se expresa 


por la fórmula: 


N (11,4) 


da, = — e odady + e “odody 


La energía potencial total de deforma- 
ción se determina por la suma de las inte- 
grales de los correspondientes trabajos 


elementales: 


peo (Ly oe nto ip rdedy - 
D 0 


Expresemos esta fórmula a través de 
los desplazamientos. Utilizando (10.4) 


obtenemos: 


po 5 aan 10) aja) - «aji? - 
D 


- 2p ) ddxdy -S u ode - z p.5. (11.6) 
S a A 
Entonces el planteamiento variacional 
de la tarea de teoría de la elasticidad 
es el siguiente: determinar el vector de 
los desplazamientos u(ey)» determinado 
y continuo en la región cerrada D junto. 
con sus derivadas de primero y segundo 
orden, el' cual satisface en la frontera 
de la región las condiciones (10.8) y 
que le comunica a la energía potencial 


de deformación (1.16) un valor mínimo. 
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12, VEMOSTRACIÓN DE LAS RELACIONES 
FUNDAMENTALES DEL MÉTODO 
.DE LOS ELEMENTOS FINITOS 
PARA LAS TAREAS DE TEORÍA 
DE LA ELASTICIDAD - 


Se tiene una región dada D, desarro- 
llada en los elementos finitos | ver 
fig. 3.1 a)]. A cada nudo (1) del ele- 
mento se encuentran vinculados los corres- 
pondientes desplazamientos horizontales 
u, Y verticales aa lo que representa que 
a diferencia de las tareas de transmisión 
del calor, cada nudo posee no uno, sino 
dos grados de libertad. El vector de los 
desplazamientos 5, para el nudo Z es 


igual a: 


(12.1) 


El vector resultante de los desplaza- 
mientos de los nudos del elemento se re- 


presenta como: 


(12.2) 


Ol ... 
pa 
SS 


La aproximación del campo bidimensio- 
nal de los desplazamientos Se realiza por 


la misma función de la forma: 


u(x,y) IA +... + Na 
u (xy) e 


V (uy) N_p HEo: to. + A 


N N O0.NM 
, 9 NM 0.80 


A E (115, 
ON ON. 0Nn 
2 rn 


+ (12.3) 


Sustituyamos la expresión resultante 
(12.3) en la fórmula (11.6) e introduzca- 


mos la notación: 


138] =14] 14) (12.4) 


de esta forma: 


e E is a 
SE TEMA, ELSE, = 
€ 


- 51m daedy = 5.1 mMqds pa BP. 
3 
€ 


(12597 
donde : 


P_: Vector de las fuerzas concentradas 


en los nudos del elemento. 


Volvamos a escribir la ecuación (12.5) 


en la forma: 


A 
SS IA +8 +8, E, + 
=T 71- 1 ¿7 =T = 
ñ e ia” ES a RIA > De a 
(12.6) 
donde: 


[XM1, = 50811011 Bldedy: Matriz rigidez 
€ del elemento 


+R 


pei 


+ E 
eo E, 


E A 
SS1B1"[D]e,dudy; 


D, 


B, =- 150 dudys E, == (m%gds 
D, ss 3, 


"> 
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La energía de toda la región es igual 
a la suma de las energías de todos los 
elementos. Entonces la energía potencial 


de toda la región: 


R* = 51M 54 52 (12.7) 
donde: 
[XK] =2 [XI,: Matriz rigidez del sis- 
4 
tema (MRS). 
R=Y Ro: Vector resultante de los 
dd términos independientes. 
Mz Múmero de elementos. 
e Vector de los desplaza- 


mientos nodales de toda 


la región. 


La condición de mínimo de la función 


12.7 es la siguiente: 


(12.8) 


Después de satisfacer las condiciones 
de borde de este sistema de ecuaciones 
algebraicas lineales se obtiene el vector 
resultante 5. Para la solución total de 
la tarea sobre la determinación del esta- 
do tensional y deformacional del sistema 
en la región dada es necesario encontrar 
los vectores U y €. Esto se realiza de 
forma separada para cada elemento de la 
región con ayuda de las fórmulas (10.4) y 


(10.5). 


_m 
" 


[4] 4, = 14J101 5, = 13] $, 


a1 
Ú 


[D] €, - €, ) = 10148] $, - €, ) 
€ € 


(12.9) 


r 


A través de los vectores o. y E. con 
ayuda de las fórmulas habituales de la 
teoría de la elasticidad pueden hallarse 
las tensiones principales, las deforma- 
ciones principales, el ángulo de inclina- 

ción de las tensiones principales y otras 
características necesarias sin ningún 


tipo de complicaciones teóricas. 


15, RELACIONES DEL MÉTODO DE ELEMENTOS 
FINITOS PARA EL ELEMENTO TRIANGULAR 
EN LA TAREA DE LA TEORÍA 
DE LA ELASTICIDAD 


El elemento de primer orden posee seis 
grados de libertad. La matriz de las 
funciones de la forma en correspondencia 
con la ecuación (12.2) se escribe de la 


siguiente forma: 


(13,11) 


Las derivadas de las coordenas L se 


determinan por las relaciones (6.1). 


Entonces la matriz [B] en la fórmula 
(12.4) toma la forma que se presenta a 


. . 
continuación : 


Para el elemento de primer oxden la 
matriz [B] no depende de las coordenadas. 
Por tanto, la expresión para la matriz 
rigidez del elemento se simplifica de la 


forma siguiente: 


141, = 181410] LE S5dedy = LA 70 [8] A 
D 
€ 

(13.3) 


Después de la multiplicación de las 


matrices, de forma definitiva obtenemos 


A E O E 
1 3 que: 
[4], = 
B b*o4 Ge? 
1 1 
, 2 3.2 
De + 
8, E 1 Pe $e SIM 
Bbb +G6co Beb +Gbe Bb «de 
SA LL =- y Ur 2 tt 2 y 


1 
= paa + 68/b; pe e +Gbb Bbe 


pe + ap 
£- 7 y e: M7 y 2 2 


BDbD +Gee Beb +*be Bbb +Gcec Bob +Ghbc. BDsG 
4,1, 12 i 3 dd 2103 > a A 2.3 d .2*S 2 3 


Bba +6eb Bac +Gbb Bbe +Geb Bcc +Gbb Bbe Be + Gb 
$ 1 3 1.3 EAS 13 AS 2 3 1 IE De a 333 E 3 3 
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En este caso tenemos que PB, =A + 26; 
es =A + G; B, =A, 


pueden ser obtenidas las matrices rigidez 


De forma semejante 


del elemento para elementos de segundo 
orden y de orden superior y en este caso 
surge la necesidad de integrar a través 
de las coordenadas L, lo cual se realiza 


a través de la fórmula (5.3). 


Calculemos los componentes del vector 
resultante de los términos independien- 
tes Ro + 

Bajo la acción de fuerzas másicas cons- 


tantes en el elemento tenemos: 


La expresión para R, la obtenemos en 
la suposición de que el vector q es cono- 
cido a través del lado 1-2 del triángulo 
(ver fig. 5.1). Supongamos que 


4 = const. y Uy = const. , entonces: 


El] 
" 


(13.6) 


£ e Longitud del lado 1-2. Los dos 

últimos ceros en el vector apare- 
cen a partir del hecho de que en 
el lado 1-2 del triángulo L, = O, 
que también se tiene en cuenta en 


la matriz [NM]. 


En el caso de influencia de la tempe- 
ratura el vector de las deformaciones ini: 
ciales por temperatura puede ser escrito 


en la formma.' »*? 


E..=f0)1 013.7) 


donde: 


f = 0%: Estado tensional plano gene- 


ralizado., 
f =(1 + Y7)4: Estado deformacional plano. 


0 (u,4): Temperatura del interior del 


elemento. 


%: Coeficiente de transmisión 
lineal de la temperatura en 


el material dado: Vos = 0, 


ya que la dilatación libre 
por temperatura no origina 
deformaciones angulares.  Te- 
niendo en cuenta que en el 
elemento la temperatura se 


aproxima por la f£óemu- 


LR SI 
> a 0 13.8 
0 (2,y) = L, > + b, > + E A ( ) 
obtenenos: 
= Tr, a 
R__ =-[B]"[0] S$f€dxdy = 
€0 
D 
€ 
7 1 
== fIB 10155 0dmy = 
O Le 


b 
6 
=- 6 A+G 
EN ep" 1/b 
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Cy (13.9 


Lás deformaciones del elemento se de- 


. 


terminan por las fórmulas: 
bu wtbur+bu 

A ES 

E CvU+eov +v 
24 TS ERE 
cubo seu s+bvseus+bw 

1,51 ii | 2-2 254 y 3 313 
(13.10) 


Las tensiones se determinan por la fór- 


mila: 


(A + 2G)e . + Ge ( 
y Y , 


ue, + a + 20, ZO ie 


il 


1 
€ 24 


(1327) 


Como se desprende de las .fórmulas 
(13.10) y (13.11) los vectores de las ten- 
siones y de las deformaciones no dependen 
de las coordenadas, es decir, que los mis- 
mos son constantes dentro de los límites 
dél elemento. Esto es consecuencia de 
que se supone una condición de Compatibi- 
lidad simplificada, esto es, que solo las 
funciones u(u,y) se exige que sean conti- 


nuas. 


Habitualmente, y teniendo en cuenta 
esto, en los cálculos prácticos se supone 
que los vectores o, y e, están aplicados 
en el centro de gravedad del elemento. 
Buenos resultados se obtienen con el 
cálculo de los valores medios de las ten- 


siones y las deformaciones en el nudo. 


2 E Y E 
- Mo € - Me 
A A (13.12) 
donde: 


(1): Número del nudo. 


(M) : Número de elementos que forman 


parte del nudo analizado. 


Ji, REALIZACIÓN DEL MÉTODO 
DE LOS ELEMENTOS FINITOS EN LH! 


La realización de los cálculos del mé- 
todo de los elementos finitos exige la 
realización de un gran volumen de opera- 
ciones y puede realizarse solamente con 
la ayuda de máquinas computadoras (IBM) « 
Los postulados presentados en los epígra- 
fes anteriores del método de los elemen- 
tos finitos es necesario analizarlos como 
el aparato teórico de solución de la 
tarea planteada, el cual requiere ser 
transformado en un programa de computa- 


ción para una máquina concreta. 


Desde el punto de vista de la computa- 
ción el método de los elementos finitos 
es uno de los métodos más voluminosos de 
solución de tareas aplicadas. Para dis- 
minuir el gasto de tiempo de máquina y 
elevar la precisión de la solución te- 
niendo en cuenta la utilización de un 


elevado número de puntos nodales, es ne- 


«cesaria la utilización de máquinas con 


una alta velocidad de operación y una su- 
ficiente capacidad de memoria. Una gran 
significación tiene la existencia de ins- 
talaciones de memoria exterior como kan- 
das magnéticas, discos, tambores maynéti- 


cos y Otras. 


En este epígrafe, examinaremos aspec- 
tos generales de realización del método 
de los elementos finitos en IBM en el 
ejemplo de las tareas de transmisión del 
calor. Para la descripción de diferentes 
aspecto del programa se utilizará el 
idioma ALGOL 60 que es el más extendido 
en nuestro país hasta la fecha (URSS) - 

En caso de ser necesario todo el procedi- 
«miento que se describirá puede ser fácil- 
mente traducido al caso del idioma 


FORTRAN. 


En el presente, tanto en nuestro país 
como en el extranjero' se poseé una gran 
cantidad de programas basados en el méto- 
do de los elementos finitos. Al genera- 
lizar la experiencia de su realización es 
necesario señalar que a pesar de que la 
solución de una tarea concreta de cual- 
quier aspecto de la mecánica tiene sus 
rasgos distintivos, a todos los programas 
se les puede aplicar una serie de aspec- 
tos generales, que son producto de las 
particularidades propias del método de 
los elementos finitos y son aplicables a 
la descripción de la estructura de todos 


estos programas. 


Limitando nuestro análisis solamente 
a las tareas lineales en el campo de la 
"mecánica, se puede proponer el siguiente 


diagrama de bloque (fig. 14.1) el cual en 


COMIENZO 


principio generaliza a los programas de 


computación señalados. 


A pesar del hecho de que cada bloque 
se realiza de una manera distinta en de- 
pendencia del tipo de tarea y de las di- 
ferentes variantes de realización de dis- 
tintos aspectos del método de los elemen- 
tos finitos, se puede afirmar que 
cualquier programa requiere seguir el 
diagrama de bloque indicado en la figu- 


ra 14,1. 


Los bloques que forman parte del es- 
cuema de la figura 14.1 toman un sentido 
concreto durante la realización de una u 


otra tarea de una rama úe la mecánica. 


Por sistema fundamental o sistema base 
se comprende en el método de los elcuen- 
tos finitos al conjunto de elementos fi- 


nitos par el cual se sustituye una región 


| CONFORMACIÓN Y ENTRADA DE LOS DATOS DE PARTIDA 


2 FORMACIÓN DEL VECTOR DE LA PARTE DERECHA (CARGAS) 


3 FORMACION DE LAS MATRICES DE LOS ELEMENTOS SEPARADOS 


4 FORMACIÓN DE LAS MATRICES DE TODO EL SISTEMA 


¿ SOLUCIÓN DEL SISTEMA DE ECUACIONES Y DETERMINACIÓN , 
DE LAS CARACTERÍSTICAS NODALES DESCONOCIDAS (INCOGNITAS) 


CALCULO DE CARACTERISTICAS COMPLEMENTARIAS EXIGIDAS 


EN CADA TAREA CONCRETA 


7SELECCION Y SALIDA DE LOS RESULTADOS 


Figura 14.1 


dada. El ejemplo de un sistema base de 
una determinada región se muestra en la 


figura 14.2. 


Cada nudo y cada elemento de la región 
requieren estar numerados. Este esquema 
deberá ser introducido en la máquina a 
través de algún medio y debido a que la 
entrada de datos en forma gráfica a la 
máquina hasta el presente no ha tenido 
una amplia aplicación, será necesario cam- 
biar el esquema por un cierto conjunto de 


números que sea capaz de caracterizarlo. 


En la solución de cualquier tarea por 
el método de los elementos finitos se 
tiene una cierta cantidad de información 
mínima la cual es imprescindible de dar 
a la máquina para lograr la descripción 
y realización de la misma. A este con- 
junto es posible relacionar las siguien- 


tes características en nuestro caso 


(tabla 14.1). 
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TABIA 14.1 


o 5 PP 


Número total de nudos M1, 
Número total de elementos 


Número de puntos en los cuales 
vienen dados los valores de la tel 


peratura M3. 


Descripción topológica de la malla 


de la región A. 


Descripción de las coordenadas de 


los nudos Xo 


Conjunto de los coeficientes de 


transmisión de los elementos Le 


Conjunto de los níímeros de los 

nudos en los cuales es conocido el 
valor de la temperatura La 
Conjunto de los valores de las tenr- 


peraturas en los bordes TH. 


Conjunto de las fuentes distribui- 


das de calor Q. 


— e 5 5 5 5 5 == 


Además, en dependencia del esquema 
concreto de realización del algoritmo del 
método de los elementos finitos que se 
utilice por el calculista, se introducen 
una serie de características complementa- 
rias y de información, tanto con el obje- 
tivo de acelerar el trabajo de cualquier 
parte del prograna, como con el de mejorar 
su organización y disminuir el volumen 


total de los datos. 


En determinados casos son posibles una 
serie de simplificaciones, obtenidas por 
ejemplo, a través de una configuración 
simple o regular de las regiones de cálcu- 
lo o por hacer constantes los valores de 


las constantes físicas. 


Analicemos las etapas fundamentales de 
preparación de la información en el ejen- 
plo de la figura 3.4. Supongaros que en 
el korde superior entre los nudos 1 y 27 
existe una temperatura dada con un valor 
de 5C, y en la frontera inferior (nudos 
del 4 al 32) la temperatura es de 8 Cc. 
El coeficiente de transmisión del calor 
en el caso general es propio para cada 
elemento. No hay fuentes distribuidas de 
calor. Además, a partir de la figura es 
evidente que los valores de las constan- 
tes, serán los siguientes: M1 = 32, 


M2 = 42, M3 = 14, 


Para describir la malla en que se di- 
vidirá la región es necesario describir 
cada uno de sus elementos. El conjunto 4 
tiene dimensiones 1; 3,1: M2 ropredane 
ta una tríada de níímeros en la forma 

(1,6,7), (2,1,7), (2,7,8) etc., donde 
cada uno de los números representa el nú- 
mero de los nudos del elemento de acuerdo 


con el orden de dichos elementos. 
El conjunto X[1:2,1: M1] representa 


los pares de números (141), (u2,42). 


(132,432) cada uno de los cuales represen- 
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ta las coordenadas de los nudos de la 
malla obtenida. De forma similar los res- 
tantes conjuntos tienen la siguiente 
forma: L[1: M2] => (A, ¡A, ).. 242); 

7[1: M3] + (1,6,5,10,15,15,21,27,4,9,14, 
20,26,32); TH 1: M3] > (5,5,5,5,5,5,5,8, 
8,8,8,8,8,8); QA 1: 1] > (0,0 ... 0) 


(42 veces). 


De esta forma nosotres podemos descri- 
bir por completo todos los datos necesa- 
rios del sistema hase. El procedimiento 
utilizado aquí es el más elemental para 
la formulación de los datos. Existen 
otros procedimientos mucho más efectivos 
para conformar los mismos, sin embargo, 
dentro de los objetivos docentes nos con- 


cretaremos al método aquí descrito. 


El siguiente paso es el relacionádo 
con la formación del vector F (bloque 2 
de la figura 14.1), el cual puede ser des- 


crito por el siguiente procedimiento: 


Procedure FORMF; tegin Q1 = QM; 
for X= N1,N2,V3 do F[K ]= £|X] + Qu 


delta/3 end. (14.1) 


Aquí Q1 es un parámetro auxiliar; 
N1,N2,N3 representa los números de los 
nudos del elemento; X, variable entera; 
delta, es área del elemento [ver fórmu- 


la (6.3)] /: número de orden del elemento. 


Pasando a la descripción de la forma- 
ción de las matrices, nos detendremos en 
el problema de realizar una correcta nu- 
meración de los nudos en la malla de ele- 
mentos finitos. Analicemos la estructura 
general de la matriz rigidez del siste- 
ma (MRS) sobre la base del ejemplo de una 
serie de elementos mostrados en la figu- 


ra 14.3 a) y bn”, , 


En el caso de una numeración arhbitra- 
ria por completo | figura 14.3 a)] la co- 


rrespondiente matriz rigidez del sistema 


c) e) 
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Figura 14.3 


[ figura (14.3 b)] resulta en una forma 
completa y lleva desde la primera hasta 
la décima fila y la existencia de ceros 
en la matriz se produce de una manera 
muy irregular. Enel segundo caso de nu- 
meración [figura 14.3 c)] se obtiene una 
matriz rigidez del sistema en forma orde- 
nada [figura 14.3 d)] donde el ancho de 
la media banda (que es la diferencia 
entre el número del elemento no nulo más 
alejado del elemento diagonal y el núme- 
ro del elemento diagonal de una misma 
fila) es igual a seis. En el tercer caso 


[figuras 14.3 e),14,3 £)]el ancho de la 


banda será igual a tres. 


Como es evidente la última forma de 
numeración resulta la más ventajosa por 
cuanto podemos establecer en la máquina 
computadora (IBM) una matriz rectangular 
HH 1: 2xd, 1:M1] donde d es el ancho má- 
ximo de semibanda para la conservación 
de la matriz rigidez del sistema sin 
tomar en cuenta la mayoría de los ceros 
que se encuentran por fuera de la banda. 


Por ello durante la numeración de la re- 
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gión es necesario esforzarse por dismi- 
nuir el valor de d, lo cual permite en 
la memoria operativa de la máquina reali- 


zar un mayor número de operaciones. 


De esa forma, por ejemplo, es necesa- 
rio numerar la región a partir de su lado 
más pequeño | figura 14.3 d)]. La magni- 
tud d se determina por la diferencia má- 
xima entre dos nudos vecinos (para la 
figura 14.2, d = 7) y esta magnituc puede 
introducirse a la máquina conjuntamente 


con los restantes datos. 


En el presente párrafo, con un fin ne- 
tamente docente analizaremos la forma de 
la matriz llena, para la matriz rigidez 
del sistema, los propios lectores pueden 
confeccionar variantes del programa te- 
niendo en cuenta un ancho finito de la 


banda . 


La formación de la matriz rigidez 
(6.2) se describe por el siguiente proce- 
dimiento: : 
Procedure ME; for k <= 1273 00” LoOX 
n = 1,2,3 do h[k,n] = lambda x (bl X]xbl n] + 
+ c[X] xd nl1)/4x% delta). (14,2) 


Tenemos que: 
hl 1:3, 1:3] en la matriz rigidez de 


elemento; Lambda = | NM]; k,n son varia- 


bles enteras. 


El tránsito de la matriz rigidez del 
elemento a la matriz rigidez del sistema 


se da por el siguiente procedimiento: 


Procedure M5; beyin X1 = K2= 0; 
fox K = N1,02,3 do 
begin K1 = K1 + 1; for n = N1,N2,N3 
do begin X2 =X2 +1; il k,n] = 

= H[k,n] + HX1,K2) 


end; K2 = 0 encdl end. (14.3) 


En este caso X1 y K2 son variables 
auxiliares y 4[ 1:41, 1:41] es la matriz 


rigidez del sistema, 


Después de formarse la matriz rigidez 
del sistema es necesario satisfacer las 
condiciones de bordes. Presentemos el 
procedimiento que sirve para realizar 


esta operación. 


Procedure GU, begin for K = 1 step1 
until “43 do begin 1 =71' [X]; form = 1 
step 1 until M1 do An,m) = 0 Hn,n] = 1 
end; for K= 1 step 1 until Mi do 
begin 85 = 0 for n = 1 step 1 

until M3.do bácia M += Anal, 55% só 4 

+ AMK,m]XT4[n]; HXk,m =0 


end, FIX] = PLA] - $8 end end.  (14,4)' 


Para la solución del sistéma defini- 
tivo de ecuaciones utilizaremos un pro- 
cedimiento estandarizado de solución de 
sistemas de ecuaciones algebraicas, al 
cual condicionadamente llamaremos método 
de Gauss. Su uso depende del tipo de pe- 
riféricos que se utilice. Como resulta- 
do se obtiene el vector resultante de la 


temperatura en los nudos Theta (0). 


Presentemos el procedimiento auxiliar 


de formación de las características gyeo- 


métricas del elemento: 


Procedure PR begin M = A11,M; V2 = 


= A[2,N]; N3 = AL3,01 5 X1 =Xx1,01Y; X2 = 
= X(1,02]; X3 = X[ 1,031; Y1=X12,011; Y2 = 
= X[2,02]; Y3 = X12,031; B[ 1] = Y2 - 

= Y3$ B [2] =Y3 - Y1; 513) = Y1 - 2; 
(411) =X3 - X25 (12) = X1 - 43; C[3] = 
X1 - X2; delta = abs (b( 11% 41 + 

B 2]x X2 + B[3luw 43)/2 Lambda = 

LW] end. (14.5) 


+ 


1" 


El programa definitivo de solución de 
la tarea estacionaria de la transmisión 


del calor puede escrisvirse asÍ: 


begin integer M414,012,M03, Ky ¿Rp Dl, 
N1,N2,N3, K1,K2 

real Xx1,42,43, Y1,2,13, delta, lanida, 
98; input (0141,142,113). 

Legin integer array Ál 1;3, 1:42) , L[ 1:43] ; 
array X[ 1:27 1:11) L, QUA], TH] 1:413), 
H 1:49, 1:41], b,C,h11:3), F, Theta 

11: M1]; 

Procedure FOHMF: procedure Mk; Procedu- 
re MS (p); 

procedure GU; procedure Pl; 

input 4,4,7,L,78,%) ; (14.6) 
for K = 1 step 1 until ¿/1 dofor mM 


Ú" 
a 


step 1 until MY do HKk, =0 


! 
O 


for K= 1 step 1 until 1/1 do P[K] 
N = 1; K1: PR FORMF; ME; 45; N = UN + 1; 
if nsdad2 

then got 11; 

GU, 

Gauss; 

Ouput (theta) 


end end. 


Aquí los operadores input y output 
agrupan las correspondientes entrada y sa- 
lida de las magnitudes, indicadas dentro 
de los paréntesis. Su forma concreta .de 
representación deperde del tipo de peri- 
férico. 


El esquema generalizado del programa 
de computación (ver fig. 14.1) es aplica- 
ble también al caso de la tarea no esta- 
cionaria de la transmisión del calor. 

Sin embargo, en este caso, el contenido 
del bloque 5 se modifica, debido a que, y 
a diferencia de la solución de la tarea 
estacionaria obtenemos un sistema de ecua- 
ciones diferenciales relativo a las tem- 


peraturas nodales. 


Analicemos los aspectos de la realiza- 
ción del método de los elementos finitos 
en el caso de una tarea no estacionaria, 
Al igual que en el caso de la tarea esta- 
cionaria, la región investigada se divide 
en elementos triangulares los cuales se 
describen a través de los mismos paráme- 
tros. La regulación del efecto del tiem- 
po se realiza a través de la magnitud del 
paso con relación al tiempo PT y con el 


número límite de pasos Ha 


Además, en lugar del conjunto L se in- 
troduce el conjunto CL [1:2,1:M2] el cual 
representa los pares de números (“, A,), 
le 


2 A,)+- para todos los elementos de la 


región. El campo inicial de temperaturas 
en los nudos de la región se da por el 
conjunto NP [1:41]. En el caso de la 
matriz diagonal de calor específico del 
elemento (“CE) es cómodo representarla en 
forma de un vector unidimensional 

TC [1/41], la formación del cual se des- 


cribe por el procedimiento: 


Procedure MIS for Kk = N1,/2,0M3 do 
1C [XK] = 7C [X] + CL 1,0] < delta/3. 


(14.7) 


Además, en el procedimiento PR(110) es 
necesario cambiar la expresión para el 
coeficiente de transmisión del calor. 

Por completo varía su sentido el procedi- 
miento de satisfacción de las condiciones 


de borde. Es cómodo introducir dos pro- 
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cedimientos: 


Procedure GU1; for K = 1 step 1 until. +13 


do begin n = I[X] for = 1 step 1 until //1 


do il mm] = 0; Mn] = 1 end; 
Procedure (GU2; begin... end. (14.5) 


El procedimiento GU1 produce la anula- 
ción de los elementos de la fila de la 


matriz rigidez del sistema A 1:41,1:41] y 


le adjudica a los elementos de la diagonal 


principal un valor igual a la unidad. 


El procedimiento GU2 determina los va- 
lores de las temperaturas de borde en el 
momento de tiempo Lt. La forma del opexra- 
dor o del conjunto de operadores en este 
procedimiento se determina por el calcu- 
lista de acuerdo con el tipo de tarea a 
resolver. Así, por ejemplo, si en los 
primeros seis puntos que componen el con- 
junto T, la temperatura varía de acuerdo 
con la ley 5sen 10£, y al mismo tiempo en 
los restantes puntos el valor de la tempe- 
ratura es constante e igual a -5 *, en- 
tonces, el procedimiento correspondiente 


toma la forma: 


procedure GU2, begin for A =1 step 
1 until 6 do theta I71XJ] = 5 x sen (10xt) 
7 step 1 until 43 do Theta | 71 4]]= 
(14.9) 


for 4 = 


= TH[7[X1)] end. 


En el último caso se supone que la ten 
peratura constante está dada por el con- 
junto 7H. Los primeros seis números del 
conjunto no tienen valor, por cuanto 
ellos se vuelven a calcular en el proce- 
dimiento (14.9), por lo que, con el fin 
de completar el conjunto a veces es con- 


veniente suponer un valor nulo para estos 


números. 


Una operación fundamental en la tarea 
no estacionaria es la operación de llevar 
la ecuación a la forma de la ecuación 


(7.10). La normalización de la ecuación 


se realiza a través del siguiente procedi- 


miento; 


Procedure NORM; for K = 1 step 1 until M1 
do begin SS = TC [XK] for n = 1 step 1 
until M1 do AlX,n] = H kn] /£S end. 


(14.10) 


La solución del sistema de ecuaciones 
se Yealiza con ayuda del procedimiento 
estandarizado conocido con el nombre de 
método de Runge-Kutta el cual condiciona- 


damente llamaremos RUNGE. 


La forma de operar con el mismo depen- 
derá del tipo de periférico que se utili- 
ce. Entonces la forma general del pro- 


grama será la siguiente: 


begin integer M1,42,M3, R,k,mM,N, N1,N2,N3, 

KV,K2,KR, 

real 41,42,43, Y1,Y2,Y3, delta, Lambda, 

PE, 58, t 

input (411,M2,4M3, PT,R); 

begin integer array A[1:3, 1:M2], 7[1:413,] 

array X[1:2,1:141] CL[ 1:47 1:M3], 

HA 1:M1, 1:41] B,C,h [11:3] NP, TC, 

Theta | 1:41]. 

Procedure ME; Procedure MS, Procedure GUl; 

Procedure GU2; Procedure PR; Procedure UTS; 
Procedure NORM; 

input (4,4,7,CL,TR NP); 
for K = 1 step 1 until 41 do for n= 1 


(14,11) 


step 1 until M1 do Ak,n] = 0 for X= 1 
step 1 until M1 do kegin Theta 

[X] = NP [X]; TC [£] = 0 end; 

WV = 1¿R1:PR; ME; MS; MTS; N =U + 

+ 1: i£ N < M2 Then gota R1; 

GUY; NORM; t = 0; KR = 1; R2: L = + 

+ PT; RUNGE; 

GU2; out put (ft, Theta); KE = KR + 1; 

if Kr <R Then gota *2; end end; KR, con- 


tador de pasos. 


Con ayuda de los programas (14.6) y 


(14.11) se puede resolver una tarea no 
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muy grande (hasta 100-120 nudos) en el 
volumen de memoria operativa de una máqui- 
na computadora de hasta 32 000 células 
(como la MINSK 32, BECM-6 y Otras). Es 
bueno señalar que a pesar de que estos 
programas tienen un carácter meramente 


docente, ellos reflejan por completo las 


'particularidades que surgen en la elabo- 


ración de programas reales. 


El programa para la solución de las ta- 
reas de la teoría de elasticidad plana es 
por su estructura y contenido totalmente 
semejante al programa (14.6). El aumento 
del número de grados de likertad hasta 
seis, aumenta como es natural las dimen- 
siones de las matrices que lo forman. La 
programación de los procedimientos de 
cálculo de las tensiones y las defoxrma- 
ciones por las fórmulas (13.10) y (13.11) 
no presenta ninguna dificultad por lo que 


no nos detendremos en su realización. 


En el presente, a disposición de los 
autores de esta publicación se tiene un 
conjunto de programas (elaborados por 
1.1. Diomin) que dan solución a las ta- 
reas planas y unidimensionales, estacio- 
narias y no estacionarias de la transmi- 
sión del calor, de la teoría de la elas- 
ticidad y de la termcelasticidad. La 
potencia de trabajo de los programas es 
de 3 000 a 4 000 ecuaciones, pero con una 
serie de modificaciones permite resolver 


sistemas con 8 000 ú 10 000 ecuaciones. 


15. EJEMPLOS DE CÁLCULO 


Analicemos los aspectos de precisión y 
de errores del método de los elementos 
finitos en una serie de ejemplos que po- 


seen una solución exacta. 


A. Calentando un medio cilindro homogéneo 


| figura 15.1 a) y b)] dentro del cual está 


Ms e O A 


dada una temperatura constante 
9, = 100*C, y en el exterior 
0, =0”C . La solución exacta de 
esta tarea se presenta en este tra- 


bajos. * 


Figura 15.1 


pre 
Ny 


Teniendo en cuenta la simetría de la 
Yegión analizaremos solamente la cuar- 
ta parte del anillo, la cual es repre- 
sentativa del anillo en su conjunto. 
FL error entre la solución exacta y la 


* solución numérica no excede del 0,05 %. 


B. Enfriamiento de una banda seminfinita 


[fig. 15.2 a) y Bl 


Figura 15.2 


La solución exacta de esta tarea se 
brinda en el epígrafe 15. La solución 
numérica solo en la etapa inicial difiere 
algo de la solución exacta pero al cabo 


de un período de tiempo (en nuestro caso 


b) 


— SOLUCION EXACTA 
000 SOLUCION MEF 


3 h) la misma coincide con la solución 
exacta. Del gráfico [ fig. 15.b)] es evi- 
dente que tanto en el caso de la matriz 
completa; como en el caso de la matriz 


diagonal se obtienen buenos resultados. 


Durante la utilización del método de 
los elementos finitos a la función de la 


forma se aplican una serie de restriccio- 


b) 


8 grados t:5h 


Solución exacta 

o Solución MEF con matriz 
llena [Cc] 

Ñ solucion MEF con matriz 
diagonal [Cc] 


2 JA e e Bo 9 O: (mM 


Figura 15.3 


nes por la selección de tipos específicos 
de funciones de interpolación, por esta 
razón el método de elementos finitos en 
general da una cierta reducción de los 
valores de la temperatura cuando se uti- 
lizan las matrices llenas de rigidez del 
elemento, lo que se ratifica por el grá- 
fico de la figura 15.2 kb) para un tiempo 
igual a2h ($ = 2 h). 


En este caso obtenemos la frontera in- 


ferior de la solución. La utilización de 


matrices [c1), en forma diagonal en algo 
suaviza la "rigidez" general del método 
de elementos finitos y permite obtener la 
frontera superior de la solución. Duran- 
te la continuación del proceso ambas so- 
luciones convergen rápidamente a la solu- 
ción exacta y en los límites del error 


adoptado coincide con ella, 


a) 


C. Distribución de la temperatura en una 
capa de suelo [fig. 15.3 a) y b)]. 


Se tiene una capa de suelo de longitud 
infinita y espesor igual a 1 MT [ figu- 
ra 15.3 a)] en la superficie de la cual 
se aplica una variación diaria de la tenm- 
peratura a través de una fuerza armónica 
del tipo E =0 


cp + A, sen wW (t + tp). 


donde: 


0. : Temperatura media diaria. 


A.: Amplitud, 


w: Frecuencia, 


dis 


y + Fase inicial de las desviaciones. 


Investiguemos el cambio del régimen de 
temperatura en suelo bajo los siguientes 
valores de los parámetros: ap = -10 *C; 
A, =20"C; W= 0,262 I/Mi ty =7,8h. 


b 
9 = Op + Agsen wit + to) 


e(x,y,0) =07 


lm 
6h 12 h mh 
-0,92 7,82 -19,08 
AER 

Az E 
E 
ME] 
== 


ATT) 


o 
e ana dl 


y comparamos los resultados con la solu- 
ción exacta obtenida por M.G. Vaniuchen-= 
kovi y NaN. Malikovi por el método de las 


funciones de origen. 


0 (uyrtd = U- 0d a +70, + 


+ LN sen wW(t + t,) + %, cos w(t + ty) + 


> na, sen NN 
+ ZA E AA 
1 n exp ( pP ) MN 
(15.1) 
donde: 
b, = A, cos k pa e, d, = 
—k Bo 
= A, sen k ne a RS E 
dy BEI y + 2 op 
atra, a nm? A, 
- sen wWt + 20 e 
a+ nin? at 
4 cuando n = 1,3,5 + 
Pa 
0 cuando n = 2,4,6 -- 
2 cuando n = 1,3,5 + 
EA 
2 cuando 1 = 2,4,6 + 


En la figura 15.3 b) se presentan los 
gráficos del campo de temperatura en es- 
tado casi estacionario. En este Caso, 
también el método de los elementos fini- 
tos da resultados que coinciden perfecta- 
mente con los resultados exactos. La 
desviación de la solución numérica a la 
exacta no excede de unas cuantás centési- 


mas de grado. 


Al igual que en el caso anterior la 
utilización de la matriz [C] implica des- 


viaciones en la media de alrededor de 
0,02 ee por debajo de su valor, mientras 
que la matriz LC], en forma diagonal da 
esa misma desviación, pero por exceso. 

El valor medio entre ellas prácticamente 
coincide exactamente con la curva de la 
solución analítica exacta. Con un aunen- 
to de la densidad de la malla, se puede, 


inclusive, eliminar este pequeño error. 


En la figura 15.4 se muestra el carác- 
ter de convergencia de la solución al ré- 
gimen cuasi estacionario. En este caso, 
po PON es el valor cuasi estacionario de la 
temperatura a una profundidad dada. Las 
diferencias entre las soluciones exactas 
y numéricas se observan solo en los pri- 
meros cuatro Ú cinco días, después de 
los cuales ambas soluciones coinciden 


totalmente. 


Si el campo inicial en la capa tiene 
temperatura de 10 %c, entonces, el esta- 
do cuasi estacionario aparece algo des- 
pués de los cinco primeros días, sin em- 
bargo, la forma de distribución del 


mismo no varía (ver fig. 15.4)» 


De los ejemplos presentados es eviden- 
te que el método de los elementos fini- 
tos da resultados suficientemente preci- 
sos durante la investigación de la carac- 
terística fundamental del proceso (en 


este caso, de la temperatura) .?* 


Analicemos algunos ejemplos donde una 
de las tareas fundamentales sea el cálcu- 


lo de los vectores o y € 
D. Compresión de una capa de suelo por 
la acción de su peso propio. 


Analicemos una capa de suelo 
[fig, 15.5 a) y b)], cuyas característi- 


cas son: 


E = 35,0 MPaz Y = 0,35 y P = 17 640 N/A. 


100 (?/a) 


—-SOLUCIÓN EXACTA 

O SOLUCIÓN DEL M.E.F CON MATRIZ 
"LLENA € 

6 SOLUCIÓN DEL MEF CON MATRIZ 
DIAGONAL Cy 
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Dentro de la capa se cumplen las condi- 
ciones del estado deformacional plano. 
La región de cálculo se divide en 750 
elementos. La solución exacta de esta 
tarea se determina por las conocidas fór- 


mulas de teoría de la elasticidad! *. 


¿od Me 


y Pi A yy (15.2) 


Del análisis de los gráficos de las 
tensiones 0, Y e [ figura 15.5 1)] se de- 
termina que el método de los elementos 
finitos da un error considerable en el 
extremo libre debido a que las fórmulas 


(13.11) permiten determinar solamente el 


valor medio en el interior del elemento. . 
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tídíos) 


Esta es la insuficiencia fundamental del 
método de los elementos finitos en forma 


de método de los desplazamientos. 


En la frontera inferior el error al- 
canza aproximadamente un valor del 2%. 
En el resto de la región las tensiones 
aplicadas en los nudos de acuerdo con la 
fórmula (13.12) coinciden con los valores 
exactos. A partir de ello, es evidente 
que aún en los casos más sencillos duran- 
te la solución de tareas prácticas, los 
elementos vinculados directamente a la 
frontera libre del elemento es necesario 
sacarlos del análisis y determinar el 
tipo de campo de tensiones por la si- 


guiente fila de elementos. 


Figura 15.5 


D. Estado tensional de un cilindro calen- 


tado (figura 15.6).. 


La solución exacta se obtiene en los 
trabajos! ?*%, Las características fÍ- 
sicas: 

a = 50 em, b = 100 cm, £ = 10MX a, 

Y = 0,3, 4 = 10"*1/grad. La tarea se re- 
suelve para la mitad del anillo (teniendo 
en cuenta las características de simetría 
del problema)» En estas condiciones, Se 
establecen las características de sime- 
tría Vaz 
desplazamientos verticales del cilindro 


= UY, = 0. Para excluir los 
co 


como cuerpo rígido se estableció la con- 
dición Vg = 0. De los gráficos [ver fi- 
gura 15.6 b) y c)] es evidente que la so- 


» 
— solucion exacta 


O Solución del métudo de elementos finitos 


lución por el método de los elementos fi- 
nitos de la tarea de termoelasticidad 
coincide adecuadamente con la solución 
exacta, en las regiones que se encuentran 
separadas del borde libre. En el Lorde 
libre la solución del método de los ele- 


mentos finitos no es confiable. 


De los ejemplos analizados se puede 
concluir que el método de los elementos 
finitos de una buena precisión en los va- 
lores de la característica: fundamental 
buscada: temperatura o desplazamientos en 
los casos analizados. El cálculo de la 
media aritmética de estos valores en los 
nudos permite obtener resultados que po- 
seen una precisión suficiente para los 


.fines prácticos. 


— SOLUCIÓN EXACTA 


(3 
O SOLUCION POR EL 
Y 
METODO DE LOS ELEMENTOS 
FINITOS 


Figura 15.6 
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"NOTA DEL TRADUCTOR 


La traducción de la monografía "El mé- 
todo de los elementos finitos en las apli- 
caciones a la teoría de las construccio- 
nes" de los profesores N. N. leontiev e 
T. I. Diomin constituye un paso más en el 
propósito de presentar al estudiante cu- 
hano de la especialidad de Ingeniería 
Civil los métodos modernos de cálculo de 
estructuras a través de los libros de 


textos de la escuela soviética. 


La presente obra estimula a su vez el 
interés por el conocimiento de la lengua 
rusa, imprescindible (al igual la lengua 
inglesa) para todo especialista que as- 


pire a mantenerse medianamente informado 
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de los logros de la ciencia y la técnica 
en la rama del conocimiento a la cual se 


dedica, en este caso, la Ingeniería Civil. 


Las pretensiones de la publicación, 
claramente expresada por los autores en 
su introducción son limitadas, pero al 
mismo tiempo, suficientes para informar 
e iniciar al lector en el conocimiento 
del método de los elementos finitos debi- 
do a que orienta sobre posibles fuentes 
bibliográficas a las cuales recurrir en 


caso de desear profundizar en su estudio. 
4 


Además, se brinda una exposición sobre 
las características de las soluciones 
analíticas y numéricas, así como de los 
modelos físico y matemático a que se re- 
duce cualquier problema de análisis es- 
tructural, que a pesar de su sencillez, 
abre horizontes, y ayuda a ordenar con- 
ceptos generales, necesarios para todo es- 


pecialista en su etapa de formación, 


Esperamos que esta traducción, reali- 
zada como una tarea extraplán, en el 
marco de la Jornada Ideológica Camilo-Che, 
sirva en tus manos para inculcar un mayor 
número de conocimientos, y ayude en su 
propósito de estar cada día más como se- 
nalara nuestro Comandante en Jefe el pa- 
sado 26 de Julio "a la altura del reto 
planteado entre el socialismo y el capi- 
talismo", convencidos de que el futuro, 


pertenece por entero al Socialismo. 


CÁLCULO DE CÁSCARAS PRISMÁTICAS 


1. INTRODUCCIÓN 


En este segundo tema se exponen los 
fundamentos de la mecánica de la cons- 
trucción de los sistemas espaciales de 


paredes delgadas. 


Se entiende como tal aquellas estruc- 
turas que están compuestas por placas 
rectangulares alargadas, unidas entre sí 
de Forma continua O por articulaciones a 


lo largo de sus lados mayores. 


iste tipo de sistema tiene una utili- 
zación tan extendida en diferentes ramas 
de la técnica que si pretendiéramos hacer 
un recuento completo de las mismas, ocu- 
paría una extensión mucho mayor que la 
que podemos permitirnos en esta monogra- 
fía, por lo que nos limitaremos a señalar 


solamente algunas de sus aplicaciones. 


Comenzaremos con un breve resumen de 
ellas a los problemas de la construcción, 
teniendo en cuenta que esta monografía 
está dirigida en primer lugar a estudian- 
tes de las especialidades de Ingeniería 
Civil y Arquitectura aunque puede ser 
útil también para la preparación de espe- 


cialistas de otras ramas técnicas. 


En este caso, los sistemas espaciales 
de paredes delgadas las encontramos, 
“sobre todo, formando parte de cubiertas 
plegadas de edificaciones civiles e in- 
dustriales | fig. 1.1 a)] para el cálculo 
de las cuales es que se desarrolla la 
teoría que aquí presentamos. Además, po- 


“¡emos encontrarlos en las llamadas cu- 
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TEMA 2 


biertas dentadas o aserradas | fig. 1.1 b))]. 
en los esquemas de cálculo de los balco- 
nes de teatros, cines y salones de confe- 
rencias [fig. 1.1 c)], los cuales repre- 
sentan estructuras plegadas. No es difí- 
cil darse cuenta también que mechas 
edificaciones industriales modernas de 
una sola planta, representan sistemas es- 
paciales de paredes delgadas, compuestas 
de placas rectangulares | fig. 1.1 d), e), 
Dn]. 


serie de estructuras de edificios altos. 


Esto es aplicable también a una 


Los paneles de los entrepisos, en muchos 
casos, pueden ser calculados con un alto 
grado de precisión como sistemas espacia- 
les de paredes delgadas | fig. 1.1 9)]. 
Además, a este grupo pueden adicionarse 
ferrocarriles 


los puentes de carreteras, 


y los viaductos | fig. 1.1 h)]. 


Como ejemplos de sistemas espaciales 
de paredes delgadas entre las construc- 
ciones hidrotécnicas se pueden señalar 
los canales de desague de las presas, sus 
vertederos, sus compuertas y los edifi- 
cios de las centrales hidroeléctricas, 


entre otros, 


En la construcción de barcos frecuen- 
temente se utilizan diques flotantes los 
cuales pueden ser calculados por esta 
teoría [ fig. 1.1 1)], además, es necesa- 


rio señalar que muchos barcos modernos 


“pueden analizarse como sistemas espacia- 


les de paredes delgadas. 


Muchos elementos de los aviones, como 
por ejemplo, las alas, el fuselaje, etc. 
| fig. 1.1 j)] se calculan en el presente 
y se proyectan como sistemas espaciales 


de paredes delgadas. 


La mayoría de los medios modernos de 
transporte vagones metálicos, del trans- 
porte ferroviario, de los metros, de tran- 
wías de alta velocidad, de la carrocería 
de los ómnibus, de los trolebuses, etC.; 


pueden ser analizados como sistemas espa- 


ciales de paredes delgadas. 


Figura 1.1 


Este pequeño y limitado resumen permi- 
te al lector convencerse de que los sis- 
temas espaciales de paredes delgadas se 
utilizan ampliamente en el presente y sin 
duda alguna, estos alcanzarán una mayor 


extensión en el futuro cercano. 
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2, HIPÓTESIS FUNDAMENTALES, - DETERMINACIÓN 
DEL NÚMERO DE FUNCIONES INCOGNITAS 
Y DESPLAZAMIENTOS POSIBLES 
DE UN PÓRTICO ELEMENTAL, FÓRMULAS 
DE TENSIONES Y SOLICITACIONES 
FUNDAMENTALES 


En la hase de la teoría técnica en ge- 
neral de las cáscaras se encuentran las 
conocidas hipótesis de Love-Kirchoff, las 
cuales constituyen una extensión al caso 
bidimensional de las conocidas hipótesis 
de la teoría elemental de la flexión de 
las vigas que se explica en los cursos de 


resistencia de materiales. 


Primera hipótesis: El elemento normal 
a la superficie media m-n (fig. 2.1) des- 
pués de la deformación de la cáscara, se 
mantiene recto y perpendicular a dicha 


superficie conservando su longitud. 


Figura 2.1 


Segunda hipótesis: Las capas elementa- 
les aisladas que conforman el elemento 
en la dirección normal a su superficie no 


presionan unas sobre otras. 


Para los sistemas prismáticos de pare- 
des delgadas analizados, los cuales están 


compuestos de placas rectangulares alar- 


gadas (l/a > 5), rígida o articuladamente 
ligadas entre sí, fueron introducidas hi- 
pótesis adicionales por V.Z. Vlasov.! >? 


La primera hipótesis se refiere a que 
después de la deformación del sistema el 
ancho de la placa no varía y la proyec- 
ción de la línea A-B [ fig. 2.2] que co- 
rresponde a la posición Z = constante, 
sobre el plano de la placa, se mantiene 


recta. / 


Esta simplificación es equivalente a 
la conocida hipótesis de las secciones 
planas, aplicada en la Teoría elemental 
de la flexión de las vigas y que exitosa- 
mente ha sido aplicada en la práctica de 


la ingeniería. 


Figura 2.2 


La placa en este caso se analiza como 


una viga, la cual se flexa en el plano de 


su mayor rigidez (sección transversal), 
pero con una diferencia esencial, mien- 
tras en la teoría elemental de la flexión 
de las vigas se desprecia la deformación 
por cizallamiento, en el análisis de los 
sistemas espaciales de paredes delgadas 


esta deformación se toma en cuenta. 


La segunda hipótesis está relacionada 


con el hecho de que las tensiones norma- 


les y tangenciales en las secciones norma- 
les se distribuyen a través del espesor 
de forma uniforme, lo que lleva a la dis- 
tribución de fuerzas interiores que se 


muestra en la figura 2.1. 


Esta hipótesis es clara, lógicamente, si 
se representa la flexión de una placa 
rectangular alargada, empotrada por su 
contorno y sometida a una carga que varía 
lentamente a lo largo de la coordenada 2 


(fig. 2.1). 


En ese caso la curvatura de la flexión 
en la dirección longitudinal y la curva- 
tura de la torsión son considerablemente 
menores que la curvatura de la flexión 
del contorno de la sección transversal, 
lo Cual justifica el hecho de que se des- 
precien los momentos flectores en la di- 
rección longitudinal y los momentos tor- 
sores en comparación con los momentos 


flectores en la sección transversal. 


Estas hipótesis permiten analizar el 
sistema prismático, como compuesto por un 


conjunto infinito de pórticos planos. 


Los desplazamientos de cada vórtico 
plano en dirección perpendicular a su 
plano, quedan determinados por completo a 
partir de los desplazamientos de los nudos 
del pórtico a lo largo del eje (0-2) (ver 
fig. 2.1). El grado de libertad de los 
nudos del pórtico en dirección perpendi- 


cular a su plano es igual a UN) donde: 
n= y.a 251) 
donde 


y: Número de nudos del pórtico inclu- 


yendo los de apoyos» 


C,: Número de ligaduras o barras de 
apoyo que se encuentran fuera del 


plano del pórtico. 


El número de desplazamientos lineales 


de los nudos del pórtico en su plano se 


detexminan por la fórmula: 
(2.2) 


¿onde 


CU: Número de barras del pórtico y de 
barras de apoyo que descansan en el 


plano del pórtico. 


Es necesario aquí destacar que debido 
al hecho de que el sistema espacial está 
compuesto de un núnero infinito de pórti- 
cos planos, las incógnitas en este Caso 
no serán números, sino serán funciones y 
por ello las ecuaciones canónicas del mé- 
todo de los desplazamientos no resultarán 
ecuaciones algebraicas, sino ecuaciones 


diferenciales ordinarias. 


Analicemos detalladamente los posibles 
desplazamientos del pórtico elemental 
plano que responde a las hipótesis adop- 


tadas. 


En el caso de los desplazamientos del 
pórtico plano en la dirección perpendicu- 
lar a su plano, el mismo se desliza a 
través de la superficie prismática del 
sistema. Sus desplazamientos se determi- 


nan por las funciones AS una de las 


cuales se muestra en la figura 2.3. 
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Los desplazamientos entre los puntos 
nodales se determinarán a través de una 
ley lineal en correspondencia con la pri- 


mera hipótesis de V.Z. Vlasov. 


La magnitud real de los desplazamientos 


longitudinales de los nudos de diferentes 


pórticos la lo largo del eje 0-2) es tam- 
bién diferente por lo que los desplaza- 

mientos longitudinales de cualquier punto 
del sistema prismático se puede represen- 


tar en la forma: 


Utz,0) =U, (e) 0, (8) +0, (2) 0, (8) + == 


n 
¿2 Ur a) 
donde: 
U, (a): Funciones buscadas de una sola 
variable. 
Para comprender lo que sigue, es con- 


veniente destacar que en el caso del des- 


plazamiento analizado (fig. 2.3) surgen 


deformaciones por cizallamiento las cua- 


les son iguales a: 


we 
= y : a 
> Y (s) (2.4) 


Llamemos también la atención del lec- 
tor sobre un heclo totalmente evidente a 
partir de los gráficos de y, (8): y es el 
hecho de que se pueden combinar nuevos 
grupos de gráficos linealmente indepen- 
dientes, como es el caso de un gráfico 


simétrico y otro antisimétrico. 


El dibujo de las deformaciones del pór- 
tico en el caso de un desplazamiento ele- 
mental en su propio plano se muestra en 


la figura 2.4. 


Sus desplazamientos en el contorno 
(las proyecciones del vector desplazamien- 
to sobre el contorno del pórtico) se re- 


presentan por la función Y E). 


Det 


- a 


A AAA > - 
3 1 
4 


Figura 2.4 


El gráfico M, (8) de los momentos flec- 
tores que responden a uno de estos esta- 
dos es cómodo construir por los métodos 
de la mecánica de la construcción de los 
sistemas de barras planas estáticamente 
indeterminadas, por el método de los des- 
blazamientos, por el método del equilibrio 
de los momentos nodales, de Bernarski- 
Cross y otros, Rara vez resulta conve- 


niente usar el método de las fuerzas, 


De la misma forma en que los desplaza- 
mientos de los diferentes pórticos que 
componen el sistema son diferentes, enton- 
ces, los desplazamientos del contorno de 
cualquier punto del sistema espacial pue- 


den sex representados en la forma: 


Viz,6) = 7, (2) Y, (8) + V, (3) Y, (8) + 
Mm 
t 5 EV y EW, (8) (2.5) 


donde : 


Y, (E): Funciones buscadas de la varia- 


ble 2, 


Utilizando las conocidas relaciones 
del caso plano de la teoría de la elasti- 
cidad y de la mecánica de la construcción 
de los sistemas de barras y también las 
ecuaciones (2,3) y (2,5) para expresar 
los desplazamientos, obtenemos las si- 
guientes fórmulas para las tensiones y 


solicitaciones fundamentales: 


n 
ou a 
0(2,8) = E pr E ¿E U, (2) es (8) 
(2.6) 
du d 
r(2,8) =.G 7 + =— ) = 
a, Ya 


GH ¿2,0, (2) p;(s) + Y o 0, (8)] 
(2.7) 


Un 


M(z,8) 


My (248) + ES Vz (2) bi, (s) 


(2.8) 


Aquí, E y G son los módulos de elasti- 
cidad en tracción y cortante respectiva- 
mente; 4, (48) son los gráficos de momen- 


tos flectores de la carga. 


La fórmula (2.6) expresa la ley de 
Hooke para un estado tensional unidimen- 


sional. 


El asterisco sobre las funciones Var 
Vos qa Y, representa derivación con 
respecto a la misma variable de la cual 


depende la correspondiente función. 


3, ECUACIONES DIFERENCIALES DE EQUILIBRIO 
EN FUNCIÓN DE LOS DESPLAZAMIENTOS 
Y ESTABLECIMIENTO DE LAS CONDICIONES 
DE BORDE 
Analicemos el equilibrio de un pórtico 
elemental (fig. 3.1) extraído del sistema 
espacial por dos secciones infinitamente 


cercanas. 


Dibujemos todas las fuerzas exteriores 
en relación con el elemento dado. Al 
elemento ds de la sección transversal 
% = constante, apliquemos los esfuerzos 
oóds y T6as, 


donde: 
0: Espesor de la placa dada, 


Mientras al elemento ds de la sección 
Z + d 3 = constante aplicamos esfuerzos 
semejantes, pero teniendo en cuenta el 
incremento de los mismos, justificado por 


el aumento de la coordenada 2. 


Además, el elemento dzd8 de la super- 
ficie media está aplicada la carga exte- 
rior dada, representada por las componen- 


tes q y Pe 


Para el establecimiento de las ecuacio- 
nes de equilibrio del pórtico dado utili- 
zaremos el principio de los desplazamien- 


tos posibles de Lagrange. 


De acuerdo con ese principio para el 
sistema que se encuentra en equilibrio el 
trabajo de todas las fuerzas en un des- 
plazamiento posible cualquiera tiene que 


ser igual a cero, obteniendo : 


$. === 02:ddes y 


0 - 5rb dado 3 + 


$8 


+ $ pdzds Y, =0 (3.1) 
6 


dad 


J e Tiida Y —¿— deds + 


+$ qdzds y, = 0 (3.2) 
g 


Cualquiera de las ecuaciones (3.1) ex- 
presa la igualdad a cero del trabajo de 
las fuerzas exteriores e interiores sobre 
los desplazamientos en la dirección per- 
pendicular a su plano, determinado por la 


función O" 


Las ecuaciones (3.2) tienen un sentido 
similar, pero los desplazamientos, sobre 
los cuales se calcula el trabajo, ocurren 


en el plano de la seceión transversal y 


«se determinan por la función Yi 


El sentido de los primeros y últimos 
sumandos en la parte izquierda de las 


ecuaciones (3.1) y (3.2) es evidente. 


El sumando central en la ecuación 
(3.1) representa el trabajo de la ten- 
sión tangencial interior T en relación 
con el pórtico elemental sobre la defor- 


mación de cizallamiento Pr 


El sumando central en la ecuación (3.2) 
determina el trabajo de los mlmomentos 


flectores transversales Y sobre los ángu- 
4, 
los de giro mutuos El ds. 


El signo menos delante del sumando 
central se debe a que la fuerza interior 
del sistema elástico se opone a la defor- 


Dd 
mación, 
Sustituyendo en las ecuaciones (3.1) 


y (3.2) las expresiones (2,6) y (2,8) y 


dividiendo sobre dá, obtenemos: 


UN UN Na 
y Basile AU 5 A Lp = 
11 tj t 131 Cit PA kg K Ge] 
Un % Ta y 1 
OS de 7 Z — q,=0 
AE AR ERRE 
(3.3) 


Los coeficientes de las ecuaciones ob- 


tenidas se determinan por las fórmulas: 


E o , 
Ya... A apo Pu y vv Ade 


Cay 7 E Ye $ds; Y kh e Ye Y, ds: 
MH 
1 Ral e 


8g 


_Y poseen la propiedad de la reciprocidad. 


Los términos independientes se deter- 
minan como el trabajo de la carga dada 
sobre los desplazamientos que responden a 


las funciones P, U Y 


A A A (3.5) 


Destaquemoy, que bajo las hipótesis 


adoptadas las funciones Y, es Ye son li- 


neales, y los espesores 5 son constantes 

dentro de los límites de cada placa, por 

lo tanto, durante el cálculo por las fór- 
mulas (3.4) y (3.5) es posible y útil 


utilizar la regla de Vershagin. 


Durante la obtención de las ecuaciones 
(3.3) fue utilizada la propiedad de oxto- 
gonalidad de los yráficos 4, en relación 
con los gráficos unitarios, la cual se 
conoce de los cursos de mecánica de la 
construcción de los sistemas de barras 


planas. 


dá 
O ads =0 (3.6) 


En el resultado de la integración del 
sistema | ver ecuación (3.3)] de ecuacio- 
nes diferenciales ordinarias (los aspec- 
tos de la integración de este sistema los 
analizaremos más abajo). Entran en la 
solución 2(1, +M,) constantes arbitra- 
rias. Estas constantes arbitrarias se 
determinan a partir de las condiciones 
de bordes en las secciones extremas del 


sistema prismático. 


Se limita en este caso al análisis de 


condiciones de borde simples. 


Comencemos con el caso cuando la sec- 
ción extrema del sistema z = O está uni- 
da con un masivo rígido y se obtiene una 


ligadura totaimente rígida (fig. 3.2). 


Entonces para este borde se puede es- 


cribir: 


U(0,8) = 0 ; V(0,s) = 0 (1857) 


Figura 3.2 


Estas condiciones serán satisfeclias 
si en correspondencia con las fórmulas 


(2.3) y (2.5) nosotros aplicamos: 
U¿(0) = 0  V,t0) =0 (3.8) 


Las condiciones analizadas recilen el 
nombre de condiciones cinemáticas o geo- 


métricas. - 


Pasemos al caso de borde lilre 2 = L, 
que está cargado por las tensiones 
0, 3 0, (1,8) 1 7 (40) (3.9) 


l ¿ 


De acuerdo con las fómmulas (2.6) y 


(2.7) se puede escribir que: 


Ni 
a , m3 % 


UN Ta 


, , Ss 


(3.10) 


Las condiciones (3.10) se pueden sa- 
tisfacer solo en casos exclusivos. La 
exigencia de que las tensiones normales 


o, en los límites de cada placa se dis- 


tribuyan de acuerdo con una ley lineal, 
mientras que las tensiones tangenciales 
son constantes, es una condición necesa- 


ria pero no suficiente. 


Por eso, en el caso general en que se 
¡da ta y va es. necesario comprender que no 
es posible satisfacer de forma precisa 


las condiciones de borde estáticas (3.10). 


Debido a ello hay que exigir que se 
satisfagan las condiciones de igualdad de 
los trabajos de las fuerzas externas e 
internas, aplicadas a la línea del con- 
torno de la sección extrema, en los ante- 
riormente seleccionados desplazamientos 


posibles. 


n 
A wi -. 
S E ¿EV (De (se, (e)Jóde = 


= $ ¿97 (6)e ,(8)6 de 


Ni 


FELIZ. (e: (8) + ES ve ALA (131 V; (e)óds = 


ul L= 


a 
19) 


= $ ¿7 ¿(61V, (8)5 ds 


(3.11) 


Utilizando la nomenclatura (3.4), a 
estas condiciones se le puede dar la forma 


siguiente: 


UN 
EE AO te) 2 0, ds 


mn ma 
je AA > es 2 pl T TY ds 


(3.12) 


Para el extremo libre no cargado es ne- 
cesario poner en la parte derecha de las 


ecuaciones (3.12) que a e O. 


Analicemos, por fin, las condiciones de 
apoyo del extremo del sistema prismático 
sobre un diafragma (fig. 3.3) absolutamen- 
te rígido en su plano, y absolutamente 


flexible en dirección perpendicular a 
dicho plano. 


En este caso se puede escribir que: 


V(0,S) = 0; 0(0,8) = 0 (3.13) 


A partir de las cuales, y después de 
utilizar las fórmulas (2,5) y (2,6) ten- 


dremos: 


V, (0) =0, U¿ (0) = 0 (3,14) 


» Figura 3,3 


listas condiciones de borde (3.13) se 
llaman mixtas, teniendo en cuenta que una 
parte de ellas expresa la ausencia de des- 
plazamientos y la otra la ausencia de ten- 


siones. 


li, INTEGRACIÓN DE LAS ECUACIONES 
DIFERENCIALES POR EL MÉTODO EXACTO 


Para mostrar la metodología de integra- 
ción de los sistemas de ecuaciones obteni- 
dos y hacer las correspondientes conclu- 
siones, es suficiente analizar un sistema 
reducido de ecuaciones. 

Yi, Ul=<b le == Y=0 
11 1 PLA 


dd 
(4.1) 


1 
1Ós AMES” E IN Y EA = 
sl | 241 1 11 11 G q, a 


q, = constante 


El sistema (4.1) de ecuaciones diferen- 
ciales ordinarias es lineal y no homogéneo 


por lo que su solución tiene la forma: 


UA = ga? aa ya, Y = yo) pra ye 
N 


(4.2) 
donde: 
ye, AE Solución general del sis- 
icua homogéneo de ecuacio- 
nes (4.1) cuando (q, = 0), 
ye pe, Solución particular del 


sistema de ecuaciones (4.1), 


El sistema de ecuaciones (4.1) tiene 
coeficientes constantes por eso su solu- 
ción (para dy = 0) tiene la forma: 


rá 


E (4.3) 


Sustituyendo las expresiones (4.3) en 
el sistema (4,1) para q, = 0 y sacando 

a rE A 

factor común e obtenemos el sistema de 


ecuaciones algebraicas homogéneas: 


va Pb JA-C rbB=0 
1L > ES! 11 


(4.4) 
COrYAÁs+ (a Y-5 )JB=0 
11 11 11 


Este sistema tiene dos tipos de solu- 


ciones: 


La primera solución en el caso 
' 
A =B= 0 resulta una solución trivial, 
la cual corresponde como es evidente a 
partir de las fórmulas (4.3) a la ausen- 


cia de desplazamientos del sistema. 


La segunda solución, en la cual 
AFOyYB%*O la obtenemos a partir de la 
condición de hacer cero el determinante 
de la matriz de los coeficientes que mul- 


tiplica al vector de A y Bo 


56. 


o 2 Mm ai = Ss.) 
ED 11 11 
(4.5) 


Desarrollando el determinante (4.5) 


obtenemos ; 


*P -2ar +b=0 (4.6) 
donde : 
2 , y *o"+YE JS Ya 
PI y 11 11 1 Ca ya 
11 11 
Bb: 18 
5 de EE (4.7) 
TE, P 


Para sistemas de paredes delgadas, su- 
ficientemente finos, el coeficiente $,, 
es extremadamente pequeño en comparación 
con todos los demás oveficientes del 
sistema de ecuaciones (4.1). 

Por eso, como regla general en el caso, 
d > b, las rálces de las ecuaciones 
(4,6) son reales y se detemifian de acuer- 


do con las fórmulas: 


r =Ya=vVa -b =* 
152 
(4.8) 
2 
y =tYa+ya -b =*B 
3,14 


. 
Corresponde a estas ralces la solución 


siguiente: 
ye? = Á, cosh az + A, senh az + 
1 1 
+ A, cosh Bz +A, senh Bz 
(9) . 
ja =p,4, senh «ñ +p,4, cosh 02 ++ 


+p,4, senh Bz + P,A, cosh Bz 


(4.9) * 


La forma general de las funciones hi- 
perbólicas cosh < y senh x se muestran 
en la figura 4.1. Ellas se determinan a 
través de la función mostrada por las 


fórmulas: 


cosh x= 


A través de estas fórmulas se pueden 
calcular los valores del cosh x y el 
senh Y para tantos valores del argumento, 
comó no pueden reflejarse en las tablas 


de funciones hiperbólicas. 


Por diferenciación directa de las fór- 


mulas (4,10) no es difícil comprobar que: 


<G—(cosh X)=senh x; sen x)=cosh x 
(4.11) 


La forma adoptada de la solución para 


(o) 


1 se obtuvo sobre la base 


la función V 
de esta representación. En cualesquiera 
de las ecuaciones (4.1) aparece una fun- 
ción con su segunda derivada y la prime- 


ra derivada de la otra función. 


Por eso, teniendo en cuenta las rela- 


ciones (4.11) y lo señalado anteriormen- 


lo) 


1 
escribirla en la forma presentada ante- 


te, la solución para V es necesario 


riormente en la ecuación (4.9). 
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Sustituyendo las expresiones (4.9) en 
las ecuaciones (4.1) para el caso q, = 0 
encontraremos que las mismas pueden sa- 


tisfacerse para valores P ¡, Como: 


va a. -b Co Q 
P mo 11 11 En 11 EZ 
; Cc a Ss -p a 
11 11 11 
va p?-b CB 


(4.12) 


En este caso no es ditiícil convencer- 
se que si el valor de « se calculó exac- 


tamente, entonces ambas fracciones en la 


expresión de P, deben coincidir. 


Eso es consecuencia de haber igualado 
a cero el determinante del steama (4,4) 
compuesto de los coriciente: j1e multi- 


plican al vector de 4 , be 


Lo mismó puede ser expresado en rela- 
. 2 - 2 
ción con las fracciones en la expresión 


para P, con relación ab. 


De esta forme, al igualarse las co- 
rrespondientes fracciones en las expre- 
siones (4.12) se obtiene un método de con- 
trol para la comprobación del cálculo 


de a yb. 


La solución en la forma (4.9) es cómo- 
da en el caso en que se investiguen sis- 

temas prismáticos de longitud finita, en 

particular, si las condiciones de bordes 

en las secciones extremas son equiva- 


lentes. 


En este último caso, la mitad de las 
constantes arbitrarias As debido a la si- 
metría del problema toman valores igua- 
les a cero, comprendiéndose que en este 
caso, el origen de coordenadas se encuen- 
tra situado a la mitad de la longitud del 


elemento prismático» 


Durante la investigación del estado 
tensional de sistemas prismáticos de lon- 
gitud infinita o semiinfinita la solución ' 


más cómoda se obtiene en la forma: 


A O ANS E de 
1 2 3 4 
AO O E 
1 10% 29 $3 
a 
4 4 


(4.13) 


Los valores de Py se determinan por la 
sustitución de las expresiones (4.13) en 
las ecuaciones (4.1) para el caso di =D 
y por comparación de los coeficientes de 
e€ en potencias equivalentes. Para siste- 
mas infinitos o semiinfinitos a partir de 
las condiciones de amortiguamiento de los 
desplazamientos y de los esfuerzos puede 
comprobarse que en la medida en que nos 
alejamos del punto de aplicación de la 


carga es necesario colocar A =A =0, 
1 3 


La solución particular del sistema de 
ecuaciones (4.1) se determina de una 
forma semejante. Ella para q, F ( tiene 


la forma: 


Y, 


GS 
11 


(2) y) 
1 


U =03 (4.14) 


Teniendo las soluciones (4.9) y (4.14) 
de acuerdo con. (4.2) obtenemos la solu- 
ción general del sistema de ecuaciones 
(4.1). 


Con el aumento del espesor de las pla” 
cas en condiciones extrictamente iguales, 
crece el papel del coeficiente S¡,+ lo 


cual nos puede llevar al caso a? < b, 


Para este caso introduciremos la si- 


guiente nomenclatura; 


d-b=zivVb-a =1tc (4.15) 


y las ralces de la ecuación algebraica 


(4.6) se puede representar en la forma: 


r = tat1f (4.16) 
15213594 


Las magnitudes 4 y f se determinan en 
correspondencia con representaciones tri- 
gonométricas de varible compleja por las 


fórmulas: 


a =yp cos —— ¡86 =P sen — 
(4.17) 
donde : 


(4.18) 


A las raices (4.16) corresponde una 


solución del tipo: 


pan = A, cosh az cos Ba + 


+ A, cosh az sen faz + A, senh az cos Ba + 


+ A, senh az senh $z 


V lo) 


j = B, cosh €z cos Pz + 


+ B, cosh az sen Bz + B, senh az cosbz + 
+ B, senh 4z sen bz (4,19) 

Las constantes 5, se expresan a través 
de las constantes A de los dos sistemas 


independientes de ecuaciones algebraicas: 


ya 


0B, +88, = —Ll20B A, + (ap? - 
11 
7 
> Ya ) Ay]; 
11 
Ya 4 4 
83. 40h, AE TA e > 


- —LsaA -208 4] (4.20) 
2 2 3 


va 
ab, - BB, = E [- 2 0BA4, + (a? "-p? - 
á A 20, JAI 
es : 
va 
BB, +aB, ed e =p? - 
b 11 
Ds +20BA,] (4.21) 


El sistema (4.20) es mejor resolverlo 
en cada caso concreto sustituyendo los 
valores numéricos de los coeficientes en 


las ecuaciones de Az y Bo 


Por último, si ocurriese, el caso par- 
ticular extremo, en el cual A =b, en- 
tonces la ecuación (4.6) tendrá las dos 


- 
ralCces. 


P = Y a =b yr ==-=ya=-b 


(4.22 


En este caso la solución. general de 
las ecuaciones homoyéneas se puede escri- 
bir en la forma: 


( 


yen = A, cosh bz + A,2 senh bz + 


+ A, senh bz + A, 2 cosh bz 


0 = B,cosh ba + B,z senh ba + 


+ Bysenh ba + B,2 cosh bz (4.23) 


El vínculo de las constantes AL Y dz 
se establece como en el caso anterior por 
la sustitución de las expresiones (4.22) 
en la primera ecuación (4.1) y por la 
comparación de los coeficientes en el 


caso de funciones equivalentes, 


El caso a? <b ya = bno es caracte- 
rístico para sistemas de paredes delga- 
das, por lo que su análisis detallado se 
sale de los objetivos de la presente pu- 


blicación. 
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5. EJEMPLOS DE CÁLCULO 


Ejemplo 1: Calculemos el sistema mos- 


trado en la figura 5.1. 


"Teniendo en cuenta el carácter antisi- 
métrico de la deformación del sistema, 
seleccionemos las funciones e, (8) y Y, (6) 


en la forma mostrada en la figura 5.2. 


c) 


Figura 5.2 


La saeta en el gráfico e, aquí y en lo 
que sigue la dirigiremos a lo largo del 
contorno en el sentido del crecimiento 


algebraico de la función Pje 


Utilizando la regla de Vershaguin y 
las fórmulas (3.4) obtenemos: 


a. =0,25 sé =0,0254; 
1 


b =1,55d=0,15 d 
11 


e =3d01d4ir =bd=0,1 P 
11 


(5.1) 


Para el cálculo del coeficiente 5 
construiremos el gráfico 4, del desplaza- 
miento unitario del elemento horizontal 
del pórtico. El cálculo lo realizaremos 


por el método de los desplazamientos. 


La incógnita auxiliar la representare- 
mos por Z. Los gráficos 4, del despla- 
zamiento unitario señalado y M, de £, = 1 


se mestran en la figura 5.3. 


sobre la base de la figura 5.3 se obtiene 


la siguiente ecuación: 


qu. PAE LR 
22 d ie q 


por lo tanto, tenemos: 


Multiplicando el gráfico de 4, por 
este valor y sumándolo con el gráfico de 
A, obtenemos el gráfico de A, que necesi- 
tamos, el cual se muestra en la figu- 
ra 5.4, y que multiplicándolo por sí 


mismo nos permite obtener que: 


YB 3 
O aa 4 El 8 
pi A -J: El d3 = 2 ; LE == 12 
(5.2) 


Este mismo valor se puede encontrar 
más fácilmente si los 63... se determinan 
como la reacción en la ligadura introdu- 
cida [fig. 5.4] a través de la suma de 
las fuerzas cortantes en los elementos 


verticales del pórtico. 


Calculemos ahora el pórtico bajo la 


acción de la carga dada por el método de 
los desplazamientos en las condiciones en 
que los desplazamientos de los nudos no 


son lineales. 


Sin detenernos en la explicación del 
proceso que ya es conocido por los estu- 
diantes presentaremos el gráfico defini- 
tivo Ma (Edg..5:25)6 


Este mismo resultado se puede encontrar 
de forma más sencilla, si se determina (, 


como la reacción en la ligadura introdu- 


240 cida a través de la suma de las fuerzas 
cortantes como se muestra en la figu- 
ra 5.6. 
q, 
d 
d 4 q 
a p 2 > 240 240 
Figura 5.6 
Para el cálculo del trabajo de la 1 16 4 7 
J 4; o d + 240 ud + 240 qua = 3 e 


carga dada q establezcamos la ecuación 
diferencial del eje aislado de la columna 
izquierda del pórtico bajo el estado uni- 


tario 4, y realicemos su integración: 


a 2% = , 
E A A 
d d 
3 
y =- 7 PAN 


1l origen de coordenadas lo colocare- 
mos en el extremo inferior de la columna: 
y el eje Oy lo dirigiremos hacia la parte 


derecha . 


Entonces, las condiciones de borde 
para la determinación de las constantes 
C, y C¿ toman la forma y (0) = 0, y (4) = 1, 


a partir de las cuales C(, = 0, y 


4 
C, 34 luego tendremos: 
yy lo rs A 
3 a 3 


A partir de la fórmula anterior se 
puede calcular q, en la forma: 
d d 3 


e ia A o L 
qa / qydz =q: 1 ESE 33 


Ls 
3 qd (5.3) 
61 


Prestemos atención, si se utiliza la 
distribución de la carga a través de los 


nudos entonces tendremos: 


y el error relativo de este procedimiento 


en el caso dado será igual a: 

—- —_————— 10 e 100 = 14 5% 
q, 7 

el cual no es permisible, 


Para los cálculos que siguen, adoptare- 


mos que: 
y PEE 8 
G 


De acuerdo con las fórmulas (4.7) cal- 


culemos: 


0,625 d* 


donde no es difícil comprender que en el 
caso dado a? > hb, por eso de acuerdo con 


las fórmulas 4.8 encontraremos: 


a =0,16d*;f = 0,884 8 d”* 


después de lo cual y de acuerdo con las 


fórmulas (4.12) tenemos: 


IS E 


: Bj == 14142. Fa 


y ahora se puede escribir la solución ge- 


neral en la forma: 


U, = A, cosh 0,16 + +4, senh 0,16 y + 


1 
+4 5h 0,884 8 g +4 Z 
¿ COS A y , Ssenh 0,884 8 q 


Z 
Y, =- 9275 dA, senh 0,16 a 


- 9,275 dA,cosh 0,16 + - 
Z 
- 1,142 3 dA,semh 0,884 8 F|- 


- 1,142 3 dA,cosh 0,884 8 ++ 


+ 1,75. 10 H2 (5.4) 
Analicemos el caso en que el sistema 
esté apoyado en dos muros delgados (figu- 
ra 5.7). El origen de coordenadas lo co- 


locaremos a la mitad de la luz. 


Debido a la simetría de la deformación 
del sistema relativo al plano Z = 0 ten- 
dremos A, = A, = 0 y la solución (5.4) 
toma el tipo que se calcula por la ecua- 


ción siguiente: 


U, = A, senh 0,16 ++ 


Z 
A 


+ A, senh 0,884 8 


y] = 26. y cosh 0,16 ++ 
"62 


Figura 5.7 


V, == 9,275 dA, cosh 0,16 7 $ a 


- 1,142 3 dA,cosh 0,884 8 =7+ 


El 
+ 1,75 » 10 e 


— 


V/ =- 1,484 A, senh 0,16 3 


= 1,010 7 A,cosh 0,884 8 


Suponiendo l = 5 dl de acuerdo con las 


> 1 
tablas del manual de matemáticas , calcu- 


lemos para = = 255 
senh 0,4 = 0,410 8,senh 2,212 = 4,512 2 


cosh 0,4 = 1,081 1,cosh 2,212 = 4,621 8 


(5.6) 


Utilizando los valores (5.6) y las 


¿ É 
fórmulas (5.5) obtenemos para 2 


VU = 0,410 8 A, + 4,512 2 A, 
Á, A, 
ts == 
UN = 0,173 0 % + 4,089 4 q 


(5.7) 


V, $ - 10,027 dA, - 5,279 5 dA,+ 1 750 ¿> 
V/ = -0,609 6 A, - 4,560 54, 


Las condiciones de borde (3.14) 
U¡ (2,5 d) =0 y V, (2,5 d) = O nos llevan 


a las ecuaciones: 
07173 0 A, + 4,089 4 A, =0 
10,027 4, + 5,279 54, = 1,75. 10 yes 


la solución del cual es: 


A, = 178 505 $; A, = - 7,551 54-L 


2 


(5.8) 


Para estos valores de A, y A, a partir 


de las fórmulas (5.5) encontramos que: 


t = [ ( 
U! (oy = 21,88 L— (5.9) 
V, (0) = 103,00 H=— (5. 10) 


Comparemos este último resultado con 
los siguientes que se obtienen de dos 


procedimientos elementales: 


a) Del análisis de la placa superior de 
la cáscara prismática como vigas que 
se flexan en el plano de su mayor ri- 


gidez (luz más corta) encontramos que: 


5 q,1* 


TAL 


VI (0) = E 

Este resultado es 5,5 veces mayor que 
el valor real obtenido en (5.10). Si, 
además, se tomase en cuenta el efecto de 
la deformación por cortante de la viga, 
la diferencia entre estos resultados 


sería aún mayor. 


b) La flecha Y, (0) en el centro del sis- 
tema se puede intentar calcular como 


el desplazamiento de un pórtico plano 


independiente, para lo cual es necesa- 
rio despreciar la primera ecuación del 


sistema (4.1). 


En la segunda ecuación se deben despre- 
ciar todos los términos diferenciales, 
los cuales representan el aporte que 
brindan a cada uno de estos pórticos, 
los restantes pórticos paralelos que 
representan la continuidad de la es- 
tructura en la dirección Z a través de 


los esfuerzos tangenciales. 


De esta forma tendremos: 


el > 
A A 


pórtico o 1 


= 1750 GE 


obteniéndose aquí una magnitud del 
desplazamiento en el punto medio del sis- 
tema, que es diecisiete veces mayor que 


la real. 


Los resultados de este análisis evi- 
dencian el hecho que los procedimientos 
elementales presentados no son aplica- 
bles al cálculo de los desplazamientos en 
los sistemas espaciales de paredes del- 


gadas. 


Comparemos el gráfico de los momentos 
Flectores en la sección transversal de 
la franja del pórtico central de la es- 
tructura analizada como un sistema espa- 
cial de paredes delgadas y el gráfico 
correspondiente de su análisis como un 
pórtico aislado de acuerdo con las fóxr- 


mulas: 


1! 


Ú 


H, + 4, V, (0) 


cs MM Y pórtico 
pórtico-plano Pp 1.1 
Estos valores se muestran en la figu- 


ra 5.8, 


Ambos gráficos son cualitativamente 


diferentes y podemos observar que los ma- 


0,0005 qd* 


M 
sistemo especial 
de paredes delgada 


Figura 5.8 


yores momentos flectores se encuentran 
en secciones diferentes, y los valores se 
diferencian en casi tres veces uno al 


Otro. 


Utilizando el resultado (5.9) y la fi- 
gura 5.2 calculemos los valores máximos 
de la tensión normal. 


o(0, a = EU! (0) 9, (>) = 10,94 q 


(55117 


Mientras que del análisis de la placa 
superior del pórtico como una viga flexa- 


da en el plano horizontal obténemos que: 


ante 


E zu 109,4 q 


M 
U=—m.= 


Esto quiere decir diez veces mayor que 


la real. 


Destaquemos, además, que si se usa la 
teoría del pórtico plano esta tensión, en 


general, sería igual a cero. 


Pot lo tanto, queda evidenciado que 
los procedimientos elementales expuestos 
anteriormente no son ventajosos para el 
cálculo de las tensiones en los sistemas 


espaciales de paredes delgadas. 


Sobre la base de los resultados (5.8) 
y (5.7) tenemos que: 


U, (2,5 d) = 39,273 6 -L; 


vr 245 d) = -74,382 + 
los cuales permiten calcular las tensio- 
nes tangenciales en el lugar de unión del 
sistema espacial con los diafragmas de 


apoyo. 


- 
! 


= GLU, (2,5. d) ette) + Vf (25 dde, (4) = 


1 


15,709 q P/(5) - 29,753 q Y, (s)) 
(5.12) 
El gráfico de las tensiones tangencia- 


los calculado de acuerdo con esta fórmula 


se muestra en la figura 5.9. 


14,044 q 


7,854q 


Figura 5.9 


Para realizar la comparación calcule- 
mos la tensión tangencial en la placa su- 


perior usando un método elemental. 


y Slemental 


E 
media > 14983: 


La pequeña diferencia entre las ten- 


siones tangenciales reales en la placa 


superior y este último valor hallado se para au ésleuto del coericiente Y es 
11 


debe a la colaboración de los pórticos necesario construir el gráfico M, del 

elementales en la absorción de la carga. desplazamiento unitario del elemento ho- 
Ejemplo 2; Como siguiente ejemplo cal- rizontal en dicha dirección. 

culamos el sistema representado en la fi- La representación de los desplazamien- 

gura 5.10. Las funciones y, (8), p/(8) y tos se muestra en la figura 5.12. En co- 

Y, ($) se muestran en la figura 5.11. rrespondencia con estos desplazamientos 


el gráfico de 4, en el sistema base se 
Figura 5.10 muestra en la figura 5.13, donde se pre- 
senta también el gráfico buscado de M4, 
obtenido a partir de 4, por el equilibrio 


de los momentos en los nudos. 


a) 


o 1,5 


Figura 5.12 


Cc) 


Figura 5.11 


Utilizando la regla de Vershaguin y 
suponiendo 6 = 0,05. d, Y == = 2 de 
acuerdo con las fórmulas (3.4) encontra- 
remos que: 


a = q = IEA 
Y a 0,1 d ; PB OIE E 


ór  L a 9705 dé (5.13) 


Multiplicando el gráfico de M, por sí Figura 5.16 0 


mismo obtenemos: 


- 5? 
GS = 64,05 1-2 = 5,337 5 E ss 
11 A a o 
+ 
4 
(5.14) A 
S oa, 
Este resultado no es difícil de obte- E A e 5 
ner a partír de las condiciones de equi- pee 2 Py 3 = 
. A d d o 
librio 2%4_ = 0 para una parte del pórtico 0,038 8889 q 0,016 6667 q S 


(fig. 5.14), donde los esfuerzos en las 
secciones fueron determinados en corres- 


pondencia con el gráfico de M, . Bhora de acuerdo con la fórmula (4.7) 


calculemos que: 
24= 1,026 683 “y? ¡; b= 0,080 062 u”* 


donde en el caso dado 4? > bh, por lo 


tanto, utilizando la fórmula (4.5%) ton- 
dremos : 
a = 0,291 585 4 * ¡$ = 0,970 395 u*' 
Figura 5.14 y de acuerdo con las fórmulas 4.12: 
Calculando el pórtico por el método de pj =.>:9,997 01 Pp,  - 2,121 43 


los desplazamientos para la curga dada en 
cd a S Ñ En este caso, la solución general se 
condiciones de independencia lineal de > , 


: dea mwede presentar en la forma: 
sus nudos, se obtiene el gráfico de mo- E : 


mentos A, que se muestra en la figura ) = A,cosh 0,291 6 3 + 
5.15 a partir de las condiciones de equi- 

d z ys re - . r 
librio £M, O de una parte del pórtico + A, senh 0,291 6 Ls 
(£ig. 5.16) encontramos que: 

a 0,637 5 qd (5.15) + A, cosh 0,970 a 


al ba 


+ A, senh 0,970 4 


2 ' 
q,= 0,6375 gd Figura 5.15 


¿ Y, + 19799701 £, senh 0,291 6 => 


2 
0038 88941 


2 
0,0111 11q4 


- 9,997 01 A,cosh 0,291 6 +- 


- 2,121 13 A,senh 0,970 4 —- 
- 2,121 13 4,¿cosh 0,970 4 — - 


d 
- 477,152 E 


2 2 
0,105 555 qU 0,005 556 q* 
4 ; bl (5.16) 


Analicemos el caso en que en el extre- - Z 
: Y, =- A, (9,997 01 senh 0,291 6 —- - 
mo Z = 0 existe una ligadura totalmente 1 A 
empotrada que impide todo tipo de movi- - 2,121 13 senh 0,970 4 FS 
miento y un extremo totalmente libre en , 

el borde Z = l. - A, (9,997 OL cosh 0,291 6 > 


Sustituyendo primero a las condiciones z 
- 9,997 Ol cosh 0,970 4 q + 


para Z2= 0 
E h 0,970 4 $) 
y, (0) = 0 + 477,752 G - COS . e 
A, Z 
y Y, (0) +0 VI == =q—(2,915 13 cosh 0,291 6 q- - 
y desarrollando estas condiciones, obte- - 2,058 34 cosh 0,970 a 
nemos las ecuaciones: " 
e —— (2,915 13 senh 0,291 6 E - 
4 XA =-0 d c 
1 3 
Z 
-9,997 01 A, - 2,121 13 A, +447,752 a = 0 - 9,701 10 semh 0,970 4 —=) - 
A partir de las cuales despejando A, y -= 463,611 — senh 0,970 4 .. 


A, y sustituyendo en la solución (5.16) 
4 y por cuanto: 


encontramos que las fórmulas para los 


desplazamientos y las primeras derivadas (senh 0,291 6) ,4= 1,449 6 ; 


adoptan la forma: (senh 0,970 4) +4 = 24,240 ; 
(cosh 0,291 6) .4 = 1,761 0 ; 
2 
V, E A, (cosh 0,291 6 OS (cosh 0,970 4) + 4 = 24,262 
7 entonces: 
-= cosh 0,970 4 7 + 
, 0, (L) = - 22,501 A, —VIZ y TOD A, + 
+ A, (senk 0,291 6 rl 
+ 5 459,70 = 
- 4,713 06 -» senh 0,970 4 > + 
IV 23,099 8: . 110,449 
d ; y O 
+ 229 235 A senh 0,970 4 ná 
A + 5 302,90 == 
U] = ——(0,291 6 semh 0,291 6 —. - 
V, (1) = 36,925 54, + 224,942 A, - 


mn Z 
- 0,970 4 senh 0,970 4 + 
141 113,5 5 


A, 
da 


(0/291 6 cosh 0,291 6 — - 


d ” __ 44,805 9 230,930, 
a a Md 


- 4,573 55 cosh 0,970 4 eS + 


+ 


- 11 237,9 + 


2 d 
+ 21 — 
9,708 G cosh 0,970 4 d (5517) las condiciones en el borde libre Z = Ll. 
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UL) = 07 O. Y «to + 0 vi (2) = 0 En el sentido horizontal del contorno 


del extremo libre del sistema espacial 

llevan a 1 j igui : 

e as ecuaciones siguientes tenemos que 7 (L,8) = -0,09 q valor que es 
- 23,099 8.4, - 110,449 A, + E A MO 
proceso de cálculo. 

p 4d 
e REG 

> 302,90 G 0 En los lados inclinados del contorno 
libre: 


-0,009 8 A, + 0,267 0.4, - 15,925 se =0 


iba = GU, (1) p/ (8) = 7,22 y 
Resolviendo encontramos que: 


A, =-47,32 LE; A, = 57,907 4 E 


1 


Aquí las tensiones tangenciales tam- 
bién debieron ser iguales a cero, sin em- 


. e ran y Lc , y t d 
Sustituyendo estos valores en las fór- PO, CUPO 00 O 


mulas (5.17) obtenemos que: las hipótesis de partida que posee la 


teoría aproximada de cálculo presentada 


en este trabajo. 


U, (0) = 0; U[ (0) = -29,389 — ; 
Los gráficos de las tensiones se mues- 
V, (0) =0 7 V/(0) = 40,543 3 L, po e iD 


0, (1) = -7,22 a y UL) = 0,17 


d 
V, (1) = 165,03 , V!(L) = 14,35 + 


Calculando ahora las tensiones norma- 


les de acuerdo con la fórmula (2.6) 


ga 0 3 EU] (0) Pomáx. = EU/ (0) .1= 
= -58,78 q > 
ax. bie EUA 1) += EUJ(L) - 1= 


" 


0,34 
¿34 q 40,54 q 


Las diferencias de este último valor 
con respecto a cero se deben a impreci- 
siones en. el cálculo, a pesar de haberse 
realizado el mismo con un alto grado de 


precisión. 


Para las tensiones tangenciales de 


acuerdo con las fórmulas (2.7) obtenemos: 


1 (0,8) = GLU(O)p/ (s) + 


+ V/(0)Y, (s)] = 40,54 q 


7 (L,s) = G[U, (L)9] (8) + V, (1)Y, (8)] PIU 


, Los gráficos de los momentos flectores 
en la sección transversal de las franjas 
pórtico del borde libre y para un pórtico 
plano aislado que no forme parte de un 
sistema espacial de paredes delgadas lo 
construimos de acuerdo con las fórmulas: 


M=M 


b + 4, V, (2) 


É = pórtico; 
pórtico A, E 


y pórtico 
j ¡ 


Y se representa en la figura 5.18, 


2 
0,016 737 7qd 


A 


2 
opt 04 qd 


0,132 3729 qd* 


a) 0,021 261 9 qa* 


0,041 7317 dq 


b) 


0,183 1903 4d” 
Figura 5.18 


Estos gráficos se diferencian cualita- 
tivamente uno del otro y sus valores ex- 


tremos difieren considerablemente., 


Para otras secciones de la estructura 
espacial de paredes delgadas, más aleja- 
das de la sección libre, estas diferen- 


cias serán mayores aún. 
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Detengámonos ahora en un asunto funda- 
mental: las posibilidades de comprobar si 
los resultados del cálculo son correctas 


o no. 


Por cuanto, en la base de la teoría 
presentada descansa el método de los des- 
plazamientos, entonces, las comprobacio- 
nes de los resultados del cálculo están 
ligadas naturalmente con las condiciones 
de equilibrio del sistema en su conjunto 
y de cada una de sus partes de forma se- 


parada. 


Nosotros utilizamos las condiciones de 
equilibrio de la placa superior, separan- 
do la misma de las placas inclinadas que 
sirven de apoyo y del masivo que le da 
empotramiento en la sección Z= 0. Si 
aplican las fuerzas correspondientes 
al efecto de estos cuerpos, que dan pro- 
yección sobre el eje 0S, y se tiene en 
cuenta que pretendemos proyectar todas 
las fuerzas en la dirección de este eje 


(£ig. 5.19). 


y e L 
q,2= q(0) - Q(1) - G8,, SV, (Zidz = 0 

2 

Q(0)= 2,027 qd 
2 
q50,637 5 qd 

65, V,(Z) 
PT 


2 
Q(Y=0,00 4 9d 


Figura 5.19 


El miembro integral expresa el aporte 
de los pórticos elementales a la resis- 
tencia total del sistema frente a la 


carga actuante. 


Para evaluar el aporte de cada uno de 
los sumandos sustituyamos los valores co- 
rrespondientes de los mismos para el sis- 


tema analizado. 


2,55 q - 2,027 q - 0,004 qd - 


- 0,518 q4* =0 


. 
[o] 


2 
2,55 qd - 2,549 qd =0 


De esta forma la comprobación ha sido 


realizada con un alto grado de precisión. 


Es evidente de todo lo expresado ante- 
riormente que es posible realizar una 
considerable cantidad de comprobaciones 
de diferentes tipos sobre la base de com- 
probar las condiciones de equilibrio de 


una parte, o de todo el sistema. 


La comparación de los resultados obte- 
nidos, con los resultados que se obten- 
drían a través de los procedimientos ele- 
mentales que analizamos en el ejemplo an- 
terior los presentamos en forma resumida 


en la tabla 5.1. 


Tal y como vimos en los resultados del 
ejemplo anterior, se verifica la no co- 
rrespondencia de los métodos elementales 
para la determinación de las tensiones y 
los desplazamientos en los sistemas espa- 


ciales de paredes delgadas. 


TABLA 5.1 
Cálculo de 

Magnitud la cáscara 

4d 
Y, (2) 330,06 —¿—= 

d 
y, (2) -14,44 Y 
Máx. O) -58,78 q 
media) 40,54 q 


6, CÁLCULO DE SISTEMAS PRISMÁTICOS 
POR EL MÉTODO DE LAS SERIES 
TRIGONOMÉTRICAS 


Los ejemplos presentados en el epígra- 
fe anterior de cálculo de cáscaras sobre 
la base de la solución exacta del siste- 
ina de ecuaciones diferenciales de equili- 
brio, en función de los desplazamientos, 
muestra la gran laboriosidad de este pro- 
ceso y la necesidad de realizar el cálcu- 
lo conservando un elevado número de 


cifras significativas. 


Es totalmente evidente que el cálculo 
manual en el caso de sistemas prismáticos 
- * . 
con un número elevado de incógnitas es 
demasiado difícil, para no decir práctica- 
mente imposible, por lo que resulta avi- 
dente la necesidad de buscar métodos 


mucho más simples de solución de la tarea. 


Analicemos en calidad de ejemplo, el 
sistema representado en la figura 5.7 
para el cual ya nosotros obtuvimos la so- 
lución sobre la base de la solución 


exacta de las ecuaciones. 


Esto nos permitirá en lo que sigue 
comparar los resultados que obtendremos, 
con los alcanzados en la solución exacta 
y apartarnos del procedimiento descrito 
anteriormente de cálculo de los coefi- 


cientes y los términos independientes de 


Cálculo del Cálculo de 
voladizo pórtico plano 
4 896 . 955,50 o 
-816 = 0 
-=612 q 0 
51 q 0 
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las ecuaciones de equilibrio en función 


de los desplazamientos. 


. 


Las condiciones de borde en los extre- 
mos, esto es, en las secciones 2 = 0 y 


Z = l toman la forma: 
Y) =0 y U/ =0 (6.1) 


A estas (y valga la pena recalcar «ue 
sólo a estas) condiciones de borde es po- 
sible satisfacer, haciendo: 


me . o ” J MONA Z »h . 
U, (2) =U0,, cosA,Z2 + U,, cos Ajó +... 


Plays VU cen A 2 Vi Sen AZ + 


(6.2) 


== (6. 3) 


Sustituyendo las expresiones (6.2) en 


el sistema de ecuaciones: 


Yao BO E Viso 
11 mur E TAE 
1 
CC Uhap Ves Y. —— =0 
11 1 E] PS | E 1 3 (6.4) 
El sistema (6.4) en el caso general no 


se satisface. 


Se analiza la parte izquierda de las 
ecuaciones, diferente de cero, como una 
carga no equilibrada sobre la franja de 


pórtico elemental de ancho de = 1, 


Se obliga al trabajo de esta carga 
sobre los desplazamientos posibles 
cos A 2 y senA,Z a hacerse igual a cero, 
o lo que es lo mismo, las ecuaciones (6.4) 
volvamos a escribirlas ahora en la si- 


quiente forma: 
¡Uva =b JU cosA Z2+ (ya - 
0 111 LESA 1 2 11 


= ¡b JU. ¿oBA Ll Hs =CE6 AV e0sA2Z- 
11 12 2 11 11 1 


A 


V cosAZ- .)] cos A 2dZ= 0 
11 2 12 2 mM 


(6.5) 


o AU senAZ-=C AU send 2-..- 
0 tt, Y Li 1 IE A 6. 2 


y 2 
= (Y AO Al) sen A, L- (r, ¿A + 
y O 11 11 


+ 5 3)VM o send £ == es+ de sen A, Zdz = 0 
IS Tf 


pa - 

Lestaquemos que en el caso úauo el nú- 
mero de ecuaciones (6.5) es Jjgual al 

doble del número de sumandos que se con- 


serven en las series (6.2). 


Prestando atención a las propiedaces 
de ortogonelidad de las ¿unciones trigo- 


métricas, esto es a la igualdzd: 


E ; 
ENE Xx. 2d4 =- 0 
=) y cos A, 0; 


L : ; 
ff senmA, 2 send ZdZ =.0 (0.0) 
0 y Mi 
cuando 1 * im también tenemos Jes ¿¡Grisulas 
siguientes: 


l sae a de é 1 
J cos! A Zdl sj send A fu = — (uv. 7) 
0 Mi 0 m de 


Ubtenemos, cambiando todos los signos 


AS 1 
por su inverso y dividiendo por q que: 


B 2 f . 4 e 
ia i A el 11 1 mn 9 
(0.8) 
4 
2 Vii 
A + (1 MES Y... = : 
e mo Am Be A Wi ú 
donde se introdujo la nomenclatura: 
1 ms 1542 1% (2) sen A ZdZ (6.9) 
Mm pre 1 m 


En caso de estar actuando una fuerza 
uniformemente distribuida a lo largo de 
la longitud de vuelo de la estructura 


agus constante tenemos que: 


al 


AS . 10 
44m mí 4, (6. 10) 


donde: 


mz Número entero e impar. 


En el caso de estar actuando una 
fuerza: concentrada P, aplicada en la 
sección Z = 2, la integral en la fórmula 
(6.9) es necesario calcularla como el 
trabajo de esta fuerza sobre el desplaza- 
miento de su punto de aplicación, como 


resultado de todo esto tenemos que: 


Pone 1 sen Am Z 


í D (6.11) 


Uam 
A partir de las ecuaciones algebraicas 


(6.8) se determinan Y n y Va" 


Para el cálculo de las tensiones en 
cualquier punto del sistema, en correspon- 


dencia con (2.6), (2.7) y (6.2) tenemos? 


a 
I 


= HU! pp! =-EQMU send Z + 
1 1 O Y: 1 


+A,U, senA,Z+.)4, (8) 


— 
Y 


-G(U p” 
1 


+ VW) =Gl(U pp! + 
1 1.1 1 | 


+AV p)cosA Z+ 
e E 1 


+ (U 9 HA V $) 
12 4 2 11 


s ñ cos A Z + ...)]) 
1 2 


. (A cos A 2 +) 46.12) 


Comparemos ahora los valores obtenidos 
aquí para las flexiones horizontales de 


la placa superior, con el valor exacto 


V, = 103,00 E 


Conservando en la expresión (6.2) un 


FR, 


solo sumando tendremos que V, =104,84 


El error es de aproximadamente un 1,8 % . 


En la magnitud U, sobre los planos 
Z=0y Z= l en comparación con el valor 
exacto U, = 39,27 q/l tendremos un error 
del orden del 4 % tomando solamente el 


primer término de las series (6.2). 


Para las tensiones normales 0% en el 
centro de la luz, nosotros obtuvimos el 


valor de comparación (etalon) 0 = 10,94 q. 
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Conservando solo el primer término en las 
series (6.2) tendremos que 0 = 11,85 q 


(donde el error supera el 8 %). 


Conservando el primero y el tercer tér- 


mino 0 = 10,72 (resultado en un 2 % menor 


que el exacto). El proceso de convergen- 
cia para las tensiones tangenciales es 

aún más lento y por ejemplo, para obtener 
valores de T con un 5 % de error es nece- 


sario conservar en la serie el primero, 


tercero, quinto, séptimo y el noveno tér- 


mino . 

Destaquemos que bajo otras cargas, por 
ejemplo, una fuerza concentrada us. la sec- 
ción 2 = 1/2 la influencia de la conver- 
gencia del proceso puede ser aún más 
lenta, lo cual exigiría conservar en 1e> 
serios (6.2) un «yor número de tórninos. 

li general, es bueno senalar (que este 
mátodo de reducir el sistima de ecuacio- 
nes diferenciales a un gistema de ecua- 


ciones algelraicas con ayuda de las se- 
ries trigonométricas resulta sumamente 


efectivo. 


En la conclusión de este párrafo des- 
taquemos que la aplicación de las funcio- 
nes fundamentales de las vibraciones pro- 
pias de las vigas o de las funciones de 
la curva estadística tropieza con dificul- 


tades: 


En primex lugar, durante su utiliza- 
ción no se obtiene un sistema independien- 
te de ecuaciones del tipo (6.8) para cada 
miembro de la serie y, en segundo lugar, 
la rapidez de la convergencia de las se- 
ries, en particular, para el caso de un 
extremo empotrado y otro libre deja mucho 


que desear. 


Por todo lo cual para el caso general 
de fijación de los extremos del sistema 
prismático es conveniente la aplicación 


de técnicas de computación y diagramas 


IBM, por ejemplo, en la forma en que esto 


se presenta en el siguiente epígrafe 7. 


7, APLICACIÓN DE LA MÁQUINA 
COMPUTADORA NATRI AL CÁLCULO 
DE SISTEMAS PRISMÁTICOS 


En calidad de ejemplo analicemos el 
sistema prismático, calculado en el epí- 
grafe 5 (fig. 5.10). Adoptemos directa- 


mente los siguientes datos de partidas: 


Ya +- 2 0 e 05800 Mi 150; 


q =ú 
En el caso de estos valores y en co- 
rrespondencia con las fórmulas (5.13-5.14) 


tenemos: 


2 

va a a OS 
11 11 11 

A * » » 

18 A A =10/001 334: 375: 
11 11 

4, 

mado” 3,187, 54 


y las ecuaciones fundamentales adoptan 


la forma: 

245 UN - 0434, 

045: 07 + 1025 vi = 0,001 334 375 Y, + 
+ 341875: =-0 


be Y) 


y las correspondientes a ellas, condicio- 


nes de borde serán: 


U (0) =03;V (0) = 0 
1 1 


U 


Uf (d) 0.7. 045 U, (1) + 1,25 vr CL) ='0 


GL-2) 
Para la integración numérica de las 
ecuaciones diferenciales (7.1) en la má- 


quina computadora NAIRÍ haciendo el cam- 
bio de variables. 


== y? = = U = 
ECxilib=l, + y=Y19=V (7.3) 
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representar las ecuaciones dadas en fonuia 


de un sistema normal: 


El = 0,12: b: + 052 y 
bi = 07 
y! =-=0,4 C + 0,001 067 5 y = 2,55 


(7.4) 


Lie 


Por cuanto la mitad de las condiciones 
de borde están dedas para / = Ll, nosotros 
delemos llevar la tarea de trontera al 
caso de la tarea de Cauchy con condicio- 


nes iniciales dadas. Para ello, resolve- 


remos dos veces el sistema homoyéneo 

(despreciando -2,55) en la tercera ccua- 
LA : SN m7 

ción, bajo las siguientes condiciones 


iniciales: 


CO) = Lo DO) "=00 3 410) = 0 73 y(0) = u 
(tarea 1) (7.5) 
y. ¿en el ¡casor 

C(0) =U ¿ 0(0) =0 7 y(0) = Ll ; g(0) = O 
(tarea 2) (7.6) 


además, resolvemos el sistema 7.4 ¡para 


condiciones iniciales nulas 


C(0) =0 7 D0(0) =0 3; y(0) =U ; g(0) = 0 


(tarea (Y) 7.7) 


Si se obtienen los valores de las un- 
ciones buscadas en el centro del interva- 
lo de integración para todas las tareas 
parciales, nosotros podemos satisfacer 


las condiciones iniciales para Z= Ll. 


Como resultado de lo cual, encontrare- 
mos el valor inicial restante C (0) = 
= U! (0) y y(0) = V*(0) con lo cual la 
tarea, resulta reducida a la tarea inicial 


de Cauchy. 


Destaquemos que en el proceso de ob- 
tención del sistema normal del tipo (7.4) 


frecuentemente resulta conveniente que se 


hagan ortogonales las funciones y, (s). y 
Y, (8) como resultado de lo cual resultan 
iguales a cero todos los coeficientes su-' 


plementarios ee y Pip" 


Antes de pasar a la realización de 
esta solución de acuerdo con los proyra- 
mas estándares para nuestro caso concre- 
to, recordemos que en calidad de símbolos 
para las funciones buscadas se pueden 
utilizar los siguientes: 2, b y Y,N, M, 
EAN CUY e ¡OE Pi Mia Ds 
Ss bir mientras que para la variable in- 


dependiente se usa la -(t). 


De acuerdo con FeoktistowW la máquina 
computadora NAIRI puede resolver los sis- 
temas de ecuaciones del tipo (7.4) hasta 
de orden 22 inclusive, utilizando un paso 
constante de integración y hasta de orden 
17 para el caso de una precisión dada como 


esto se describe en lo que sigue: 


La llamada a los programas y la entra- 
da de los sistemas (7.4) en la máquina 


se realiza de la siguiente forma (para la 


tarea 1). 

66 A 

U,UUU 1 (Lk) 
78 n 

0,5 (bk) 
20 (bk) 
U (bk) 
1 (Lk) 
0 (bx) 
0 (bk) 
U (bk) 
W 

dy 


El segundo número de arriba (0,000 1) 
determina el grado de precisión estable- 
cido por nosotros. En la quinta fila 
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establecemos el paso de la integración 
(0,5). 


dican el final y el comienzo del interva- 


Los dos números que le siguen in- 
lo de integración. Los restantes nfímeros 
indican las condiciones iniciales corres- 
El símbolo (bk) 


indica que es necesario dar retroceso al 


pondientes a la tarea 1, 


carrete. 


Después de la selección de dy la má- 
«quina misma da el retroceso al carrete y 
pide informe sobre el orden del sistema 


introducido, escribiendo n= . 


Es necesario escribir el orden del 


sistema (en nuestro caso 4) y dar el re- 


troceso al carrete, después de lo cual la 


máquina escribe ('' A 


Es necesario indicar solre el tecla, 
la parte derecha de lu primera ecuación 
y dar retroceso al carrete, luego, tan 


4 


pronto la máquina escrile h' = es ne- 


cesario introducir la parte derecha de 
la segunda ecuación, y así continuar 
hasta tanto no sean introducidas todas 


las ¡unciones del sistema. 


En respuesta a lo anterior, la máqui- 
na pasa a realizar la integración escri- 
biendo los resultados después de cada 


paso. 


La continuación de la escritura inte- 


rrumpida tiene la forma: 


n=4 (Lk) 
e* =-0,12 b + 0,2 y (Lk) 
DiE=e (Lk) 
y' = 0,4 + 0,001 067 5 (kk) 
g*=y (LK) 


Nótese que en la tercera ecuación fue 


eliminado el término independiente (-2,55). 


En lo que sigue la máquina escribe 


sola y al final del proceso ella escri- 


birá: 
t = 20,000 000 00 


Co = 


25,792 805 43 


b = 131,710 519 79 


y = -53,930 904 86 


gy = -262,650 505 06 
y se detiene. 


La solución de la tarea 2 se realiza 
de la misma forma estando la única dife- 
rencia en las condiciones iniciales (7.6), 


lo cual lleva a los siguientes resultados. 


20,000 000 00 


hb = 131, 325 263 97 
y = -52,301 402 56 
yg = -215,533 357 1 


Por fin, la solución de la tarea Y 
tendrá solo dos diferencias con las ta- 


reas anteriores. 


Petmera diferencia: todas las condi- 
ciones iniciales son nulas. 

Segunda «dtjfereneta: el sistema (7.4) 
debe ser introducido de manera completa, 
esto quiere decir con el .sumando (-2,55) 
en la tercera ecuación, 

Al final del proceso la máquina in- 


prime: 


t 


20,000 000 06 


(él 


n 


-334,879 413 60 


= 
u 


-1 489,032 760 62 
5 499,610 084 53 


1 842,444 046 02 


a] 
n 


Ahora es posible satisfacer las condi- 
ciones de borde en el extremo libre ! 


Z=1l= 20. Para esto es necesario. en 


15 


odrrespondencia con las dos últimas con- 
diciones (7.2), componer y resolver las 


ecuaciones: 


C (2)U'*(0) + C (L)V*(0) + C (1) =0 
1 1 2 1 4 
[Cc bl) «rr y (23) U*(O) + [COC b (Ll) + 
11 3 11% 1 11 2 
+ 2, Y, (DIV, (o) + Ey + 


+ Pg = 0 (7.8) 


El índice delante de las funciones 
e(1), b(1) y y(L) indica la nomenclatura 
condicional de las tareas parciales re- 


sueltas anteriormente. 


Escribamos las ecuaciones (7,8) «un 
forma nunérica: 
25,792 805 43 U* (0) + 26,965 455 29 vr to)- 
- 334,879 413 60 = 0 


(0,5 + 131,710 519 79 - 
- 1,25 + 53,930 904 86) Ur to) + 
+ (0,5 » 131,325 263 97 - 


- 1,25 * 52,301 401 56) Y; (0) ', 
- 0,5 + 1 489,032 760 62 + 


+ 1,25 *» 549,610 084 53= 0 


Realizando las operaciones aritméticas 


obtenemos 
25,792 805 43 Ut (o) + 
+ 26,965 455 29 ye (0) = 


= 334,879 413 60 -1,558 371 14d (0) + 


+ 0,285 876 Y (0) = 57 503 80 
que son las ecuaciones obtenidas en forma 


numérica, 


Aunque estas últimas ecuaciones no son 
difíciles de resolver todo, se puede re- 
solver en la NAIRI por sus programas es- 


tándares. 


Para ello los coeficientes y los tér- 
minos independientes son necesarios in- 


troducirlos por filas, por ejemplo: 
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25,792 805 43 (pk) 


26,965 455 29 (1k) 
334,879 413 60 (bk) 
-1,558 371 00 (bk) 
0,285 876 00 (bk) 
57,503 800 00 (bk) 
W 

Cy 

Después, la máquina pide el orden del 


sistema introducido escribiendo -. = des- 
pués de lo cual es necesario teciear el 
orden del sistema (en nuestro caso, 2) y 


dar retroceso al carrete. 


El programa permite resolver sistemas 
de ecuaciones lineales hasta de orden 28 


(incluido). 


Si a la máquina se le dan sistemas de 
orden superior a 28, automáticamente es- 


cribe "no puedo" y la máquina se detiene. 


Si se le da un sistema no qompatible, 
entonces escribe "no tiene solución". Si 
la solución no es Única, entonces, impri- 


me "no es capaz de resolverlo", 


El resultado de la solución del siste- 


ma de ecuaciones algebraicas lineales es: 
4, =U!(0) = -29,453 620 ¡ 


= yA = 
o v, (0) = 40,591 569 


Ahora, por última vez, resolvemos el 
sistema de ecuaciones (7.6) en la máquina, 


bajo las siguientes condiciones iniciales, 


-29,453 620; b= 0 ; y = 40,591 569 ; 


' 


e 
(7.9) 


No es difícil comprender que para la 
realización de esto es necesario marcar 


sobre el teclado del teletipo. 


66 h (Ek) 
0,000 1 
u 
78 n (Lk) 
0,5 (Lk) 
20 (bk) 
Ú (Lk) 
-29,453 620  (Lk) 
0 (Lk) 
40,591 569 (Lk) 
0 (pk) 
Ww 
dl 
n=4 (bk) 
ce" ="0,12b+0,2y  (bk) 
b'=e 
y' =-0,4 € + 0,001 067 59 -2,55 — (Ek) 
g". y (bk) 


Después de lo cual la propia máquina 
escribe los valores numéricos a través 
del intervalo fijado anteriormente 
(h = 0,5) o divide aún más el mismo, si 
con él no se puede alcanzar la precisión 


exigida de (€ = 0,000 1). 


Teniendo en cuenta los valores numéri- 


cos de las funciones csU, b=U, , y=V/ y 


g= an no es difícil construir los yráfi- 
cos de la flecha Ls de las tensiones nor- 
males a lo largo de los bordes 6 = 


SUE. nes pi + y de las tensiones tan- 


genciales T = Guo + Y Y). 


La comprobación de los resultados del 
cálculo se realiza tal y como se hizo al 


final del epígrafe 5. 


La Única diferencia consiste en que la 


integral en la ecuación de comprobación. 


(5.18) no puede ser tomada en forma cerra-: 


da, así como que la máquina da solo valo- 

res numéricos de V en diferentes puntos 
1 

del intervalo de integración y no su 


función. 


Pero en este caso es necesario cambiar 
la integración por una suma, la cual se 
puede realizar, utilizando la regla del 
trapecio, la de Simpson, u otras. 

El algoritmo descrito es aplicable 
también al cálculo de sistemas espaciales 
de paredes delgauas con ayuda de otras 
máquinas computadoras, en las cuales el 
programa de solución de la tarea inicial 
(u original) para el sistema de ecuacioó- 
nes diferenciales ordinarias no aparece 


coripilado en la propia máquina. 


En el presente, para todas las máqui- 
nas IBM se tiene una biblioteca suficien- 
temente amplia de programas estándares, 
entre los cuales es posible encontrar el : 
programa de solución de la tarea de Cau- 
chy para ecuaciones diferenciales ordina- 


rias. 


S, PEQUEÑA INTRODUCCIÓN HISTÓRICA 
Y DEL ESTADO ACTUAL DE LA TEORÍA 
DE LOS SISTEMAS ESPACIALES DE PAREDES 
DELGADAS 


El origen del desarrollo de la teoría 


presentada en la literatura nacional se 


encuentra en los trabajos!» * de V.Z. Vla- 


sov y de P.L. Pasternak. *> * 


En el trabajo de Vlasov! , las ceforma- 
ciones y los esfuerzos en el elemento 
del sistema plegado (placas rectangula- 
res aisladas) así como la carga exterior 
se representan en forma de series de 
Fourier. 


Por primera vez en comparación con los 
autores de los trabajos de las referen- 
cias 5, 6, 7 y B, se presenta un análi- 
sis de la aplicación a los sistemas ple- 
gados de los métodos de las fuerzas, de 
los desplazamientos y del método mixto, 
mostrándose que al método mixto corres- 
ponden ecuaciones algebraicas de ocho 
miembros de estructura simple, 

sas ecuaciones fueron obtenidas des- 
pués del análisis de los estados elemen- 


tales del sistema base del mótodo mixto. 


im diforencia con V.Z. Vlasov y eviden- 
temente de forma independiente a él, 
qm NN E 
P.L, Fastemak , así como E, Grúber , 
partieron en su análisis del sistema base 


del método de las fuerzas. 


En calidad de incógnitas fueron adop- 
tados los momentos flectores en la direc- 
ción transversal y las fuerzas deslizan- 


tes a lo largo de los bordes del sistema. 


A diferencia del trabajo primero de 
E. Grúber donde sobre la base del método 
de la fuerza fueron establecidas primera- 
mente las ecuaciones diferenciales ordi- 
narias, y solo luego con ayuda de las se- 
ries trigonométricas fue realizado el 
paso a las ecuaciones algebraicas, Pibe 
Pasternak, desde el mismo comienzo, repre- 
sentaba todas las funciones a través de 
series de Fourier y obtiene para su solu- 
ción un sistema de ecuaciones algebraicas 
de doce miembros mayor que el de V. Z. 
Vlasov para la solución del cual recomien- 


da utilizar el método de Gauss. __. 


A pesar, de ello, luego el propio 
P.L. Pasternak explicaba a sus alumnos 
las ecuaciones de ocho miembros de 
V.Z. Vlasov (por ejemplo, en su cu 
especial de estructuras de hormigón ar- 


mado). 
De los trabajos de autores extranjeros 


correspondientes 4 este período, resulta 
justo destacar también los trabajos co- 


rrespondientes a las referencias 9 y 10. 


En el libro de V.Z. Vlasov ,? se da una 
expresión mucho más completa del método, 
que se complementa con ejemplos y se com- 
paran los resultados alcanzados con otros 
más precisos obtenidos a partir de la 


teoría de la flexión. 


Aquí por primera vez se presentan la" 
ecuaciones diferenciales de ocho micmbros 


(no algebraicas) .* 


Eso es conveniente para el cálculo de 
sistemas plegados bajo condiciones arbi- 
trarias de apoyo en las seccione: uxLre- 
mas y no solo en el caso de apoyo de los 
extremos sobre diafragmas delgados de 
condiciones ideales. 


% Gio un de- 


En el trabajo V.Z2. Vlasov 
sarrollo consecuente a sus ideas. Él 
desarrolló el proceso de la integración 
analítica exacta del sistema de ecuacio- 
nes diferenciales ordinarias de ocho 
miembros y mostró su gran laboriosidad e 


inconveniencia. 


Para la solución aproximada de las 
ecuaciones diferenciales V.Z. Vlasov pro- 
puso aplicar las funciones fundamentales 
de las vibraciones transversales de las 
vigás; las cuales transforman el sistana, 


en un sistema de ecuaciones algebraicas. 


Se analizó, además, el estallecimiento 
de las ecuaciones de ocho miembros bajo 
diferentes condiciones en los bordes lon- 


gitudinales rectos. 


Aquí se analizan, además, otras posi- 
bilidades de simplificación (por ejemplo, 
las que se producen de considerar inde- 
lormalble el contorno de la sección trans- 
versal), así como posibilidades de per- 
feccionamiento y elevación de la preci- 
sión en otros casos ligados a tener en 
cuenta el efecto de los momentos flecto- 


ves en la dirección longitudinal, 


dl a 
Al lector le resultara interesante co- 
nocer que por la escritura uel libro le 
fue otorgado a V.Z. Vlasov el yredo cien- 


tíftaico de doctor en ciencias técnicas. 


A diferencia de las puld caciones des- 


critas hasta aquí, destinadas al cálculo 
de culier tás plézgadas donde no se Cuna 
en cuenta el elcuto de la defurdtue ¿ón a 
vizallanmiento en el plano medio, es el 


trabajo de Viasov,'” ¡or 1rimera vez se 


tuvo en cuenta ul electr que provoca 


esta deformación. 

Aquí Ve. vlasov partió del método de 
los desplezamientos, y por primera vez 
desarrolló un nuevo método que consiste 
en separar del sistema dos secciones 
transversales extremadamente cercana una 
a la otra, y de esta forma se extrae una 


franja de pórtico, 

Sobre ella, con ayuda del principio de 
los desplazamientos posilles de Lagrange, 
se establecen las ecuaciones diferencia- 
les ordinarias de equilibrio en función 
de los desplazamientos (ver el epígrafe 3 


de esta publicación). 


Se resuelven, además, ejemplos concre- 
tos, y se presentan las ecuaciones diná- 
micas teniendo en cuenta el efecto de las 


fuerzas tangenciales de inercia. 


Se establecen ecuaciones equivalentes 
de equilibrio en función de los desplaza- 
mientos para la tarea plana de la teoría 


de la elasticidad. 


En el trabajo capital de V.Z. Vlasov,'” 


y sus modificaciones sustanciales en las 
diferentes ediciones, se presenta el 
cálculo de los sistemas plegados sobre la 
Lase del sistema de ecuaciones de ocho 
miembros y del método variacional de los 
desplazamientos para sistemas espaciales 
de paredes delgadas con sección transver- 
sal cerrada para la tarea plana de la 
Teoría de la elasticidad, de la flexión 
de las placas y de los sistemas de placas 


de Lhorde no mixto. 


Como consecuencia de la aplicación del 
método variacional de los desplazamientos 
a las tareas plana y espacial de la Teo- 
ría de la elasticidad se obtuvieron mode- 
los de cálculo de bases elásticas de una 


y tres capas. 


Estas son ampliamente utilizadas en la 
elaboración de trabajos de disertación 


para la oltención de grados científicos. 


Fn el libro'* además, se analizan los 
sistemas plegados estáticamente determi- 
nados, 


vibraciones con ocho miembros teniendo en 


se obtiene las ecuaciones de las 


cuenta las masas del sistema sobre las 


aristas o bordes longitudinales. 


Se muestran las ecuaciones diferen- 
ciales de equilibrio en función de los 
desplazamientos para un sistema prismáti- 
co piramidal (o lo que es lo mismo) sua- 
vemente cónico con el contorno no defor- 


mable de la sección transversal. 


Se presenta una tarea muy interesante 
sobre el cálculo del sistema, teniendo en 
cuenta su proceso constructivo y se dan 


los resultados del mismo. 


Se desarrolla la teoría de la estabi- 
lidad de los sistemas espaciales de pare- 
des delgadas, basada en la hipótesis 


sobre la distribución plana de las ten- 
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siones normales en las secciones transver- 


sales hasta la pérdida de la estabilidad. 


Se tratan temas referentes a la úeter- 
minación de las tensiones provocadas por 
los cambios de temperatura y el cálculo 
de sistemas pretensados. Se dan solución 
a una serie de nuevas e importantes tareas 


de significación práctica. 


En el trabajo de 1.E. Mileikoskovo'* 
el método de V.Z. Vlasov se extendió a 4 
las tareas tridimensionales de la Teoría 
de la elasticidad y se obtuvo tanto un 
sistema de ecuaciones diferenciales en 
derivadas parciales de dos variables, 
como un sistema de ecuaciones diferencia- 


les ordinarias. 


Aquí se exponen los aspectos de la 
elaboración de la hipótesis de las sec- 
ciones planas, y también el desarrollo 
del sistema de ecuaciones en ecuaciones 


independientes. 


En la base del libro de I.E. Mileikov- 
kovo > se establecen las hipótesis sobre 
ausencia de alargamiento del contorno y 
de la deformación por cizallamiento en la 
superficie media, pero teniendo en cuenta 
todas las fuerzas interiores, sobre la 
Lase de la teoría de la flexión y además 
una pequeña curvatura en la dirección lon- 


gitudinal. 


Se propone una selección original de las 
funciones de distribución de los desplaza- 
mientos vinculada con los perfiles de pare- 
des delgadas, y se obtiene un sistema de 
ecuaciones diferenciales ordinarias en fun- 
ción de los desplazamientos, teniendo en 
cuenta las posibles simplificaciones del 
mi smo » Ñ 

Aquí también se muestra la elaloración 
de la teoría de las secciones planas y se 


expone el desarrollo del sistema de ecua- 


ciones originales (de partida) en ecuacio- 
nes independientes de cuaxrto:grado a tra- 
vés de la solución de la ecuación de 
"estado". 


Además, de la utilización de las funcio- 
nes fundamentales de la vibración de la 
viga se propone utilizar, para el cálculo 
aproximado, funciones de la flecha estáti- 
ca de las vigas de la carga correspon- 


diente. 


En las ecuaciones de las cáscaras cortas 
se recomienda despreciar los momentos flea- 
tores transversales con la excspción de 
la zona en que se une la cáscara con el 


elemento de borde. 


Se analiza el cálculo de las fisuras y 
de las deformaciones plásticas del hormi- 
gón para cáscaras de mediana longitud 
sobre la base de la teoría de formación 


de fisuras de V.I. Murasheva. 


A este trabajo se dedica también el 
libro de B.S. Vasilikov,'* El corre spon- 
diente desarrollo de las ideas anteriores 
en el cálculo de las cáscaras de hormigón 
amado, teniendo en cuenta la formación 
de fisuras se presenta en el libro de 
I.F. Mileikoskovo.'* 


La comprobación experimental de los re- 
sultados teóricos, realizada en el Institu- 
to Central de Investigaciones Científicas 
de Obras Estructurales, de nombre Kuche- 
renko, demostró lo adecuado de dichos re- 
sultados., 


En los libros indicados anteriormente 
se presenta un elevado número de aplicacio- 


nes a tareas concretas de la construcción, 


En ellos se presentan ejemplos de cálcu- 
lo de cubiertas de edificios, de presas de 
arco aligeradas y otras Obras; ejemplos de 
cálculo bajo cargas dinámicas y estáticas, 


y también bajo cambios de temperatura de 
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compuertas hidrotécnicas planas se presen- 


tan en llos trabajos de Pastuchijin.'* 


Un gran papel en el desarrollo de la teo- 
ría de los sistemas espaciales de paredes 
delgadas y de aplicación de la misma al 
cálculo de elementos de aviones de paredes 
delgadas lo constituyeron los trabajos de 
A.A. Umanskiv,17s 1 


Pero su autor no se limitó solo a la 
elaboración de publicaciones científicas, 
sino que, además, escribió el primer libro 
de texto,'”'sobre la teoría de los elenen- 
tos de paredes delgadas de los aviones , 
anplianente utilizado, además , por espe- 
Ccialistas de nunerosas ramas de la téc- 
nica. . 

Esto ha desempeñado un papel conside- 
rable en la preparación de ingenieros y 


científicos. 
lo desempeña, el libro de A.A. Umanski ,'* 


En el presente, este papel 


el cual recoge, tanto los resultados de 
las investigaciones científicas del autor, 
como una serie de logros alcanzados pos- 


teriormente en la teoría de las cáscaras, 


Los trabajos, realizados por los colakbo- 
radores del Instituto de Aviación de Moscú 
I.F. Obrazsovi, A.N. Ecpatievski, B.A. Ko- 
novalovim, G.G. Onanovim y otros, ligados 
con el desarrollo de la teoría de V.Z, Vla- 
sov y de su aplicación a los aspectos de la 
mecánica de la construcción de los aviones, 
aparecen en las publicaciones de I.F. 


Obrasov.*?+?0,21,22 


En el trabajo Obrasov, evidentemente 
se investiga por primera vez el llamado 
"tipo de flecha del ala", medio empotrado 
en el plano centroidal absolutamente 
rígido. 

Una representación más completa de la 
aplicación de esta teoría a la mecánica 
de la construcción de los aviones se pre- 


senta en la referencia 22 donde se expo- 


nen ño solo los trabajos anteriores del 
autor sino también de otros científicos, 


in parte $e presentan las investigacio= 
nes de cáscaras de pequeña conicidad 
(A.N, Elpetievekid, BA. Konovalov) la ob= 
tención de la expresión de la energía po- 
tencial total del sistema de ecuaciones di= 
ferenciales ordinarias, teniendo en cuenta 
el efecto de la temperatura (A.N, Elpa- 
tievskid), se tratan además, los aspeéctos 
de las tensiones de temperatura y las vi- 
raciones propias del cajón de pección 


transversal reutamgular (Astl. Blpatievskid)a 


He propone la solución en forma mixta, cuan- 


do en las ecuaciones de equilibrio, además 
de aparecer los desplazamientos, aparecen 


los momentos. 


Bn este caño, en cadá ecuación aparecen 
solamente las funciones buscadas y sus de- 
tivadas de primer ordena sto Crea las 
condicjones para la integración numérica 
inmediata del sistema de ecuaciones dife- 
rencialeb ordinarias en la máquina computa= 


dora 1BM (NiAs Conoloseba)a 


Aquí se presentan, además, las ecuacio= 
nes de equilibrio en función de los despla- 
Z2amientos en eoordenadas no reectangularos 
(G+G. Ononov) para cáscaras prismáticas y 
piramidales, 


Además, se analiza un ejemplo conereto y 
se presenta su comparación con el experi= 
mento, Be expone, además, el ediculo de 


alas en forma de plagas en voladizo. 


En el artículo de Pastuchijiny? ee 
propone un sistema de ecuaciones diferen= 
ciales ordinarias, obtenidas a partir de 
la expresión de la energía poteneial total. 


Esta tiene en cuenta, tante todos los 
factores de las fuerzas interiores, earac= 
terísticos en la teoría técnica de las cás=- 
caras como la defomación por eisgallamien= 
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to en la superficie media, la cual puede 
exigirse, por ejemplo, en el cálculo de 
puentes relativamente anchos en las ciu- 


dades. 


Á partir de este pequeno resullen, no es 
difícil comprender que la teoría de cálculo 
de los sistemas espáciales de paredes del- 
yadas, aleanzó un considerable desarrollo 
en huestro país y ha tenido innumerables y 


diversas aplicaciones. 


Una visión más completa de esto puede 
ser encontrada por el lector en los artícu- 
los resfímenes de O.La Ondachvili y Als Pa= 
novko 4,2% donde be presenta una amplia Li- 
Lliograffa sobre este tema, aunque aún es 


14d tadas 


kn los Úáltimos años hemos tenido conoci-= 
tdento de los trabajos realizados en la 
aplicación del método de los elementos fi- 
nitos al cálculo de los sistemas espaciales 


de paredes delgadasa 


Con la aplicación de potentes máquinas 
computadoras, este método puede resultar 


considerablemente efectivos. 
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NOTA DEL TRADUCTOR 

El. trabajo que hoy presentamos a consi- 
deración del especialista cubano en rela- 
ción con el cálculo de los sistemas espa- 
ciales de paredes delgadas presenta las hi- 
pótesis fundamentales y algunos de los 
métodos fundamentales de cálculo de este 


tipo de estructuras. 


Se tienen en cuenta el trabajo espacial 
de las mismas, así como las diferencias 
cualitativas de los estados de tensión de- 
formación que en los mismos se producen de 
aquellos que se obtienen de algunos esque- 
mas simplificados (y en la mayoría de los 
casos, esquemas planos) aún de frecuente 


uso en nuestro país. 


Es de destacar, además, las multifacéti- 
cas posibilidades de utilización de esta 
teoría en diferentes ramas de la técnica, 
aunque la presente publicación está diri- 
:gida fundamentalmente al análisis de es- 
tructuras prisiticas, en particular de ; 
aquellas que conforman las naves industria- 


les con cubiertas plegadas. 


En este sentido, es necesario destacar 
la estrecha vinculación que debe existir 
entre la solución constructiva adoptada y 
las características del esquema de análisis 


a utilizar. 


- Si bien resulta incorrecto utilizar es- 
quemas de análisis planos para calcular 
estructuras espaciales la adopción de es- 
quemas de análisis espaciales (como los que 
se presentan en este trabajo) para el aná- 
lisis de estructuras de cubierta y cerra- 
mento que no sean capaces de garantizar 
estas condiciones de trabajo puede implicar 


el colapso de la misma. 


Sin enbargo, conociendo esta teoría, sus 
hipótesis de cálculo y los métodos que uti- 
liza (los cuales están muy cercanos a los 


métodos tradicionales que conocen nuestros 


ingenieros, método de las fuerzas, método 
de los desplazamientos y método mixto) es 
mucho más fácil orientarse sobre la conve- 


niencia de garantizar constructivamente 


los vínculos que permitan el trabajo espa- 


cial de la misma, 


Eso sucede en algunos casos, a través 
de las uniones de los elementos prefabri- 
cados u otras soluciones constructivas o 
la conveniencia de aislar los elementos 
provocando que el comportamiento de la es- 
tructura se acerque más a los esquemas 


planos de trabajo» 


En general, los propósitos de este tra- 


bajo, no ván más allá de divulgar la exis- 


tencia y características fundamentales de 
estos métodos, y de motivar su estudio y 
aprendizaje a través de las diferentes 
fuentes bibliográficas que se destacan en 
el contenido de la publicación, con vistas 
a que estos conocimientos se reviertan en 
un mejor aprovechamiento de nuestros recur 
sos, en un mayor grado de desarrollo técni- 
co en nuestra rama y, por lo tanto, contri- 
buyan aunque sea en forma discreta al de- 


sarrollo de nuestra economía socialista. 


